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ABSTRACT 


A  method  was  derived  to  determine  the  approximate  hydrograph 
shape  for  ungauged  watersheds  located  in  Alberta  and  southern  Yukon 
Territory.  The  published  discharge  data  for  each  of  70  gauged  watersheds 
was  analysed  by  computer  to  determine  the  most  rapid  rates  of  change  of 
average  daily  flow  at  every  discharge  level  within  the  measured  range. 

The  most  rapid  rates  of  rise  and  fall  were  plotted  as  logarithims  of  dis¬ 
charge  versus  time  in  the  form  of  single-peak,  composite  hydrographs.  The 
hydrograph  limbs  were  approximated  by  straight  lines  and  a  realistic  base- 
.flow  discharge  determined  from  each  graph. 

A  number  of  topographical  and  meteorological  parameters  were 
determined  for  each  of  the  70  watersheds,  and  regression  equations 
established  using,  in  turn,  the  maximum  rates  of  rise  and  fall  and  baseflow 
discharge  as  dependent  variables,  in  a  step-wise  multiple  linear  regression 
analysis . 

An  idealised  hydrograph  was  constructed  for  each  of  the  70 
watersheds,  based  on  the  maximum  rates  of  rise  and  fall  and  baseflow 
discharge  previously  determined,  and  representing  a  particular  volume  of 
direct  runoff.  This  led  to  the  derivation  of  an  equation  relating  peak 
flow  to  the  depth  of  excess  precipitation  over  the  basin  and  to  the  idealised- 
hydrograph  characteristics.  The  equation  was  used  to  compute  an  expected 
peak  flow  for  each  of  the  study  watersheds,  for  depths  of  excess 
precipitation  ranging  from  0.1  to  10  ins.  Regression  equations  were  obtained, 
relating  peak  flow  for  each  depth  of  excess  precipitation  to  the  topographical 
and  meteorological  characteristics  of  the  watersheds. 

The  equations  were  simplified  into  one  expression  relating  peak 
flow  to  depth  of  excess  precipitation  and  drainage  basin  area.  This 
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expression  was  then  refined  to  take  account  of  hydrograph  shape.  The 
expressions  derived  from  the  study  allow  the  determination  of  the 
approximate  shape  of  the  hydrograph  resulting  from  a  given  depth  of 
excess  precipitation  for  an  ungauged  watershed.  A  knowledge  of  certain 
readily  available  basin  characteristics  is  required  and  the  method  is 
applicable  only  to  the  region  from  which  the  sample  was  derived,  and  for 
basins  up  to  15,000  sq .  miles  in  area. 
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CHAPTER  I 


INTRODUCTION 


1 . 1  General 

Before  hydrographs  of  ungauged  watersheds  can  be  determined, 
the  hydrographs  of  gauged  watersheds  must  first  be  quantified,  and 
relationships  established  between  the  parameters  which  describe 
hydrograph  shape  and  other  parameters  which  would  be  available  for 
ungauged  watersheds,  for  example,  the  physical  characteristics  of  the 
basin. 

As  early  as  1930,  attempts  were  made  to  represent  a  flood 
hydrograph  of  a  particular  river  by  one  graph  (referred  to  as  the 
T-hour  unit  hydrograph)  with  excess  precipitation  as  the  only  variable, 
and  it  was  accepted  that  the  effect  of  all  physical  characteristics  of 
the  drainage  basin  were  incorporated  in  this  graph.  Later,  the  ins¬ 
tantaneous  unit  hydrograph  was  developed  in  order  to  eliminate  the 
parameter  T  (the  duration  of  the  uniform  excess  precipitation) . 

The  unit  hydrograph  is  used  to  determine  flood  hydrographs 
resulting  from  severe  storms  through  the  simple  procedure  of  super¬ 
position,  in  spite  of  the  fact  that  the  assumptions  on  which  it  is 
based  are  questionable.  A  limitation  of  the  use  of  the  unit  hydrograph 
is  that  preparation  requires  detailed  knowledge  of  the  duration  of  the 
precipitation  which  caused  the  direct  runoff.  For  many  rivers,  this 
information  is  not  available  and  the  unit  hydrograph  cannot  be  derived. 
Usually,  the  snowmelt  hydrograph  is  different  from  the  rainfall  hydro¬ 
graph. 

There  have  also  been  attempts  by  Nash  (1957) ,  De  Coursey 
(1966)  and  Gupta  (1974)  to  develop  and  test  a  mathematical  equation 
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to  describe  the  shape  of  a  surface  runoff  hydrograph.  The  attempts 
have  achieved  considerable  success,  but  in  order  to  evaluate  the 
parameters  of  the  equation  for  a  particular  basin,  usually,  at  least 
one  representative  hydrograph  from  that  basin  is  necessary.  Thus, 
these  mathematical  representations  of  hydrograph  shape  cannot  be 
derived  for  ungauged  watersheds.  It  would  be  advantageous  if 
characteristic  watershed  hydrographs  could  be  represented  by  one  or 
more  numerical  parameters  to  allow  correlation  with  pertinent  drainage 
basin  characteristics  through  a  regression  analysis,  so  that  the 
parameters  could  be  derived  for  ungauged  watersheds . 

1 . 2  Objectives 

The  objectives  of  the  study  were  as  follows: 

1)  To  write  a  computer  program  to  determine  the  greatest  rates  of 
change  of  discharge  from  a  period  of  hydrometric  record,  and  to  plot 
a  single  peak  hydrograph  composed  of  the  steepest  sections  so  found. 

2)  To  plot  composite  hydrographs  for  a  number  of  selected  gauged 
watersheds  located  in  Alberta  and  southern  Yukon  Territory. 

3)  To  approximate  the  rising  and  falling  limbs  of  each  composite 
hydrograph  by  straight  lines.  (The  slopes  of  such  lines  could  be  used 
as  a  watershed  classification  system.) 

4)  To  compute  regression  equations  relating  the  slopes  of  the  rising 
and  falling  composite  hydrograph  limbs  to  certain  physical  characteris¬ 
tics  of  the  watershed. 

5)  To  prepare  idealized  straight  line  hydrographs  for  each  sample 
watershed,  to  determine  a  relationship  between  the  peak  discharge  and 
the  volume  of  direct  runoff,  and  to  correlate  the  peak  flow  with  the 
physical  characteristics  of  the  watersheds. 
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6)  Using  the  relationships  of  4  and  5}  to  derive  a  method  to  determine 
basic  hydrograph  shape,  for  a  given  depth  of  excess  precipitation,  for 
ungauged  watersheds. 

1 . 3  Factors  Affecting  Hydrograph  Shape 

The  shape  of  a  hydrograph  resulting  from  a  single,  short- 
duration  storm  distributed  uniformly  over  the  drainage  area  follows  a 
general  pattern;  a  rising  portion,  a  crest  segment  and  a  falling 
portion  (Fig.  1) .  The  shape  is  influenced  by  climatic  factors  and  by 
the  topographical  and  geological  features  of  the  basin.  If  the  input 
were  distributed  uniformly  over  the  basin,  then  the  shape  of  the  rising 
limb  would  be  dependent  upon  the  intensity  and  duration  of  the  input. 
However,  it  is  generally  assumed  that  the  recession  limb  would  be 
virtually  independent  of  the  storm  characteristics  since  it  represents 
the  withdrawal  of  water  from  storage  in  the  basin,  after  input  has 
ceased.  The  shape  of  the  recession  limb  is  almost  solely  dependent 
upon  the  physical  features  of  the  basin  (Gray  1970) .  Ideally,  for  a 
particular  watershed,  if  no  further  input  takes  place  after  the  cessation 
of  uniformly  areally-distributed  input,  the  recession  portion  of  the 
hydrograph,  from  any  particular  discharge  value,  will  be  the  same, 
regardless  of  the  shape  of  the  rising  limb.  The  time  base  of  the  rising 
limb  is  definitely  influenced  by  the  duration  of  the  input.  (Linsley  1949). 

For  large  basins  in  particular,  it  would  be  most  unlikely  for 
input  to  be  distributed  uniformly  over  the  basin,  therefore  the 
hydrograph  shape  is  influenced  by  the  areal  distribution.  In  general, 
input  concentrated  near  the  basin  outlet  will  produce  steeper  rising  and 
recession  limbs  and  a  greater  peak  flow  for  a  given  volume  of  excess 
precipitation,  than  if  the  same  volume  were  uniformly,  areally  distributed. 
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The  rising  limb  of  a  typical  hydrograph  is  concave  upwards, 
due  to  the  decreasing  proportion  of  precipitation  'lost'  to  infiltration, 
interception  and  depression  storage  as  the  storm  progresses.  In  addition, 
drainage  basins  tend  to  be  approximately  pear-shaped  so  that  the  area 
between  isocrones  is  greatest  in  the  middle  and  upper  portions  of  the 
basin.  Thus,  uniform  input  over  a  basin  produces  an  increasing  rate 
of  discharge. 

The  recession  limb  is  also  typically  concave  upwards  since  it 
represents  the  exponential  decay  of  the  flood  wave. 

The  area  under  a  hydrograph  between  two  points  in  time,  t^ 
and  (Fig*  1)  represents  the  volume  of  runoff  in  that  time  interval. 

Generally,  such  a  runoff  volume  has  two  components;  direct  runoff  and 
base  flow  runoff.  Direct  runoff  is  basically  water  which  reaches  the 
stream  channel  by  flowing  over,  or  more  usually  through,  the  soil  mantle; 
while  base  flow  is  water  supplied  from  the  saturated  area  below  the 
groundwater  table. 

Base  flow  discharges  remain  relatively  constant  compared  to 
direct  runoff  discharges, for  the  duration  of  direct  runoff,  and  it  is 
useful  to  separate  direct  runoff  from  baseflow.  Many  separation  methods 

are  in  current  use  but  all  involve  a  large  measure  of  subjectivity. 

There  can  be  no  exact  separation  since  the  physical  change  from  direct 
flow  to  baseflow  is  transitional.  However,  a  separation  line  on  a 
hydrograph  can  usually  be  drawn  to  produce  a  sensible  division  into 
direct  runoff  and  baseflow. 

It  is  apparent  from  Fig.  1  that  for  a  given  volume  of  direct 
runoff,  the  greatest  peak  discharge  will  occur  when  the  rising  and  falling 
limbs  are  steepest.  A  steep  rising  limb  will  result  from  intense  input, 
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particularly  if  it  is  concentrated  near  the  basin  outlet,  and  a  steep 
falling  limb  will  be  produced  when  there  is  no  further  input  during 
the  recession  period. 

1 . 4  Maximum  Hydrograph  Slopes 

It  is  generally  accepted  that  there  is  a  limit  to  the  intensity 
and  duration  of  rainfall  which  can  possibly  occur  over  a  particular 
area.  The  limiting  combination  produces  the  Maximum  Probable  Storm. 

The  assumption  that  there  is  a  limit  to  the  combination  of  meteorological 
factors  which  produce  rain,  can  be  extended  to  cover  depth  of  snowpack 
and  melt-rate.  Thus,  we  can  assume  that  there  exists  a  maximum  possible 
rate  of  input  at  any  particular  site.  The  Maximum  Probable  Flood, 
produced  as  a  result  of  these  extreme  input  rates  over  a  drainage  basin, 
theoretically  has  a  return  period  of  infinity  (Koelzer  and  Bitoun  1964) , 
and  therefore  will  not  have  occurred  during  any  period  of  stream  gauge 
records.  However,  it  is  quite  possible  that  maximum  intensities  have 
been  reached  for  short  durations,  or  have  occurred  over  portions  of  basins, 
for  example,  near  the  outlets,  so  that  maximum  rates  of  rise  may  have 
been  achieved,  or  at  least  approached,  for  short  periods  during  different 
flood  events .  By  combining  the  maximum  rates  of  rise  of  portions  of 
various  hydrographs,  the  rate  of  rise  corresponding  to  a  rare  flood  event 
may  be  approximated. 

A  number  of  methods  exist  to  determine  the  steepest  limbs  of 
hydrographs  from  gauged  basins  but  they  are  not  completely  satisfactory. 

The  simplest  method  is  to  take  the  rising  and  falling  limbs 
directly  from  a  single  large  flood  hydrograph,  preferably  one  resulting 
from  an  intense  input.  An  improved  method  is  to  determine  the  slopes 
of  the  limbs  of  several  flood  hydrographs  and  select  the  steepest  one. 


7. 


The  limitations  of  such  a  method  are:  - 

i)  that  it  is  unlikely  that  input  conditions  remained 
constant  for  the  duration  of  the  rising  limb,  and 
ii)  it  cannot  be  known  whether  the  slopes  are  the  steepest 
possible . 

An  improved  method  is  to  build  up  a  hydrograph  covering  the 
entire  discharge  range  of  the  record,  by  selecting  the  steepest  slope 
which  has  occurred  at  every  discharge  value.  The  resulting  hydrograph 
limbs  are  thus  a  synthesis  of  parts  of  different  flood  hydrographs,  and 
consist  of  those  portions  occurring  when  conditions  producing  rapid  rates 
of  rise  and  fall  were  most  extreme.  A  problem  with  this  method  is  that 
to  obtain  the  slopes  accurately,  the  hydrograph  must  be  drawn  to  a  large 
time  scale  resulting  in  extremely  long  and  cumbersome  plots. 

A  statistical  analysis  of  the  falling  limb  values  (Verschuren 
1970)  was  presented  in  the  form  of  plots  of  discharge  on  one  day  vs. 
discharge  on  the  following  day  Qn  +  Data  for  these  plots  was  obtained 
where  the  hydrograph  crossed  preselected  levels  of  discharge  Q^,  such 
that  .  For  the  data  at  each  discharge  level,  the  gradient 

of  the  graph  was  computed  and  the  and  Qn+^  values  corresponding  to 
the  median,  the  tenth  percentile  and  the  nintieth  percentile  of  the 
gradient  were  plotted. 

It  was  recognized  that  the  computer  analysis  was  not  satisfactory 
in  many  cases  because  some  detail  was  missed  by  using  preselected  discharge 
levels,  and  more  important,  the  output  did  not  allow  a  subjective  look 
at  the  hydrograph  shape. 

The  first  objective  of  this  study  was  to  write  a  computer 
program  that  would  calculate  the  co-ordinates  of  a  hydrograph,  composed 
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of  the  steepest  rising  and  falling  portions  of  a  discharge  record  at  a 
station,  and  which  covered  the  entire  discharge  range.  The  hydrograph 
could  then  be  plotted  to  any  desired  scale  so  that  realistic  slopes  of 
the  limbs  could  be  assessed.  A  description  of  the  method  used  to  obtain 
such  composite  hydrographs  is  given  in  Chapter  2. 

1 . 5  Watersheds  Included  in  Study 

In  Alberta,  as  in  most  other  areas,  a  disproportionate  number  of 
large  watersheds  are  permanently  gauged,  resulting  in  a  lack  of  information 
on  smaller  watersheds  (less  than  200  sq.  miles  in  area) .  A  knowledge  of 
expected  peak  discharges  on  small  watersheds  is  essential  for  the  econo¬ 
mical  design  of  an  increasing  number  of  stream  crossings  and  other 
fluvial  structures.  The  absence  of  direct  streamflow  measurements  could 
be  largely  overcome  if  hydrograph  parameters  from  gauged  watersheds 
could  be  correlated  with  their  respective  basin  physical  characteristics. 
The  resulting  relationships  could  then  be  used  to  predict  hydrograph 
parameters  from  a  knowledge  of  the  basin  characteristics  alone. 

Discharge  data  published  by  Water  Survey  of  Canada  is  in  the 
form  of  mean-of-the-day  values.  These  are  evaluated  by  determining  the 
average  gauge  height  for  each  day  and  selecting,  from  the  station  rating 
curve,  the  discharge  associated  with  that  gauge  height.  For  a  large 
watershed,  where  the  rates  of  change  of  discharge  are  small,  the  hydro¬ 
graph  prepared  from  mean-of-the-day  values  is  close  to  that  prepared  from 
instantaneous  discharge  values.  For  a  small  watershed,  however  where 
discharges  are  rapidly  changing,  the  hydrograph  resulting  from  mean-of- 
the-day  values  differs  considerably  from  the  hydrograph  of  instantaneous 
values,  particularly  in  the  region  of  discharge  peaks  and  troughs.  Thus, 
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completely  erroneous  hydrograph  parameters  can  result  from  analysis  of 
mean-of-the-day  discharge  hydrographs  for  small  watersheds. 

Since  hydrographs  derived  from  instantaneous  data  are  not 
readily  available,  the  present  study  is  restricted  to  watersheds  greater 
than  500  sq.  miles.  This  figure  was  selected  after  comparisons  were 
made  of  instantaneous  and  mean-of-the-day  discharge  hydrographs,  for 
the  same  flood  events,  on  basins  ranging  from  1  to  2,000  sq .  miles.  For 
basins  of  500  sq.  miles  and  greater,  the  slopes  of  rising  and  falling 
limbs  of  the  mean-of-the-day  discharge  hydrograph  closely  approximated 
those  of  the  instantaneous  hydrographs. 

As  stated  in  Section  1.2,  one  objective  of  this  study  was  to 
correlate  parameters  describing  hydrograph  shape  with  drainage  basin 
characteristics.  While  appreciating  that  no  drainage  basin,  however 
small,  is  entirely  homogeneous  in  terms  of  slope,  aspect,  permeability, 
etc. ,  it  is  apparent  that  the  degree  of  homogeneity  decreases  with 
increasing  basin  area.  It  was  considered  desirable  to  use  basins  which, 
as  far  as  possible,  were  located  within  one  vegetation  zone,  e.g. , 
prairie  or  boreal  forest,  so  that  such  characteristics  as  Basin  Latitude 
and  60  Minute  Rainfall  (see  Section  4.1)  were  reasonably  representative. 
An  upper  time  limit  to  the  size  of  basins  included  in  the  study  was  set 
at  15,000  sq.  miles,  which  is  a  compromise  between  two  considerations:  - 
i)  to  keep  the  basins  small  and  therefore  homogeneous 
ii)  to  include  in  the  regression  analysis,  as  many  basins 
greater  than  500  sq.  miles,  as  possible. 

Sixty  two  watersheds  were  included  in  the  analysis  and  each 
satisfied  the  following  requirements:  - 
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i)  Discharge  from  the  watershed  had  been  measured  for  at 
least  10  years,  either  continuously  or  seasonally.  (In  areas  where 
gauged  watersheds  were  sparse,  7  years  of  record  was  accepted).* 

ii)  The  watershed  area  was  between  500  and  15,000  sq.  miles, 
and  iii)  The  location  was  such  that  it  gave  reasonably  even 
distribution  of  gauging  stations  throughout  the  region  of  study. 

In  regions  that  required  elimination  of  some  stations  in 
accordance  with  iii)  above,  preference  was  given  to  a  watershed  not 
situated  within  a  larger  one,  already  selected  for  study. 

Forty-two  watersheds  in  Alberta  met  the  above  requirements 
and  were  included  in  the  study.  The  scope  and  accuracy  of  the  analysis 
was  enhanced  by  the  inclusion  of  20  watersheds  (which  also  fulfilled 
the  above  requirements)  situated  in  the  Yukon  Territory  and  for  which 
the  necessary  data  existed  (Spence  1969).  Thus,  a  total  sample  of 
62  watersheds,  for  which  reasonable  periods  of  discharge  records 
existed,  were  scrutinized  and  included  in  the  regression  analysis. 

(For  station  locations,  see  Appendix  IV). 


*  The  last  year  of  record  available  for  inclusion  in  the  study  was  1972. 


CHAPTER  2 


COMPOSITE  HYDROGRAPHS 


2 .1  Method  of  Preparing  Composite  Hydrographs 

A  package  of  programs  was  written  to  produce  the  steepest 
hydrograph  possible  from  the  available  hydrometric  records  of  a 
particular  gauging  station,  and  plot  the  result  on  semilogarithmic 
co-ordinates.  The  whole  of  the  hydrometric  record  may  be  used  as 
input  or  it  can  be  divided  in  3  ways:  - 

i)  chronologically:  N  consecutive  years  of  data  may 
be  used,  where  N  can  range  from  1  to  the  total 
number  of  years  of  record, 
ii)  by  months:  n  consecutive  months  of  data  from 

each  year  may  be  included,  where  n  can  range  from 
1  to  12 , 

or  iii)  a  combination  of  i  and  ii,  i.e.,  part  of  the  year 

for  part  of  the  period  of  record. 

The  first  case  is  useful  where  the  period  of  record  is  long 
and  computer  costs  are  to  be  minimized.  The  second  case  is  useful 
when  rates  of  rise  and  fall  are  known  to  be  small,  i.e.,  during  winter 
months.  In  such  a  case  the  winter  data  can  be  omitted  from  the 
computations,  again  resulting  in  a  considerable  saving  in  computer 
costs . 

The  central  program  in  the  package  is  "SLOPE"  which  reads 
discharge  values  from  magnetic  tape  (such  as  published  by  Water  Survey 
of  Canada) .  It  is  essential  that  the  discharge  data  are  equally  spaced 
in  time,  for  example  daily  or  hourly.  On  an  input  file  "Data",  the 
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operator  specifies  the  length  of  record  to  be  used,  and  the  first  and 
last  months  of  each  year  (n^  and  n^  for  which  data  is  to  be  included. 

Missing  discharge  values  and  discontinuities  in  the  record 
are  handled  by  the  program  so  that  no  visual  inspection  of  the  input 
data  is  required. 

Since  the  CPU  time  used  by  the  program  "SLOPE"  (and  hence 
the  cost)  is  nearly  proportional  to  the  square  of  the  number  of  discharge 
values  included,  it  has  been  found  economical  to  divide  periods  of 
record  of  greater  than,  about  7  years  into  3  or  4  consecutive  periods. 

The  program  "SLOPE"  is  called  for  each  successive  sub-period  and 
produces  the  steepest  hydrograph  for  that  sub-period. 

The  3  or  4  resulting  hydrographs  are  combined  by  program  "ADD", 
the  output  of  which  becomes  the  input  for  program  "SLOPE"again ,  which 
then  produces  the  required  hydrograph  co-ordinates  for  the  whole  period. 

The  co-ordinates  of  the  final  composite  hydrograph  are  read  by 
program  "DRAW"  which  selects  the  logarithmic  ordinate  scale  range, 
curtails  the  abscissor  values,  if  necessary,  to  ensure  that  the  plot 
length  does  not  exceed  36  ins.  and  restricts  the  discharge  axis  length  to 
5  log-cycles.  A  plot  is  produced  of  a  single  hydrograph  at  an  ordinate 
scale  of  5.2  ins.  per  log-cycle  and  an  arithmetic  time  scale  of  1  ins. 
per  day.  The  above  ordinate  scale  was  chosen  so  that  the  maximum  dis¬ 
charge  range  which  occurred  on  any  one  basin  (5  log-cycles)  just  covered 
the  maximum  possible  plot  height  of  26  ins. 

2.2  Program  "SLOPE" 

The  complete  listing  for  program  "SLOPE"  is  given  in  Appendix 

I.  Discharge  values  on  successive  days  are  read  by  the  program.  If  the 
discharge  for  a  particular  day  is  missing,  the  value  is  replaced  by 
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-999.0.  (Water  Survey  of  Canada  magnetic  tapes  are  compiled  using  this 
system).  To  illustrate  the  method,  the  23-day  discharge  record  for 
which  the  steepest  composite  hydrograph  is  required,  is  plotted  on  semi- 
logarithmic  axes  (Fig.  2(a)). (In  practice,  the  input  data  would  not  be 
plotted).  It  can  be  seen  that  discharge  values  for  days  4,  5,  7  and 
17  are  missing  (i.e.,  -999.0). 

The  required  composite  hydrograph  (Fig.  2(b))  is  composed  of 
the  heavy  lines  in  Fig.  2(a)  and  was  prepared  as  follows:  - 

The  discharge  values  for  the  specified  period  are  read  into  an 
array  and  their  common  logarithms  calculated  (from  this  point  on,  the 
program  deals  entirely  with  logarithms  of  discharge  values) .  These 
values  are  arranged  in  descending  order (points  9,  11,  12,  13,  10,  etc.). 
Alongside  each  value,  the  discharge  value  occurring  on  the  following 
day  (when  it  exists)  is  listed  (points  10,  12,  13,  14,  11,  etc.,  res¬ 
pectively)  .  This  2- dimensional  array  is  separated  into  two  further 
arrays : - 

i)  QP,  containing  pairs  of  values  with  a  positive  gradient 
e.g.,  point  10  -  point  11,  and 
ii)  QPAIR,  containing  pairs  of  values  with  a  negative 
gradient  e.g.,  point  9  -  point  10  and  point  11  - 
point  12. 

A  search  procedure  is  then  carried  out,  separately  for  the 
rising  and  falling  limbs,  commencing  in  each  case  with  the  maximum 
discharge  value  (point  9),  and  proceeding  in  descending  order.  For  each 
successive  value,  those  pairs  of  values  from  QP  or  QPAIR  which  span 
the  discharge  value  under  consideration  are  listed  (e.g.,  for  the  rising 
limb,  point  13,  is  spanned  by  lines  10  -  11  and  8-9).  From  each 
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FIG.  2  ILLUSTRATION  OF  METHOD  OF  HYDROGRAPH  COMPOSITION 
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list  the  steepest  gradient  (line  8-9  for  point  13)  is  selected  and 
written  into  a  further  array  together  with  the  current  spanned  value 
(point  13) . 

A  further  gradient-comparing  procedure  results  in  2  arrays; 
one  for  the  rising  limb,  consisting  of  discharge  values  in  ascending 
order,  together  with  the  steepest  positive  gradient  preceding  each 
value;  and  for  the  falling  limb,  an  array  consisting  of  discharge 
values  in  descending  order  together  with  the  steepest  negative  gradient 
following  each  value.  These  arrays  contain  the  necessary  data  to 
construct  the  rising  and  falling  limbs  of  the  composite  hydrograph,  as 
follows:  - 

The  smallest  discharge  value  on  the  rising  limb,  for  which 
there  exists  a  preceding  linej is  point  2.  This  value  is  given  a  time 
co-ordinate  of  t^=1.0.  The  time  co-ordinate  for  the  second  smallest 
discharge  value  on  the  rising  limb  (point  3)  is  determined  by  proceeding 
from  point  2,  at  the  gradient  preceding  point  3,  until  the  value  of 
point  3  is  reached.  The  time  interval  so  covered  is  added  to  t^  to 
obtain  t2«  The  process  is  repeated  until  the  maximum  discharge  is 
reached.  The  discharge  value  at  t  =0.0  is  obtained  by  descending,  for 
one  day,  at  the  gradient  preceding  point  2. 

A  similar  process  to  that  described  above  is  used  to  compute 
the  co-ordinates  of  the  falling  limb,  starting  from  the  peak  value.  The 
result  is  a  list  of  co-ordinates  composed  of  the  steepest  parts  of  the 
hydrograph  which  have  occurred  during  the  period  of  record.  Finally, 
antilogarithms  are  taken  of  the  discharge  values  in  readiness  for  their 
being  plotted  on  a  logarithmic-scale,  ordinate  axis. 
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It  can  be  seen  from  Fig.  2(a)  that  only  part  of  the  line 
18-  19,  for  example,  is  heavy,  and  as  such  only  the  upper  part  appears 
in  the  hydrograph  of  Fig.  2(b)  .  This  is  because,  below  the  discharge 
value  associated  with  point  22,  the  line  22  -  23  is  steeper  than  line 
18  -  19  and  therefore  appears  in  the  final  graph. 

A  situation  may  arise  where  the  discharge  record  does  not 
cover  a  certain  discharge  range,  for  example,  in  Fig.  2(a),  the  dis¬ 
charge  range  12  to  14  cfs  is  not  covered  by  the  rising  limb  of  the 
input  record.  In  this  situation,  in  order  to  simplify  plotting  of  the 
output  hydrograph,  the  line  immediately  above  the  14  cfs.  level 
(preceding  point  20)  is  continued  at  the  same  slope  across  the  gap,  to 
12  cfs.  (point  3)  as  can  be  seen  in  Fig.  2(b). 

Because  of  the  above  feature  of  the  program,  when  a  period  of 
record  is  to  be  split  by  months,  care  should  be  taken  to  ensure  that  the 
discharge  range  covered  by  each  limb  is  likely  to  be  continuous .  A 
typical  annual  hydrograph  for  a  river  with  a  peak  in  June  should  not, 
if  possible,  be  divided  into  the  two  periods  January  to  May  and  June  to 
December.  The  hydrograph  for  the  first  half-year  is  predominantly 
rising  while  that  for  the  second  half  is  predominantly  falling.  This 
could  lead  to  a  discontinuous  input  record  and  resulting  artificial 
slopes  on  the  output  hydrograph.  A  better  division  would  be  April  to 
September  and  October  to  March. 

2.3  Program  "ADD" 

A  considerable  saving  in  computer  costs  can  result  by  dividing 
a  long  period  of  record  into  3  or  4  (or  5)  consecutive  sub-periods.  In 
such  cases,  calling  program  "SLOPE"  for  each  sub-period  results  in  the 
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output  of  the  co-ordinates  of  the  steepest  hydrograph  for  each  sub-period. 
These  co-ordinates  are  then  required  to  be  read  by  program  "SLOPE”  to 
produce  one  final  hydrograph.  However,  except  for  an  extremely  trivial 
case,  the  co-ordinates  of  the  sub-period  hydrographs  will  not  be  equally 
spaced  in  time,  so  that  they  cannot  be  directly  read  as  input  by  program 
"SLOPE".  This  problem  is  overcome  by  calling  program  "ADD",  which 
interpolates  between  the  existing  hydrograph  co-ordinates,  to  produce 
a  new  set  of  discharge  values  spaced  at  half-day  intervals.  Some  accuracy 
is  lost  by  this  process  particularly  near  the  hydrograph  peak,  but  since 
the  rising  and  falling  limbs  are  reasonably  straight,  the  interpolated 
co-ordinates  adequately  describe  the  original  curve. 

Program  "ADD"  operates  on  each  sub-period  in  turn  and  joins 
the  interpolated  co-ordinates  into  a  continuous  record  which  can  be  read 
as  input  by  program  "SLOPE",  to  produce  the  co-ordinates  for  the  single, 
final  hydrograph  for  the  entire  period.  When  program  "SLOPE"  is  required 
to  read  discharge  values  produced  by  program  "ADD",  the  program  "SLOPE" 
input  parameter  ADD,  must  be  set  equal  to  1.0  (otherwise  ADD=0.0). 

2 . 4  Hydrograph  Parameters 

The  foregoing  method  was  used  to  produce  a  composite  hydrograph 
for  each  of  the  62  watersheds  included  in  the  study.  The  entire  period 
of  record  for  each  gauging  station  was  fed  into  the  computer  which  plotted 
the  required  composite  hydrographs  to  ordinate  scales  of  5.2  ins/log- 
cycle  and  abscissor  scales  of  1  ins/day.  All  hydrographs  were  plotted 
to  the  same  scale  to  permit  easy  visual  comparison. 

One  of  the  objectives  of  the  study  was  to  find  the  regression 
equations  between  the  watershed  characteristics  as  independent  variables 
and  the  maximum  rate  of  rise,  maximum  rate  of  fall  and  base  flow  dis- 
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charge,  in  turn,  as  the  dependent  variables;  the  dependent  variables 
being  determined  from  the  computer-plotted  composite  hydrographs  for  each 
river  gauging  station  included  in  the  study. 

This  study  is  an  attempt  to  demonstrate  that  maximum  rates  of 
rise  and  fall  of  hydrographs  can  sometimes  be  approximated  by  an  equation 
of  the  form 


where  Q  =  discharge  at  time  t 

Qq=  discharge  at  time  t=0 
and  K  =  a  constant 

If  such  a  relationship  exists,  then  a  plot  of  time  versus  the 
logarithm  of  discharge  (or  time  vs  discharge  plotted  on  semi-logarithmic 
paper)  results  in  a  straight  line,  whose  slope  can  be  described  by  the 
constant  K  (For  fuller  analysis,  see  Section  5.1). 

An  inspection  of  the  62  composite  hydrographs  showed  that 

the  rising  and  falling  limbs  could,  in  general,  be  approximated  by  an 

t 

equation  of  the  form  Q  =  Q^K  ,  since  the  rising  and  falling  limbs  were 
often  straight  over  a  range  of  discharges  of  several  log  cycles,  and  in 
some  instances  were  practically  straight  over  the  entire  discharge 
range.  In  many  cases,  where  the  lines  were  not  straight,  there  were 
sections  of  the  line  distributed  throughout  the  discharge  range  which 
were  virtually  at  the  same,  maximum  slope.  It  will  be  recalled  that  the 
composite  hydrographs  do  not  represent  one  flood  event,  but  rather,  are 
comprised  of  sections  of  several  flood  events  which  had  occurred  during 
the  (often  limited)  period  of  record.  If  a  limiting  slope  were  to  exist, 
then  it  would  be  expected  that,  as  a  longer  period  of  record  became 
available,  more  sections  of  the  line  would  tend  towards  that  maximum, 
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limiting  slope. 

A  factor  which  could  contribute  to  non-linearity  of  the  rising 
and  falling  limbs,  is  the  fact  that  in  the  lower  discharge  range,  floods 
of  that  magnitude  are  more  numerous,  providing  a  greater  sample  from 
which  to  choose  the  steepest  slope.  If  a  limiting  slope  exists,  it 
will  probably  occur  in  the  lower  and  intermediate  discharge  ranges. 
Unless  the  observed  maximum  discharge  rate  occurred  more  than  once 
during  the  period  of  record,  the  peak  of  the  composite  hydrograph  will 
be  an  exact  replica  of  the  peak  of  that  maximum  flood.  (This  was  shown 
to  be  the  case  in  many  instances) . 

Two  measures  of  slope,  for  each  limb,  were  determined  from 
the  composite  hydrograph  of  each  study  watershed: 

i)  R'Cmax)  =  the  steepest  slope  which  occurred  for  a  signi¬ 
ficant  discharge  range  (e.g.,  more  than  1 
log-cycle) 

and  ii)  R'(av)  =  an  average  slope  between  R'(max)  and  a  steep 

slope  which  recurs  in  several  discharge  ranges. 
Fig.  3  shows  a  typical  rising  limb  of  a  composite  hydrograph  and  illus¬ 
trates  the  method.  Section  CD  has  the  steepest  slope  over  a  significant 
discharge  range  and  hence  its  slope  is  taken  as  Rv(max)  .  Sections  AB , 

DE  and  FG  have  approximately  the  same  slope  so  that  R*(av)  is  taken  as  an 
average  between  R'(max)  and  the  slope  of  section  DE.  It  is  assumed  that 
if  the  period  of  record  were  longer,  the  slopes  of  sections  BC  and  EF 
would  approximate  R'(av)  .  The  method  is  highly  subjective  and  many 
possible  values  of  R’ could  be  selected.  However,  the  above  procedure 
allowed  the  selection  of  "an  average  steepest  slope"  with  a  reasonable 
degree  of  consistency. 
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FIG.  3 


TYPICAL  RISING  LIMB  OF  COMPOSITE  HYDROGRAPH 
SHOWING  TWO  MEASURES  OF  ITS  STEEPEST  SLOPE 
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For  the  falling  limb,  the  two  slopes  K^max)  and  K'(av)  were 
selected  in  a  similar  manner. 

An  unexpected  feature  exhibited  by  many  of  the  composite 
hydrographs  was  an  initial,  almost  horizontal  portion  of  the  rising 
limb  followed  by  an  abrupt  upward  turn  to  the  approximately  straight 
line  portion.  The  recession  limb  exhibited  a  similar  trend  but  the 
change  from  the  general  recession  slope  to  the  final  horizontal 
section  was  more  gradual.  However,  the  discharge  values  of  the  two 
horizontal  portions  were  usually  very  close.  This  feature  indicated 
that  there  was  a  real  distinction  between  baseflow  and  direct  runoff. 

In  other  words,  that  without  further  precipitation  input,  the  discharge 
from  any  particular  level  would  decline  quite  rapidly  to  an  almost 
constant  value  which  presumably  would  be  base  flow.  The  discharge 
associated  with  the  change  from  the  horizontal  to  the  steeply  sloping 
section  of  the  graph  was  used  as  the  final  dependent  variable  in  the 
regression  analysis  (BF) . 

For  those  stations  where  the  minimum  recorded  discharge  was 
zero  (i.e.,  those  rivers  that  intermittently  dry  up),  there  was  no 
transition  between  zero  flow  and  the  straight  line  portion:  the 
straight  line  continued  at  the  same  slope  right  down  to  zero  -  even 
steepening  on  occasions  as  it  approached. 

A  list  of  the  hydrograph  slopes  and  baseflow  values  is  given 
in  Appendix  II,  and  a  description  of  the  regression  analyses,  with  the 
watershed  characteristics  as  independent  variables,  is  given  in  Section 
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CHAPTER  3 


THEORY  OF  CORRELATION  &  REGRESSION 

A  brief  overview  of  linear  regression  and  correlation  is 
given  together  with  definitions  of  some  of  the  more  important  parameters 
calculated  in  the  following  sections:  - 

3 . 1  Correlation  Between  Two  Variables 

t  Consider  a  set  of  observations  on  two  variables,  X  and  Y, 

then  the  simple  correlation  coefficient  (r  )  is  defined  as  the  covariance 

xy 

divided  by  the  product  of  the  standard  deviations  of  the  two  variables 

(-1<  r  <+l) .  The  simple  correlation  coefficient  is  a  measure  of  the 
xy 

linear  association  between  the  two  variables.  If  X  can  be  expressed  as 

an  exact  linear  function  of  Y,  then  r  =il .  The  closer  the  absolute 

xy 

value  of  rXyj  approaches  zero,  the  further  the  observations  depart  from 
a  linear  relationship . 

Consider  now  a  set  of  observations  on  n  variables  X.,  X„  .  .  .X  . 

12  n 

A  high  simple  correlation  coefficient  (r^)  maY  exist  between  X^  and 
X^  as  defined  above,  which  may  be  because:  - 

i)  an  increase  or  decrease  in  X^  causes  an 
increase  in  X^ 

or  ii)  an  increase  in  X^  and  X^  is  caused  by  an 

increase  in  one  or  more  of  X....X  . 

3  n 

Thus,  r^  may  give  a  false  impression  of  the  linear  association 
between  X^  and  X^ .  Therefore,  it  is  of  interest  to  be  able  to  calculate 
the  correlation  coefficient  between  two  variables  with  the  effect  of  one 
or  more  other  variables  removed.  Such  a  correlation  coefficient  is 


known  as  the  partial  correlation  coefficient  and  is  written  r 


12 .34  .  .  ,n 
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and  measures  the  linear  association  between  and  with  the  effect 

of  variables  X-...X  removed.  This  is  an  (n-2)nd  order  coefficient. 

3  n 


3 . 2  Linear  Regression 

The  subject  of  linear  regression  is  concerned  with  the 
dependence  of  one  variable  (the  dependent  variable)  on  a  linear  combi¬ 
nation  of  a  set  of  other  variables  (the  independent  variables) . 

It  is  assumed  that  there  is  a  linear  relationship  of  the  form 


Y=Bq  +  Bx  X  +  B2  X2  + 


.+  B  X  +  u 
P  P 


between  a  dependent  variable  Y  and  the  independent  variables  X^...X^. 

The  parameters  Bq,  B^....B  are  unknown  population  regression  coefficients 

and  u  is  an  unknown  random  variable,  which  is  a  measure  of  the  departure 

of  Y  from  an  exact  linear  relationship  with  the  p  independent  variables. 

Certain  assumptions  regarding  the  distribution  of  the  random  variable  u 

need  to  be  made  in  order  to  apply  statistical  tests  to  the  various 

parameters,  but  they  will  not  be  discussed  here  (K.W.  Smillie,  1965). 

Consider  n  sets  of  observations  of  the  random  variables  Y, 

X  ,  X  . . . .X  and  that  b  ,  b.  .  .  .  .b  are  estimates  of  the  population 
12p  o  Id 


regression  coefficients  Bq,  B^....B  . 

Then  Y.  =  b  +  b.  X.,  +  b„  X._ _ b  X.  +  e. 

l  o  1  ll  2  i2  p  ip  l, 


Eq.  3.1 


th 


where  e.  is  an  estimate  of  u.  and  is  called  the  i  residual.  An  estimate 
l  a 

of  Y^  can  be  made  using  the  regression  equation,  thus 


Y’  =  b  +  bnX. 1  +b0X.0  +  _ +b  X. 

l  o  1  ll  2  i2  p  ip 

th 


Eq.  3.2 


It  will  be  seen  that  the  i  residual  e.  is  the  difference  between  the 

i 

observed  and  estimated  value  of  Y. . 

l 

The  values  of  the  regression  coefficients  b  ,  . . . .b  can  be 

°  o  p 

determined  by  minimizing  the  sum  of  the  squares  of  residuals,  which 


24. 


involves  the  solution  of  a  set  of  simultaneous  equations.  When  the 
regression  coefficients  have  been  determined,  values  of  Y  may  be 
estimated  via  the  regression  equation  and  the  n  residuals  computed. 

It  can  be  shown  that 

Z  (Y  -Y) 2  =Z(Y^  -  Y’)2  +Z(Yi-Yp2  Eq  .  3.3 

Where  Z  =  the  sum  from  i  =  1  to  n 

Y  =  the  mean  of  the  observed  Y. 


Y^  =  the  estimated  value  of  Y_^  using  the  regression  equation 


t  __ 


=  the  mean  of  the  estimated  values. 

Equation  3.3  shows  that  the  total  sum  of  squares  of  deviations 

from  the  mean  of  the  dependent  variable,  has  two  components:  - 

_  2 

i)  Z  (Y^  -  Y’)  =  total  sum  of  squares  of  deviation 


and 


from  the  mean  of  the  estimated  values  of  Y 

2 

ii)  Z  (Y  -  Y_I)  =  the  sum  of  squares  of  deviation  of 
the  estimated  value  from  the  observed  value. 


Or  the  total  variation  in  the  dependent  variable 

=  variation  due  to  regression  +  residual  variation 
=  explained  variation  +  unexplained  variation 


3 . 3  Measures  of  Multiple  Linear  Correlation 

There  are  several  parameters  which  measure  the  degree  of 
multiple  linear  correlation  between  variables,  one  of  which  is  the 
multiple  correlation  coefficient  (R)  defined  as  the  simple  correlation 
coefficient  between  the  observed  values  of  the  independent  variables 
and  the  values  estimated  from  the  multiple  regression  function.  If  the 
observed  and  estimated  values  are  equal,  then  the  multiple  correlation 
coefficient  will  be  1. 
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It  can  be  shown  that 

YL  (Y  '  -  Y ' )  2 

^2  _  explained  variation  _  i _ 

total  variation  ZI  (Y  Y)^ 


Eq.  3.4 


Another  measure  of  multiple  linear  correlation  is  the  residual 
2 

mean  square  (S  )  which  measures  the  closeness  with  which  the  estimated 

values  approach  the  observed  values.  To  obtain  the  unbiased  estimate 
2 

of  S  ,  the  residual  mean  square  must  be  adjusted  for  the  number  of 
degrees  of  freedom,  thus 


S2  =  — —  7"  .  Z(Y.  -  Y!)2  Eq.  3.5 

n-p-1  ii 

where  p  is  the  number  of  independent  variables  in  the  regression  equation. 

A  third  measure  of  multiple  linear  correlation  is  the  standard 

error  of  estimate  and  is  defined  as  the  square  root  of  the  residual 

mean  square  as  defined  above. 

As  stated  previously,  the  sample  regression  coefficients 

(bQ,  b^,...b  )  are  merely  estimates  of  the  population  coefficients 

(B  ,  B1 , . . . .B  )  and  are  therefore  random  variables  themselves, 
o  1  p 

It  may  be  shown  that  any  multiple  regression  coefficient  b  is  normally 

distributed  with  a  mean  of  B.  and  that  the  unbiased  standard  error  of 

J 

estimate  of  b .  is 
3 


(n-p-1)  .  (1-R.  ) 


Eq.  3.6 


where  S  =  the  standard  error  of  estimate  (biased) 

s.  =  the  standard  deviation  of  the  marginal  distribution  of 
3 

variable  X. 

3 

and  R.  =  multiple  correlation  coefficient  of  X.  with  respect 

3  3 

to  all  other  independent  variables. 
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3 . 4  The  "Independent"  Variables 

From  the  definition  of  the  multiple  correlation  coefficient, 

we  know  that  if  the  so-called  independent  variables  have  no  linear 

association  amongst  themselves,  then  R  =  0  and  the  standard  deviation 

of  any  regression  coefficient  b  will  be  a  minimum  (Eq.  3.6).  As  the 

degree  of  linear  association  amongst  the  "independent"  variables 

increases,  so  K.  increases  and  .  increases,  so  that  the  regression 

coefficients  become  less  reliable.  In  other  words,  multiple  linear 

regression  analysis  is  valid  no  matter  how  much  the  variables  are  inter- 

correlated,  but  significance  tests  should  be  continuously  applied  to 

the  regression  coefficients  to  avoid  their  standard  deviations  becoming 

so  large  that  the  coefficients  become  worthless. 

To  prevent  the  likelihood  of  the  above  pitfall,  the  independent 

variables  used  may  be  screened,  by  calculating  a  (p  x  p)  matrix  of  the 

(p-2)nd  order  partial  correlation  coefficients  i.e.,  r  etc., 

J.Z  ♦  *  •  •  •  p 

and  eliminating  those  independent  variables  which  have  high  partial 
correlation  coefficients,  (>0.5  for  example). 

The  (p-2)nd  order  partial  correlation  coefficients  may  be 
easily  calculated  from  the  formula 


ij.1,2, .  i-1,  i+1. 


j-l,  j+1. 


-*c.  . 

/* C. . .*c. . 
V  n  33 


Eq.  3.7 


where  *c„  i-s  t*ie  ijth  element  in  the  inverse  matrix  of  simple  (zero- 
order)  correlation  coefficients. 

These  partial  correlation  coefficients  provide  a  true 
measure  of  the  intercorrelation  between  the  independent  variables  and 
should  be  used  in  preference  to  the  simple  correlation  coefficients 
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for  selecting  variables  to  include  in  the  analysis. 

3 . 5  Step  Wise  Multiple  Linear  Regression 

A  method  exists  which  introduces  the  variables  into  the 
regression  equation  one  by  one,  according  to  their  contribution  to  the 
reduction  of  the  sum  of  squares  of  deviations  of  the  dependent  variables. 
(K.W.  Smillie,  1965).  This  allows  significance  tests  to  be  applied  to 
the  regression  coefficients,  and  a  check  to  be  made  on  the  unbiased 
standard  error  of  the  estimate  at  each  step,  so  that  the  introduction 
of  new  variables  can  be  stopped  when  the  unbiased  standard  error  becomes 
a  minimum  or  the  regression  coefficient  of  the  last  variable  introduced 
becomes  non-significant. 

The  selection  of  the  variables  to  be  introduced  proceeds  as 

follows:  - 

The  squares  of  the  simple  correlation  coefficients  between  the 

2 

dependent  and  each  of  the  independent  variables  (r  )  are  calculated. 

The  first  variable  to  be  introduced  into  the  regression  equation 
Y  =  b  b-^^X.  is  the  one  with  the  largest  correlation  coefficient 

o  1  3 

squared,  i.e.,  the  one  that  explains  the  largest  portion  of  the  variation 

of  the  dependent  variable,  X  for  example. 

Now,  the  squares  of  the  first-order  partial  correlation 

coefficients  r  .  .  for  all  i  ^  A,  give  the  proportion  of  the  variation 
yj.A, 

of  Y,  unaccounted  for  by  X  ,  which  may  be  explained  by  each  of  the 

remaining  independent  variables.  The  variable  which  is  introduced  next 

into  the  regression  function  is  that  one,  X^  for  example,  which  accounts 

for  the  largest  portion  of  the  remaining  variation.  The  regression 

(2)  (2)  (2) 

equation  then  becomes  Y  =  b^  '  +  b^  ^A  +  ^2 

Subsequently,  variables  are  introduced  into  the  equation  on 
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the  basis  of  the  largest  square  of  the  partial  correlation  coefficient 
of  Y  on  the  independent  variables,  with  the  effects  of  the  variables 
already  in  the  equation  (X^,  XL^,  etc.)  partialled  out. 

The  basic  criteria  used  to  determine  which  variable  is 
introduced  into  the  equation,  is  that  the  variable  should  contribute 
most  to  the  reduction  of  "unexplained  variation"  at_  each  step.  This 
does  not  imply  that  successive  steps  reduce  the  unexplained  variation 
by  decreasing  amounts.  For  example, 

Let  V_^  =  portion  of  variance  explained  at  step  i 

2 

then  for  i  =  1,  V.  =  r  . 

1  yA 

and  for  i  =  2,  V  =  (1  -  r2  )  .  r2 

2  yA  yB.A 

Let  r2  =  0.2  >  r2 
yA  yB 


2 

and  let  r  _  .  =0.3 

yB  .A 

2 

then  V.  =  r  .  =  0.2 

1  yA 

and  V2  =  (1  -  0.2)  x  0.3  =  0.24 

i.e  VV1 

Thus  the  portion  of  the  "unexplained  variation"  reduced  in  the  first 
step  is  0.2,  and  the  portion  reduced  in  the  second  step  is  0.24. 

If  the  independent  variables  are  highly  correlated  amongst 
themselves,  then  their  rank,  on  the  basis  of  each  successive  order  of 
partial  correlation  coefficients,  can  change  drastically  from  step  to 
step,  as  one  more  variable  is  partialled  out. 

It  is  possible  that  a  variable  may  rank  high  according  to  its 
(p-2)nd  order  partial  correlation  coefficient  with  Y,  but  not  be  introduced 
into  the  equation  on  the  basis  of  its  intermediate  order  partial  correla¬ 
tion  coefficients.  Thus,  if  £  is  the  value  of  the  multiple  correlation 
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coefficient  after  the  introduction  into  the  regression  equation  of  q 
variables  from  p  independent  variables,  then  it  may  be  possible  to 
introduce  a  different  set  of  q  variables  into  the  equation  which  will 
yield  a  higher  value  of  R  .  (Yevjevich,  1972).  This  is  a  further 
reason  to  screen  the  independent  variables  for  intercorrelation,  before 
applying  the  step-wise  multiple  regression  procedure. 


3 . 6  Confidence  Intervals  On  Estimated  Values 

If  the  multiple  linear  regression  equation  is  used  to  predict 
a  single  value  (Y1)  of  the  dependent  variable  for  a  set  of  values  of  the 
p  independent  variables,  then  it  can  be  shown  (Brownless,  1960)  that 
the  variance  of  Y'  is  given  by:  - 


Var  (Y')  =  cr 


1  P  o  P"1 

—  +  X-,  (x.  -  x.)  c..  +  2  £ 
n  i=l  l  l  li  i=l 


(x.  -  x.)  c. . 
J  3  il 


P 

•  (x- 
j=i+l  i 


x±) 


Eq.  3.8 


where  c..  is  the  iith  element  in  the  variance-covariance  matrix 
1 3 

(i.e.,  the  inverse  matrix  of  the  normal  equation  coefficients), 
o' is  the  standard  deviation  of  the  population  random  components  and 
can  be  replaced  by  s,  the  unbiased  standard  error  of  the  estimate. 
This  leads  to  the  100(1  -c<)  confidence  intervals:  - 


YT  +  t1_(x/2  S  J  Var  (Yf)  Eq  .  3.9 

It  can  be  seen  from  equations  3.8  and  3.9  that  the  confidence 

limits  for  a  particular  estimated  value  of  Y’  are  not  unique.  A 
particular  value  of  Y’  can  result  from  an  infinite  number  of  different 
combinations  of  values  of  the  independent  variables  in  the  regression 
equation.  Each  combination,  while  yielding  the  same  mean  value  of  Y’ 
will  produce  different  values  for  the  confidence  limits. 


CHAPTER  4 


REGRESSION  ANALYSIS 

4 . 1  Watershed  Characteristics  included  in  Regression  Analysis 

The  drainage  basin  characteristics  included  in  the  regression 
analysis  follow  closely  those  proposed  by  Spence  (1971)  and  adopted  by 
Verschuren  (1973).  The  number  of  parameters  included  in  this  study 
was  reduced  and  the  method  of  determination  of  others  modified  as  a 
result  of  the  experience  gained  in  the  previous  works . 

The  parameters  included  met  the  following  requirements :  - 
i)  A  strong  possibility  that  each  parameter  had  a  direct 
physical  influence  on  the  rate  or  volume  of  runoff, 
and  ii)  Each  parameter  could  be  easily  and  accurately  determined 

from  1:250,000  topographic  maps.  (The  largest  scale 
maps  which  provide  a  complete  coverage  of  the  entire 
study  area  are  those  at  1:250,000). 

A  list  of  characteristics  and  the  abreviations  used  is  given 
in  Table  1. 

The  topographic  drainage  area  (TDA)  is  that  given  by  Water 
Survey  of  Canada  and  is  the  gross  drainage  area.  The  study  by  Spence 
indicated  that  the  refinement  of  subtracting  areas  which  are  internally 
draining,  to  yield  the  net  drainage  area,  is  unwarranted.  The  area  of 
the  watershed  has  a  strong  influence  on  the  volume  of  runoff  and  on  the 
peak  discharge. 

The  area  of  lakes  (LA)  was  calculated  as  a  percentage  of  the 
drainage  basin  area,  and  was  included  because  lakes  have  a  damping 
effect  on  a  flood  wave.  The  distribution  of  the  lakes  throughout  the 
basin  also  affects  damping,  but  a  single  factor  to  measure  the  distri- 
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bution  is  not  easy  to  determine,  and  not  included. 

Vegetative  cover  and  land  use  also  affect  the  volume  and  rate 
of  runoff  in  as  much  as  they  have  a  strong  influence  on  infiltration 
rates  and  evapotranspiration  losses.  However,  detailed  information  is 
not  readily  available,  so  the  only  vegetation  parameter  included  was  the 
area  occupied  by  forests  or  orchards  (FA)  expressed  as  a  percentage  of 
drainage  area.  This  information  was  measured  from  1:250,000  scale 
maps  and  is  consequently  imprecise. 

The  main  channel  length  (MCL)  was  defined  as  the  distance 
from  the  gauging  station  to  the  end  of  the  most  remote  tributary, 
measured  along  the  river  channel.  The  distance  was  determined  using  a  map 
measurer.  Meanders  too  small  to  be  measured  with  this  instrument 

were  ignored.  (The  limited  accuracy  with  which  the  MCL  could  be  measured 
did  not  warrant  the  introduction  of  the  closely  related  parameter  basin 
length  advocated  by  Spence) . 

The  basin  shape  factor  (BSF)  was  defined  as  the  square  of  the 
main  channel  length  divided  by  the  topographic  drainage  area.  It  provides 
an  indication  of  the  way  in  which  the  basin  area  is  distributed  in 
relation  to  the  main  channel^ which  in  turn  influences  the  rate  of  concen¬ 
tration  of  runoff  at  the  gauging  station  from  the  various  basin  parts. 
Using  the  above  definition,  two  basins  with  the  same  areal  shape,  slope 
and  main  valley  pattern,  will  have  different  values  of  basin  shape 
factor  according  to  the  degree  of  meandering  of  the  channel  within  the 
main  valley.  This  is  in  agreement  with  the  different  hydrograph  shapes 
expected  from  an  identical  input  over  two  such  basins. 

Basin  aspect  is  the  predominant  orientation  of  the  basin, 
measured  as  one  of  the  eight  points  of  the  compass.  It  is  introduced 
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as  two  parameters  to  take  account  of  its  cyclicity.  BA(X)  is  a  measure 
of  the  East/West  component  of  basin  aspect  and  BA(Y)  is  a  measure 
of  the  North/South  component.  Both  parameters  have  been  assigned 
arbitrary  ranges  of  0.0  to  +2.0.  Basin  aspect  has  an  effect  upon 
snowmelt  rates,  soil  moisture  conditions  and  vegetation  types,  which 
in  turn  influence  hydrograph  shape. 

Four  parameters  which  measure  watershed  elevation  are 
included  since  elevation  has  a  strong  influence  on  the  amount  and 
nature  of  precipitation  incident  on  an  area.  Gauge  elevation  (GE) 
and  maximum  elevation  (ME)  are  self-explanatory.  The  divide  elevation 
(DE)  is  the  elevation  of  the  point  where  the  main  channel  extended, 
intersects  the  watershed  boundary. 

The  elevation  of  two  points  were  determined;  10%  and  85% 
of  the  distance  along  the  main  channel,  measured  from  the  gauging 
station.  The  basin  elevation  (BE)  is  defined  as  the  mean  of  these 
two  elevations,  and  the  main  channel  slope  (MCS)  is  defined  as  the 
difference  between  these  two  elevations  divided  by  the  length  of 
channel  between  them. 

Basin  relief  (RE)  is  defined  as  the  difference  between  the 
gauge  elevation  and  the  maximum  elevation. 

The  foregoing  thirteen  parameters  were  initially  used  as 
independent  variables  in  the  regression  analysis.  As  the  analysis 
progressed,  it  became  apparent  that  possibly  some  improvement  could  be 
made  by  the  inclusion  of  further  watershed  characteristics.  The  follow¬ 
ing  two  parameters  were  therefore  determined:  - 

i)  the  latitude  of  the  centre  of  area  of  the  watershed 
expressed  in  degrees,  used  as  a  measure  of  the 
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geographic  location  of  the  basin,  (LAT) 
ii)  the  depth  of  rainfall  occurring  in  one  hour  which 
will  be  equalled  or  exceeded  on  average  once  in 
25  years  (RAIN) .  This  characteristic  is  not 
intended  to  represent  the  input  which  would  produce 
the  composite  hydrographs  described  in  Section  2, 
but  rather  to  group  and  separate  basins  into  diffe¬ 
rent  climatic  regions . 

The  fifteen  drainage  basin  characteristics  described  above 
were  determined  for  each  watershed  and  are  listed  in  Appendix  III. 

4 . 2  Regression  Analysis  and  Results 

An  attempt  was  made  to  determine  the  best  linear  relationship 
between  each  of  the  dependent  variables  (R(av) ,  R(max) ,  K(av) ,  K(max)  and 
BF)  and  the  selected  independent  variables,  using  a  stepwise  linear 
regression  program  (STPRG  step-wise  regression) .  Basically,  the 
program  operates  as  follows :  - 

In  the  first  step,  the  regression  equation  is  found  between 
the  dependent  and  one  independent  variable;  in  each  subsequent  step  one 
variable  is  added  and  the  new  regression  equation  determined.  (For 
fuller  explanation,  see  Section  3.5).  The  multiple  correlation  coeffi¬ 
cient  is  determined  at  each  step  together  with  its  unbiased  standard 
error  of  estimate,  the  regression  equation  coefficients  and  their 
corresponding  standard  errors  of  estimate.  In  general,  independent 
variables  are  introduced  into  the  analysis  until  the  unbiased  standard 
error  of  the  estimate  reaches  a  minimum  value.  However,  if  the  correla¬ 
tion  between  the  dependent  variable  and  the  last  independent  variable 
is  weak,  then  the  associated  regression  coefficient  may  be  statistically 
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insignificant,  in  which  case,  the  last  step  should  be  deleted  and  the 
analysis  terminated.  This  will  result  in  a  legitimate  value  of  the 
multiple  correlation  coefficient  somewhat  lower  than  the  value  obtained 
in  the  final  step. 

The  technique  is  purely  an  arithmetic  one  and  does  not  assess 
whether  the  relationship  given  by  the  regression  equation  is  physically 
realistic. 

In  order  to  obtain  the  highest  significant  value  of  the 
multiple  correlation  coefficient  (i.e.,  the  best  linear  relationship), 
several  separate  runs  using  program  SSPLIB  were  carried  out  for  each 
of  the  five  dependent  variables  (R(av) ,  R(max) ,  K(av) ,  K(max)  and  BF) . 

To  illustrate  this  sequence,  the  following  is  a  brief  description  of 
each  run  using  the  average  rate  of  rise  R(av)  as  the  dependent  variable. 
A  complete  list  of  regression  equations  for  R(av)  is  given  in  Table  2. 

The  first  run  included  data  from  all  62  stations  which  met 
the  requirements  of  Section  1.4.  The  first  13  variables  of  Appendix 
III,  i.e.,  those  discussed  in  Section  4.1,  were  made  available  for 
inclusion  into  the  regression  equation,  step  by  step. 

Significance  tests  indicated  that  the  first  four  regression 
coefficients  only  were  significant  at  the  5%  level.  The  resulting 
equation  is  Equation  4.1  and  the  corresponding  value  of  the  multiple 
correlation  coefficient  (R)  is  0.65. 

For  the  second  run,  common  logarithms  were  taken  of  all  the 
variables,  and  again,  all  the  independent  variables  were  available  for 
selection.  It  can  be  seen  from  Table  2  that  the  multiple  correlation 
coefficient  increased  to  0.71.  To  enable  logarithms  to  be  taken,  all 
the  independent  variables  were  screened  and  any  zero  values  were 
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arbitrarily  set  to  0.1.  This  was  well  within  the  accuracy  with  which 
the  variables  had  been  determined. 

Prior  to  making  the  third  run,  thirteen  separate  runs  were 
made  using  logarithms  of  variables,  each  run  having  one  independent 
variable  available  for  selection.  This  resulted  in  thirteen  separate 
equations  of  the  form:  - 

B 

R(av)  =  A  x  X  n 
n  n 

with  n  =  1,2 . . .13, 


X^  representing  the  independent  variables, 
and  A  and  B  -  constants. 

For  the  third  run,  the  independent  variables  were  raised  to 
the  respective  B^  powers.  The  result  is  Eq.  4.3  and  the  value  of  the 
multiple  correlation  coefficient  is  0.74. 

In  an  attempt  to  obtain  a  regression  equation  which  was 
physically  realistic,  the  number  of  independent  variables  available  for 
selection  in  the  fourth  run  was  restricted  to  eight.  They  were  the 
ones  thought  most  likely  to  affect  hydrograph  shape  and  which  had,  in 
addition,  relatively  low  zero-order  correlation  coefficients  between 


themselves  ( 

r.  .1  <  0.44)  .  The  variables 
iJ1 

available 

1. 

Topographic  Drainage  Area 

(TDA) 

2. 

Lake  Area 

(LA) 

3. 

Forest  Area 

(PA) 

4. 

Basin  Shape  Factor 

(BSF) 

5. 

Basin  Aspect  (X  direction) 

(BA(X) ) 

6. 

Basin  Aspect  (Y  direction) 

(BA(Y)) 

7. 

Basin  Relief 

(RE) 

8. 

Main  Channel  Slope 

(MCS) 

■ 


39. 


the  regression  equation  resulting  from  this  run  is: 

R(av)  =  0.18  +.000009  (TDA)  +.00004(RE)  -.013(BSF)  -.0067(MCS) 

+.0016 (FA)  Eq.  4.4 

with  multiple  correlation  coefficient  of  0.69.  Noting  that  a  large 
value  of  rate  or  rise  R(av)  indicates  a  flat  slope  for  the  rising  limb 
of  the  hydrograph  and  that  basin  shape  factor  is  defined  as  main  channel 
length  divided  by  topographic  drainage  area,  it  can  be  seen  that  the 
sign  of  the  coefficients  of  TDA,  MCS  and  FA  appear  physically  realistic, 
but  those  of  RE  and  BSF  do  not. 

The  fifth  run  was  similar  to  that  just  described  except  that 
logarithms  of  variables  were  used.  The  multiple  correlation  coefficient 
increased  slightly  to  0.71,  MCS  was  not  introduced  into  the  equation,  and 
the  sign  of  the  coefficients  of  BSF  and  RE  remained  "wrong".  Thereafter, 
attempts  to  keep  the  equation  physically  realistic  were  abandoned. 

It  appeared  from  the  results  of  the  first  five  runs  that  at 
most,  six  independent  variables  contributed  significantly  to  the  reduction 
in  variation  of  the  dependent  variable,  so  for  the  next  run,  the  seven 
independent  variables  which  had  the  lowest  zero-order  correlation  coeffi¬ 
cients  between  themselves  were  included.  A  low  value  of  the  multiple 
correlation  coefficient  of  0.64  resulted. 

It  was  apparent  that  zero-order  correlation  coefficients  did 
not  give  a  true  indication  of  the  correlation  between  variables.  There¬ 
fore,  a  small  program  was  written  to  calculate  the  matrix  of  twelfth 

order  partial  correlation  coefficients  (see  Section  3.1)  which  gave  a 
measure  of  the  true  correlation  between  two  variables  with  the  effects 

of  all  other  variables  removed. 

The  seventh  run  was  made  with  eight  independent  variables 
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available  for  selection.  The  selected  variables  had  high  twelfth 
order  partial  correlation  coefficients  with  the  dependent  variable 
and  were  not  highly  correlated  to  each  other.  Five  variables  which 
did  not  meet  these  criteria  were  eliminated.  The  run  was  made  with 
arithmetic  variables,  and  a  multiple  correlation  coefficient  value 
of  0.72  resulted. 

For  the  next  run,  common  logarithms  were  taken  of  the  variables 
and  the  twelfth  order  partial  correlation  matrix  calculated.  From  this, 
the  eight  variables  were  selected  for  inclusion  using  the  same  criteria 
as  for  the  previous  run.  The  multiple  correlation  coefficient  in¬ 
creased  to  0.74.  Three  variables  from  the  previous  equation  also 
appeared  in  this  result. 

The  ninth  run  was  made  using  seven  independent  variables  raised 
to  powers  (as  discussed  for  run  No.  3) .  The  seven  variables  included 
were  selected  from  a  new  matrix  of  twelfth  order  partial  correlation 
coefficients  using  the  same  criteria  as  before.  The  multiple  correlation 
coefficient  reached  a  high  of  0.75. 

It  can  be  seen  from  Table  2  that  TDA  was  the  first  variable  to 
be  included  in  the  regression  equation  for  each  of  the  nine  runs.  This 
is  so  because,  of  the  particular  set  of  independent  variables  included 
in  each  run,  the  relevant  function  of  TDA  had  the  highest  zero  order 
correlation  coefficient  with  the  dependent  variable  (see  Section  3.5) . 

The  regression  equation  resulting  from  the  ninth  run  was  the 
best  that  had  so  far  occurred,  in  as  much  as  the  multiple  correlation 
coefficient  was  the  highest,  and  it  appeared  that  little  improvement 
could  be  made  using  the  independent  variables  already  available.  Con¬ 
sequently,  a  plot  was  made  of  observed  (i.e.,  obtained  from  the  composite 
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hydrographs)  values  of  R(av)  versus  values  calculated  using  Equation  4.9. 
There  was  considerable  scatter  about  a  line  bisecting  the  axes  and  passing 
through  the  origin.  In  order  to  see  if  any  pattern  existed,  the  62 
stations  were  divided  into  five  groups  on  the  basis  of  their  geographic 
location,  vis:  foothills,  prairies,  aspen  parkland,  boreal  forest, 

Yukon  and  N.W.T.  A  symbol  was  selected  to  represent  each  group  and  the 
symbol  appropriate  to  each  station  location  was  marked  adjacent  to 
the  points  on  the  plot.  It  was  apparent  that  the  values  calculated 
using  the  regression  equation  were  too  high  for  those  stations  in  the 
prairies,  aspen  parkland  and  foothills,  and  to  a  lesser  extent,  too 
low  for  those  stations  in  the  boreal  forest.  The  northern  stations 
were  widely  scattered. 

In  order  to  improve  the  regression  equation,  a  parameter  or 
parameters,  were  required  that  would  connect  the  stations  of  the  prairies, 
aspen  parkland  and  foothills,  and  distinguish  them  from  the  stations  of 
the  boreal  forest.  Station  latitude  was  chosen  for  this  purpose,  although 
it  was  appreciated  that  the  foothills  spanned  a  considerable  latitude 
range,  and  that  the  scatter  of  the  Yukon  and  N.W.T.  stations  on  the  plot 
would  not  be  reduced  since  the  station  latitude  values  would  all  fall 
within  a  narrow  band.  In  addition,  it  was  realized  that  latitude 
probably  had  no  direct  effect  on  hydrograph  shape,  but  that  some  topo¬ 
graphical  or  climatological  condition  may  be  associated  with  latitude, 
which  would  have  an  influence. 

A  parameter  which  more  clearly  divided  the  stations  in  the 
desired  manner  was  found  to  be  the  60  minute  rainfall  depth  having  a 
return  period  of  25  years.  The  inclusion  of  this  parameter  is  justified 
because  the  hydrograph  parameters  being  considered  are  for  runoff  events 
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of  low  probability  and  therefore  25  years  return  period  rainfall  could 
be  relevant.  Also,  many  of  the  isohyetal  maps  for  different  durations 
and  return  periods  had  similar  general  shapes  and  therefore  the  choice 
was  not  critical. 

The  values  of  the  two  additional  independent  variables  (LAT  and 
RAIN)  are  given  in  Appendix  III  for  each  station. 

The  tenth  run  was  made  using  variables  raised  to  powers. 

Using  all  fifteen  independent  variables,  the  matrix  of  fourteenth  order 
partial  correlation  coefficients  was  determined.  Station  latitude  was 
not  well  correlated  with  the  dependent  variable  and  was  not  included 
amongst  the  seven  independent  variables  available  in  the  regression 
analysis.  The  resulting  equation  (Eq.  4.10  Table  2)  has  the  main 
channel  slope  as  the  first  variable,  the  new  rainfall  parameter  (RAIN) 
next  and  the  drainage  area  last.  The  multiple  correlation  coefficient 
of  0.76  is  a  slight  improvement  over  the  previous  highest. 

Another  plot  was  made  of  "observed"  versus  calculated  values 
of  R(av)  using  the  new  equation.  A  slight  readjustment  of  the  grouping 
of  the  stations  occurred  but  there  was  still  considerable  scatter. 

It  was  noticed  that  the  two  stations  which  had  the  greatest 
difference  between  "observed"  and  calculated  values  (Stations  32  and 
42)  were  at  or  near  the  outlet  of  large  lakes.  Three  other  stations 
(44,  45  and  46)  were  found  to  be  similarly  situated,  and  in  each  case 
the  "observed"  value  of  R(av)  was  considerably  greater  than  the  calculated 
value.  This  is  due  to  the  reduction  of  the  flood  peak  as  it  is  routed 
through  the  lake,  which  results  in  a  flatter  hydrograph  and  consequently 
higher  values  of  R(av) .  It  was  considered  that  the  stations  close  to  lake 
outlets  belonged  to  a  different  population,  and  therefore  it  was  legiti- 
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mate  to  exclude  these  stations  from  the  analysis. 

The  final  run  in  this  series,  run  eleven  was  made  using  only 

the  data  from  the  57  non-lake-outlet  stations.  The  independent  variables 

were  raised  to  powers  as  in  the  previous  run.  The  resulting  equation  is:- 

R(av)  =  -.48+.026(RAIN)-1'46+.001(TDA) ‘6°+.025(FA) ,35-.0081(MCL) *64 
50  -  7? 

+.0036(RE)  +.063(MCS)  Eq.  4.11 

and  the  multiple  correlation  coefficient  is  0.80;  the  best  value 

/ 

obtained. 

Initially,  runs  were  made  with  R(av) ,  R(max) ,  K(av)  and  K(max) 
as  dependent  variables.  However,  it  became  apparent  that  the  values  of 
the  multiple  correlation  coefficients,  associated  with  R(av)  and  R(max) 
were  very  similar  for  each  run.  The  same  applied  to  K(av)  and  K(max) . 

In  order  to  avoid  duplication,  R(max)  and  K(max)  were  eliminated  from 
the  analysis. 

An  analysis  similar  to  that  described  above  was  carried  out  for 
the  two  independent  variables;  average  rate  of  fall  K(av)  and  baseflow 
BF.  In  general,  at  each  stage  of  the  analysis,  the  multiple  correlation 
coefficient  associated  with  the  K(av)  equation  was  somewhat  less  than 
that  corresponding  to  R(av) ,  while  that  associated  with  BF  was  somewhat 
greater.  The  final  equations  are  given  in  Table  3. 

It  is  noteworthy  that  the  rainfall  parameter  (RAIN)  and 
drainage  area  (TDA)  occur  in  the  final  equations  for  R(av) ,  K(av)  and 
BF,  and  that  basin  relief  (RE)  occurs  in  BF  and  R(av) . 

4 . 3  Map  Coefficients 

In  order  to  make  equations  4.11,  4.12  and  4.13  useful  for 
predictive  purposes,  maps  were  prepared  showing  lines  of  equal  values 
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of  the  regression  coefficient  b  ,  or  intercept.  The  method  is  similar 
for  each  parameter  so  the  case  of  R(av)  only  will  be  discussed. 

Consider  a  location  on  an  ungauged  river  within  the  geographic 
region  of  the  sample  stations,  and  with  a  drainage  area  between  500  and 
15,000  square  miles.  Given  Equation  4.11  and  the  values  of  the  five 
parameters  contained  in  it  (determined  from  meteorological  and  topo¬ 
graphical  maps) ,  the  value  of  R(av)  can  be  determined  for  the  ungauged 
station.  However,  the  result  is  not  very  reliable  due  to  the  only 
moderate  value  of  the  multiple  correlation  coefficient  (0.80).  The 
value  of  R(av)  can  be  improved  by  using  a  map  coefficient  (c)  in  place 

of  the  intercept  (b  ) .  The  map  coefficient  is  an  amended  value  of  b  . 

o  o 

The  amendment  is  made  to  account  for  the  effect  of  parameters  not 
included  in  the  analysis,  e.g.,  soil  properties,  ratio  of  rainfall  to 
snowfall,  etc.  A  map  showing  lines  of  equal  values  of  c  is  prepared 
as  follows:  - 

On  an  outline  map  of  the  geographic  region  of  the  sample 

stations,  the  station  number  is  marked  at  the  centre  of  its  drainage 

basin.  (The  hydrograph  parameters  reflect  the  characteristics  of  the 

entire  basin  and  not  just  those  of  the  gauging  site)  .  Using  Equation 

4.11,  the  values  of  R(av)  are  calculated  for  each  station,  using  the 

previously  determined  parameters,  and  subtracted  from  the  "observed" 

R(av)  values  of  Appendix  II.  This  difference  is  then  added  to  the 

respective  b  value  to  produce  the  map  coefficient  c,  which  is  then 

marked  on  the  map  adjacent  to  the  station  number.  Lines  of  equal  c 

values  are  then  drawn,  based  on  the  values  at  each  station.  A  c  value 

for  any  ungauged  station  can  then  be  read  from  the  map  and  substituted 

for  b  in  the  regression  equation.  Maps  of  Alberta  and  Yukon  Territory 
o 
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showing  map  coefficients  for  R(av) ,  K(av)  and  BF  are  given  in  Appendix 
IV.  If  the  contours  could  be  drawn  exactly  as  dictated  by  the  spot  • 
values,  then  the  calculated  values  of  R(av)  for  the  sample  stations  would 
equal  the  "observed"  values.  However,  in  order  to  obtain  reasonably 
smooth  contours  some  deviation  from  the  calculated  c  values  is  necessary 
which  accounts  for  the  scatter  on  the  plots  of  observed  versus  calculated 
values  of  R(av) ,  K(av)  and  BF  (Figs.  4,  5  and  6). 

4 . 4  Regional  Analysis 

For  this  part  of  the  study,  the  42  Alberta  sample  stations 
were  sub-divided  on  a  geographical  basis  into  4  groups  -  prairies, 
foothills,  aspen  parkland  and  boreal  forest  (See  Appendix  III).  This 
resulted  in  groups  of  8,  12,  7  and  15  stations  respectively.  Stepwise- 
multiple  linear-regression  analyses  were  carried  out  on  certain  groups 
and  combinations  of  groups  in  turn,  with  R(av)  as  the  dependent  variable 
and  the  15  independent  variables  shown  in  Table  1  available  for 
selection  into  the  regression  equation.  The  results  are  given  in 
Table  4.  While  the  multiple  correlation  coefficient  of  0.98  obtained 
for  the  prairie  stations  appears  very  high,  it  must  be  remembered  that 
the  sample  consisted  of  only  eight  data  points.  Similarly,  the  foothills 
samples  contained  only  12  stations. 

It  would  be  unwise  to  use  the  regression  equations  for 
predictive  purposes  because  the  confidence  bands  are  very  broad  due  to 
the  small  sample  size.  However,  the  indications  are  that,  if  a  large 
sample  from  one  small  geographic  region  could  be  obtained,  then  good 
correlations  might  be  established. 
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OBSERVED  VALUES 


FIG.  4 


R  (AV):  OBSERVED  VS.  ESTIMATED  VALUES  BY  REGRESSION 
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OBSERVED  VALUES 


FIG.  5 


K  (AV):  OBSERVED  VS.  ESTIMATED  VALUES  BY  REGRESSION 


ESTIMATED  VALUES 
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OBSERVED  VALUES 


FIG.  6 


B/FLW 


OBSERVED  VS.  ESTIMATED  VALUES  BY  REGRESSION 
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CHAPTER  5 


PEAK  DISCHARGE 


5.1  Theory 

It  was  demonstrated  in  Section  2.4  that  the  composite 
hydrographs,  produced  by  combining  the  steepest  portions  of  recorded 
floods  hydrographs,  when  ploted  as  log  Q  versus  time,  resulted  in 
straight  lines  over  considerable  discharge  ranges.  In  addition,  there 
often  appeared  a  quite  distinct,  almost  horizontal  portion,  in  the 
lowest  discharge  range  which  has  been  called  baseflow. 

Consider  the  idealized  hydrograph  of  Fig.  7.  Fig.  7(a) 
has  an  arithmetic  time  scale  and  logarithmic  discharge  scale.  Let  the 
rising  and  falling  limbs  be  straight  lines  and  terminate  at  a  constant 
discharge  value  Q^.  Let  the  peak  discharge  occur  at  t=0.  Fig.  7(b) 
shows  the  same  hydrograph  drawn  to  an  arithmetic  scale  on  both  axes. 

Consider  the  general  equation  of  a  straight  line  on  semi- 
logarithmic  axes. 


log  Q  =  m.t  +  c  Eq.  5.1 

where  Q  and  t  are  the  variables 
c  is  the  intercept 
and  m  is  the  gradient 
Referring  to  Fig.  7(a), 
c  =  log  Qq 

-(log  Qn  -  log  Qn  +  p 

and  m  =  - 

t  , ,  -  t 
n+1  n 

Substituting  for  m  in  Eq .  5.1  gives 
log  Q  =  t (log  Qn+1  -  log  Qn)  +  log  Qq 
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(a) 


(b) 

FIG.  7  IDEALISED  HYDROGRAPH  SHAPE 
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or 


Q  =  Q, 


<n+l 


n 


Eq.  5.2 


For  the  falling  limb  of  a  hydrograph,  i.e.,  negative  gradient, 

Qn+i 

define  — —  =  K  =  rate  of  drop  of  discharge  =  -  ^77- 
n 


then  Q  =  Q  (K) 
o 


Eq.  5.3 


For  the  rising  limb  of  a  hydrograph,  i.e.,  positive  gradient, 


define 


R  -  rate  of  rise  of  discharge  = 


1_ 

R' 


then  Q  =  Qq  (R)  1  Eq.  5.4 

As  defined  above,  K  is  the  ratio  of  the  discharge  on  a  par¬ 
ticular  day  to  the  discharge  on  the  previous  day.  R  is  the  ratio  of 
discharge  on  a  particular  day  to  the  discharge  on  the  following  day. 

The  range  of  possible  values  of  R  and  K  is  0.0  to  1.0. 

The  area  under  a  hydrograph  in  a  particular  time  interval, 
represents  the  volume  of  runoff  in  that  time.  The  volume  of  direct 
runoff  is  obtained  by  subtracting  the  baseflow  volume  from  the  total 
volume.  Referring  to  Fig.  7  (b) ,  the  area  A’B'F'  represents  the  volume 
of  direct  runoff  in  the  time  interval  t ^  -  t^. 

In  order  tb  obtain  an  expression  for  the  volume  of  direct 
runoff,  consider  first  the  volume  of  direct  runoff  during  the  rising 
period  of  the  hydrograph. 

The  volume  of  direct  runoff  in  the  time  interval  t  =  t^  to 
t  =  0  is  equal  to  the  area  under  the  rising  limb  of  the  hydrograph 
and  bounded  by  the  line  Q  =  . 
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The  rising  limb  equation  is  Q  =  Qq.(R) 


-t 


i .  e . ,  In  Q  =  In  Q  -  t .  In  R 

o 


i.e. ,  t  = 


In  Q  In  Q 
o 


In  R  In  R 
The  volume  of  direct  runoff  is  given  by 


V 


(tro) 


1 1  .dQ 


In  Q 
_ jc 

In  R 


dQ  + 


In  Q 
In  R 


dQ 


“  Q .In  Q 

o 

+ 

Q .  In  Q  -  Q 

In  R 

n 

In  R 

which  leads  to 
1 


V 


(t][-0)  In  R 


Qx  (In  Qq  -  In  Q±)  +  Qx  -  Qq 


Similarly,  the  volume  of  direct  runoff  during  the  recession 
period,  i.e.,  in  the  time  interval  t  =  0  to  t  =  t^  is  given  by 


^Qx  (In  Qq  -  In  Q ±)  +  -  Qq 


t2)  In  K 

Thus,  the  total  volume  of  direct  runoff 


v  v(t1-o)  +  v(o-t2) 


i.e.,  V  = 


Q,  (In  Q  -  In  Q-)  +  Q_  -  Q 
1  o  1  1  o 


1  +  1 


Eq.  5.5 


In  K  '  In  R_ 

The  direct  runoff  volume  V  can  be  calculated  for  a  watershed 
in  terms  of  an  average  depth  of  exess  precipitation  per  unit  area  (d) 
and  the  drainage  basin  area  (TDA) 


Eq.  5.6 
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i.e.,  V  =  d.(TDA) 

Combining  Equations  5.5  and  5.6,  and  rearranging 


Q  .In  Q  -  Q  =  d.(TDA) .(In  R  x  In  K)  +  Q  .(In  Q  -  1)  Eq .  5.7 
°  °  (In  R  +  In  K)  1  1 


The  above  is  a  non-linear  equation  expressing  the  peak  dis¬ 
charge  Qq  in  terms  of  the  average  depth  of  excess  precipitation,  basin 
area,  the  rates  of  rise  and  fall  of  the  hydrograph  and  the  value  of 
baseflow  discharge  Q^. 


5 . 2  Regression  Analysis  with  Peak  Flow  as  the  Dependent  Variable 

For  each  of  the  57  non-lake-outlet  study  watersheds,  a  value 
of  peak  flow  was  calculated  using  Eq .  5.7,  for  an  average  depth  of 
excess  precipitation  from  the  basin  of  1  ins.  (Qq  (1")) .  It  should 
be  noted  that  the  volume  of  runoff,  i.e.,  d  x  TDA  must  be  expressed  in 
units  of  cf  s .  days,  i.e.,  V  £cfs.  days^J  =  26.8  x  d  x  TDA  where  d  Tins 
and  TDA  [sq.  mis 

Qq  (1m)  calculated  using  Eq .  5.7  is  referred  to  hereafter  as 
the  'observed'  value  of  peak  flow  for  1  ins.  of  excess  precipitation. 
It  should  not  be  confused  with  the  actual  recorded  peak  flow  of  the 
composite  hydrograph.  It  is  a  calculated  discharge  value  based  on 
the  observed  steepest  rising  and  falling  limbs  of  the  composite  hydro¬ 
graph  . 

Qq  (1")  was  used  as  the  dependent  variable  in  a  stepwise 
regression  analysis.,  with  the  basin  characteristics  of  Appendix  III  as 
independent  variables . 

The  equation  which  gave  the  highest  value  of  the  multiple 
correlation  coefficient  (R=0.88)  is:  - 

Qo(lM)  =  44.0  x  (TDA)0,85  x  BA(Y)°*2  x  (RAIN)0*63  Eq .  5.8 
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The  author  analyzed  hydrographs  based  on  three  years  of 
record,  from  8  small  (less  than  165  sq .  mis.)  watersheds  in  the  vicinity 
of  Watson  Lake,  Y.T.  (Verschuren  and  Bristol  1974).  The  hydrographs 
were  based  on  instantaneous  discharges  and  were  analyzed  manually  to 
determine  the  maximum  rates  of  rise  and  fall.  No  account  was  taken 
of  base  flow  in  the  study,  but  it  is  known  that  winter  flows  are 
extremely  small  and  the  effect  of  the  omission  is  probably  negligible. 

Peak  discharge  resulting  from  1"  of  excess  precipitation 
was  calculated  for  the  8  small  basins  and  plotted  against  basin  area, 
on  logarithmic  graph  paper.  The  equation  of  the  resulting  straight  line 
envelope  was  found  to  be: 

Qo(l”)  =  48(TDA)*72  Eq.  5.9 

The  similarity  between  this  equation  and  Equation  5.8,  derived  for 
much  larger  basins,  is  noteworthy. 

To  enlarge  the  range  and  size  of  the  sample,  the  drainage  basin 
characteristics  of  these  8  small  basins  were  calculated  (Station  Nos. 

63  -  70  of  Appendix  III)  and  another  regression  analysis  was  carried  out 
on  the  sample  of  65  stations  (lake  outlet  stations  once  again  being 
omitted) .  The  results  of  the  analyses  for  various  depths  of  excess 
precipitation  are  given  in  Table  5,  together  with  values  of  the  multi¬ 
ple  correlation  coefficients. 

Table  6  gives  the  regression  equations  obtained  at  the  first 
step  of  the  stepwise  analysis.  It  will  be  noted  that  topographic 
drainage  area  is  the  parameter  included  and  that  for  depths  of  excess 
precipitation  of  one  inch  or  greater,  it  occurs  raised  to  the  power  of 
0.88.  In  addition  ,  it  should  be  noted  that  the  multiple  correlation 
coefficients  are  at  most  0.03  less  than  the  values  obtained  in  Table  5. 
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TABLE  5  REGRESSION  EQUATIONS  OF  PEAK  DISCHARGE  FOR  VARIOUS  DEPTHS  OF  EXCESS 

PRECIPITATION 
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Fig.  8  shows  plots  of  the  intercept  of  the  regression  equation 
versus  depth  of  excess  precipitation  for  the  equations  of  Tables  5  and  6. 

It  can  be  seen  that  the  intercept  increases  virtually  linearly  with 
excess  precipitation. 

Fig.  9  shows  a  plot  of  'observed'  Q  (1")  values  (i.e.,  cal¬ 
culated  using  Eq.  5.7)  versus  values  calculated  using  the  regression 
equation  given  in  Table  5  (Eq.  5.12).  Because  logarithms  of  variables 
were  used  in  the  regression  analysis,  there  is  considerable  scatter  of 
the  points  in  spite  of  the  high  multiple  correlation  coefficient  of 
0.97.  In  order  to  make  the  regression  equation  more  useful  for  peak 
flow  prediction,  map  coefficients  have  been  determined  (Appendix  IV) . 

In  order  to  utilize  one  set  of  map  coefficients  for  the  six  regression 
equations  derived  for  the  six  depths  of  excess  precipitation  shown  in 
Table  5,  the  'contour'  values  for  Qq  shown  in  Appendix  IV  are  factors  by 
which  the  regression  equations  are  to  be  multiplied  (rather  than  intercept 
values  as  in  the  case  of  R(av) ,  K(av)  and  BF  map  coefficients) .  A  plot 
of  'observed'  versus  calculated  values  of  Q  (1")  adjusted  using  the  map 
coefficients  is  given  in  Fig.  10. 

5 . 3  Peak  Discharge  vs  Topographic  Drainage  Area 

To  assess  the  value  of  Equation  5.18  in  Table  6,  logarithms  of 
peak  flow  due  to  1  ins.  of  excess  precipitation  (Q  (1"))  was  plotted 
against  logarithms  of  drainage  basin  area  (TDA)  (Fig.  11).  In  spite  of 
the  high  multiple  correlation  coefficient  (0.94),  the  predicted  value  of 
peak  flow  varies  by  as  much  as  a  factor  of  10  for  a  given  drainage  basin 
area.  The  values  of  the  random  variable  plotted  in  Fig.  11  were 
calculated  using  Eq .  5.7,  and  therefore,  for  a  given  drainage  area,  effects 
on  peak  flow  due  to  differences  in  hydrograph  shape  did  not  manifest 
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FIG.  8 


REGRESSION  INTERCEPT  VS.  DEPTH  OF  EXCESS  PRECIPITATION 
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OBSERVED  VALUES 


FIG.  9 


PEAK  FLOW  FROM  I  INS  EXCESS  PPT.:  OBSERVED  VS.  ESTIMATED 
VALUES  BY  REGRESSION 
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OBSERVED  VALUES 


FIG.  10 


PEAK  FLOW  FROM  I  INS  EXCESS  PPT.:  OBSERVED  VS.  ESTIMATED 
VALUES  BY  REGRESSION,  ADJUSTED  USING  MAP  COEFFICIENTS 
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FIG.  II 


ESTIMATED  PEAK  FLOW  FROM  I  INS.  EXCESS  PPT.  VS.  BASIN  AREA 
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themselves.  The  spread  in  values,  for  a  given  basin  area  and  depth 
of  excess  precipitation,  must  be  attributable  to  the  rates  of  rise  and 
drop  of  discharge  and  to  baseflow  (R,K,  and  BF) . 

To  investigate  the  effect  of  different  R  and  K  values, 

Eq.  5.12  was  used  to  calculate  Q  (1")  values  for  each  of  the  study  water¬ 
sheds,  using  the  known  basin  characteristics.  The  Q  (1")  values  so 

o 

calculated  are  called  the  estimated  values.  The  ratio  was  found  between 
these  estimated  Q^Cl")  values,  and  the  Q  (1")  values  determined  using 
Eq.  5.7  (the  ’observed'  values).  The  stations  for  which  the  ratio 
was  greater  than  1.2  or  less  than  0.8  were  noted.  Next,  a  table  was 
prepared  (Table  7) ,  showing  the  bi-variate  distribution  of  R  and  K  for 
each  of  the  sample  stations.  Stations  for  which  the  ratio  of  estimated 
to  observed  Q  (1”)  values  was  less  than  0.8  were  underlined  and  those  for 
which  the  ratio  was  greater  than  1.2  were  double  underlined.  Table  7 
shows  a  strong  trend;  that  for  those  stations  having  low  values  of  R  and 
K,  i.e. ,  steep  rising  and  falling  hydrograph  limbs,  Eq.  5.12  underestimates 
the  value  of  Q  .  Similarly,  Eq.  5.12  overestimates  the  value  of  Qq  for 
those  stations  with  flat  hydrographs . 

On  the  basis  of  Table  7,  the  sample  stations  were  separated 
into  3  groups;  group  1  consisting  of  23  stations  with  low  values  of  both 
R  and  K,  group  3  consisting  of  12  stations  with  high  values  of  R  and  K, 
and  group  2  containing  the  remaining  30  stations.  Regression  analyses 
were  then  carried  out  on  each  group  using  basin  area  as  the  independent 
variable.  The  results  are  given  in  Table  8  and  plotted  in  Fig.  12.  The 
three  multiple  correlation  coefficients,  each  in  excess  of  0.97,  indicate 
the  high  degree  of  correlation  between  peak  flow  and  drainage  basin  area. 
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TABLE  7  BI-VARIATE  DISTRIBUTION  OF  R  AND  K  FOR  65  SAMPLE 

GAUGING  STATIONS 
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TABLE  8  REGRESSION  EQUATIONS  FOR  PEAK  FLOW  VS.  TOPOGRAPHIC  DRAINAGE  AREA,  FOR  RANGES 

A  P-  n  /  \  1  ..  A  /  \ 
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FIG.  12 


PEAK  FLOW  FROM  I  INS  OF  EXCESS  PPT.  VS.  TOPOGRAPHIC 
DRAINAGE  AREA,  FOR  RANGES  OF  R  AND  K 
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The  equations  of  Table  8  can  be  extended  and  simplified  to 

Q  =  C-  x  d  x  A  Eq .  5.25 

o  1 

where  d  =  average  depth  of  excess  precipitation  (ins.) 

C1  =30.5  for  basins  with  steep  hydrographs 
=  14.3  for  basins  with  mixed  hydrographs 
=  5.5  for  basins  with  flat  hydrographs 

5 . 4  Comparison  with  Existing  Peak  Flow  Formulae 

Many  well  known  peak  flow  formulae  are  of  the  form 

Q  =  mAn  Eq .  5.26 

where  Q  =  the  peak  discharge 
A  =  the  basin  area 
M,  n  =  are  coefficients 

In  order  to  compare  the  form  of  Eq .  5.25  with  that  of  Eq . 

5.26,  it  was  necessary  to  eliminate  the  parameter  d  (the  maximum  depth  of 
excess  precipitation  occurring  over  the  drainage  basin  area) .  Meteorolo¬ 
gical  records  for  Alberta  were  searched  to  find  the  most  severe  Maximum 
Depth  Area  Curve  for  the  region.  The  greatest  rainfall  depth  resulted 
from  a  storm  centred  50  miles  southeast  of  Banff  from  June  19  -  29,  1969. 
The  depth  of  precipitation  for  the  total  storm  duration  is  plotted  on 
logarithmic  axes  in  Fig.  13.  It  can  be  seen  that  the  curve  can  be 
approximated,  for  basin  areas  less  than  10,000  sq .  miles,  by  the  equation 

P  =  14.5  a"0,08  Eq.  5.27 

For  basin  areas  between  10,000  sq .  miles  and  100,000  sq .  miles 

the  depth  of  precipitation  is  given  by  the  equation 

-  41 

P  =  320  A 


Eq.  5.28 
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The  depth  of  excess  precipitation  (d)  equals  the  depth  of 
precipitation  (P)  multiplied  by  a  runoff  coefficient  C.  C  may  be  assumed 
approximately  equal  to  0.5 

Then  d  =  0.5  x  P 

Substituting  for  P  in  Eq .  5.27,  for  basins  less  than  10,000  sq .  miles, 

d  =  7.25A-0'08  Eq.  5.29 

and  substituting  for  P  in  Eq .  5.28,  for  basins  greater  than  10,000  sq. 
miles 

d  =  160A-0*41  Eq.  5.30 


Substituting  for  d  in  Eq .  5.23,  for  basins  less  than  10,000 

sq .  miles,  whose  hydrographs  have  intermediate  slopes 

rv  T  /  o  -7  A  *97 

Q  =  14.3  x  7.25A  x  A 

o 


i.e. ,  Q  -  100A 

o 


.9 


Eq.  5.31 


Similarly,  for  basins  greater  than  10,000  sq.  miles 

-  41  97 

Q  =  14.3  x  160A  x  A 
o 


i.e.,  Q  —  2 ,300  A 
o 


.56 


Eq.  5.32 


CHAPTER  6 


DISCUSSION  &  CONCLUSIONS 

6 . 1  Discussion 

The  result  of  this  study  is  a  procedure  which  enables  the 
determination  of  the  main  features  of  the  hydrographs  of  rare  flood 
events,  for  ungauged  watersheds  from  1  to  15,000  sq.  miles  in  area, 
located  in  Alberta  and  southern  Yukon  Territory. 

An  important  feature  of  this  study  is  that  hydrometric 
records  of  the  sample  watersheds  were  analyzed  without  the  need  for 
information  about  the  input  during  the  period  of  record.  Such  informa¬ 
tion  is  not  usually  readily  available.  The  equations  derived  allow  the 
determination  of  the  approximate  hydrograph  shape  for  different  volumes 
of  direct  runoff.  In  order  to  apply  the  method  to  ungauged  watersheds, 
the  type  of  input  likely  to  produce  extreme  floods  on  those  basins  must 
be  determined.  A  depth  of  excess  precipitation  (rain  or  snowmelt)  must 
also  be  selected.  The  frequency  with  which  the  selected  excess  preci¬ 
pitation  will  occur,  gives  an  indication  of  the  frequency  of  the  peak 
flow  determined  from  the  formulae. 

The  method  is  based  on  a  statistical  analysis  of  existing 
hydrometric  records  and  on  two  major  assumptions,  the  merits  and 
implications  of  which  will  be  discussed. 

The  method  assumes  that  there  is  a  limit  to  the  gradient 
which  can  be  attained  by  the  rising  and  falling  limbs  of  the  hydrograph 
of  a  particular  watershed.  It  is  generally  agreed  that  the  assumption 
is  approximately  true  for  the  falling  limb,  since  this  represents  the 
rate  of  withdrawal  of  water  from  storage  within  the  basin,  after  input 
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has  ceased.  The  slope  of  the  rising  limb,  however,  is  dependent  upon 
the  intensity  and  duration  of  the  input. 

Inspection  of  the  composite  hydrographs  showed  that  there 
were  sections  of  the  rising  limb,  having  a  nearly  constant  slope  over 
considerable  discharge  ranges.  This  slope  was  virtually  the  steepest 
to  occur,  and  indicated  that  a  limiting  slope  might  exist.  The 
occurrence  of  a  limiting  slope  over  a  wide  discharge  range  might 
indicate  that  the  maximum  possible  input  rate  occurred  or  was  approached 
quite  often,  for  short  durations. 

In  larger  basins,  where  channel  flow  is  predominant,  the 
ultimate  rate  of  rise  of  stream  flow  may  be  limited  by  the  increasing 
amounts  of  water  required  in  storage  to  produce  the  necessary  hydraulic 
head  to  move  the  water  through  the  basin,  combined  with  the  greater 
frictional  resistance  associated  with  overbank  flow.  Even  if  very  rare 
floods  have  rising  limbs  somewhat  steeper  than  that  indicated  by  the 
available  period  of  record,  it  can  be  seen  from  the  equations  of  Table 
8  that  peak  flow  is  not  highly  sensitive  to  variations  in  hydrograph 
slope . 

The  second  assumption  made,  is  that  the  limiting  rates  of  rise 
and  fall  of  the  hydrograph  limbs  are  virtually  straight  when  plotted  as 
logarithms  of  discharge  versus  time.  This  implies  that  the  discharge  on 
any  day  is  a  fixed  proportion  of  the  discharge  on  the  preceding  day.  The 
composite  hydrographs  produced  during  the  course  of  this  study  indicate 
that  this  is  approximately  true  for  considerable  ranges  of  discharge. 

It  is  interesting  to  compare  the  implications  of  these  assump¬ 
tions  with  those  of  unit  hydrograph  theory.  It  will  be  recalled  that  a 
unit  hydrograph  has  a  fixed  time  base,  i.e.,  the  time  to  peak  flow  and 
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the  time  for  recession  are  constant  for  a  fixed  input  duration,  regardless 
of  the  volume  of  excess  precipitation.  This  implies  that,  for  increasing 
runoff,  the  rates  of  rise  and  drop  also  increase.  However,  it  is  generally 
accepted  that  there  is  a  limit  to  the  recession  rate,  governed  by  the 
physical  characteristics  of  the  basin.  Thus,  if  the  volume  of  runoff  is 
taken  sufficiently  large,  the  recession  limb  becomes  steeper  than  is 
physically  possible.  It  seems  more  realistic  that  a  larger  volume  of 
water  in  storage  will  take  longer  to  drain,  once  input  has  ceased,  than 
will  a  smaller  volume,  so  that  the  recession  period  should  increase  with 
increasing  excess  precipitation. 

As  to  the  rising  limb,  neither  a  fixed  time  to  peak  flow  as 
in  unit  hydrograph  theory,  nor  a  constant  rate  of  rise,  can  be  justified 
theoretically.  The  implications  of  a  constant  rate  of  rise  are:  - 

a)  that  the  time  to  peak  flow  increases  as  the  volume 
of  runoff  increases, 

and  b)  that  at  a  particular  time  after  the  beginning  of 

runoff,  the  discharge  will  be  same  for  any  flood 
considered . 

With  reference  to  a)  above,  since  the  mean  value  of  the  rate 
of  rise  (R)  for  the  62  basins  studied  is  approximately  0.25,  the  discharge 
on  one  day  is  about  four  times  that  on  the  previous  day.  This  means 
that  for  increasing  runoff  volumes,  the  peak  discharge  increases 
rapidly  but  the  time  to  peak  increases  slowly.  A  slight  increase  in  the 
time  to  peak,  with  increasing  runoff,  has  been  noted  by  several  authors. 
(Linsley,  1958)  . 

With  reference  to  b)  above;  considering  that  the  procedure 
is  limited  to  'extreme*  flood  events,  the  concept  of  a  fixed  discharge 
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after  a  certain  time  is  not  unreasonable.  Possibly,  the  real  relation¬ 
ship  between  discharge  and  time,  during  the  rising  portion  of  the 
hydrograph,  lies  between  these  two  extremes.  However,  the  concept  of  a 
fixed  rate  of  rise  is  a  good  approximation  and  has  the  merit  of  simple 
mathematical  representation. 

An  interesting  feature  of  composite  hydrographs  produced  from 
existing  records,  was  the  occurrence  of  a  quite  definite  base  flow 
discharge.  It  appears  that,  at  the  cessation  of  direct  runoff,  river 
discharge  is  maintained  at  a  constant  level  due  to  the  influx  of  ground- 
water.  It  would  require  an  extremely  rare  drought  situation  to  affect 
the  baseflow  discharge.  For  small  basins  in  general,  and  for  those  of 
the  prairies  in  particular,  the  baseflow  discharges  are  zero. 

The  two  assumptions  discussed  previously,  together  with  the 
occurrence  of  a  distinct  baseflow  discharge,  lead  to  the  idealized 
hydrograph  of  Fig.  7  (triangular  when  plotted  to  semilogarithraic  scale) . 
A  feature  of  this  hydrograph,  and  a  source  of  error,  was  the  elimination 
of  a  crest-section.  Consideration  was  given  to  determining  a  crest 
width  from  the  composite  hydrographs  for  each  study  basin,  and  carrying 
out  a  regression  analysis  against  the  basin  characteristics.  However, 
the  crest  section  of  each  composite  hydrograph  was  derived  from  a  single 
storm  (see  Section  2.4)  and  was  not  necessarily  representative.  An 

/y 

analysis  of  many  crest  segments  for  each  basin  could  be  required  to 
accurately  assess  the  crest  shape.  Such  an  analysis  was  beyond  the 
scope  of  this  study. 

Regression  analysis  showed  that,  as  expected,  the  selected 
hydrograph  parameters  are  related  to  watershed  characteristics.  The 
characteristics  included  in  the  analysis  were  limited  to  those  for  which 
data  exists  or  can  readily  be  determined.  This  precluded  such  factors 
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as  soil  distribution  and  properties,  and  details  of  vegetation  cover. 

The  regression  analysis  on  the  hydrograph  shape  parameters  was  moderately 
productive;  regression  equations  were  produced  whose  multiple  correlation 
coefficients  ranged  from  0.72  to  0.86.  Later  analysis  showed  this  to  be 
sufficient . 

Equation  5.7  gives  a  relationship  between  peak  flow  (Q  )  and 
average  depth  of  excess  precipitation  (d)  rates  of  rise  and  fall  (R  and  K) 
baseflow  discharge  (Q^)  and  basin  area  (TDA) .  The  equation  is  non¬ 
linear  and  theoretically  will  produce  results  which  differ  from  the  unit 
hydrograph  concept  of  peak  flow  being  directly  proportional  to  depth  of 
excess  precipitation.  However,  for  a  particular  watershed,  the  rates  of 
rise  and  drop  are  assumed  constant  and  the  baseflow  discharge  is  very 
much  smaller  than  peak  flow  for  any  practical  depth  of  runoff.  Therefore, 
the  first  and  last  terms  of  Equation  5.7  become  insignificant  and  an 
almost  linear  relationship  between  peak  flow  and  runoff  depth  results. 
(Table  5  and  Fig.  8) . 

A  comparison  of  the  equations  in  Tables  5  and  6  shows  that 
basin  area  explains  a  large  portion  of  the  variation  between  peak  flow 
and  the  basin  characteristics.  This  is  because,  not  only  does  basin 
area  appear  directly  in  Equation  5.7,  but  it  contributes  significantly 
to  the  explained  variation  of  the  rates  of  rise  and  drop. 

It  might  appear  from  Table  6  that  basin  area  alone  is  suffi¬ 
cient  to  predict  peak  flow  accurately,  but  the  effects  of  hydrograph 
shape  cannot  be  ignored.  Consequently,  three  equations  have  been 
presented  (Table  8)  for  different  ranges  of  R  and  K,  allowing  the 
determination  of  peak  flow  from  basin  area,  for  a  given  depth  of  excess 
precipitation. 
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Analysis  of  individual  storms  over  the  Canadian  prairies 
(Gray,  1970)  yields  equations,  relating  peak  flow  with  basin  area,  of 


the  form 


Q  =  C  x  A 
o 


0.5 


where  C  ranges  from  300  to  2,300 

and  A  ranges  from  0.1  to  1,000  sq.  mis. 

A  comparison  of  the  above  with  Eq.  5.31  shows  that,  for  basins 

less  than  about  10  sq.  mis.,  the  formulae  derived  in  this  study  give  values 

of  peak  flow  less  than  those  experienced  on  the  prairies,  but  for  basins 

up  to  1,000  sq.  miles,  this  study  indicates  somewhat  higher  flows. 

A  general  formula  applicable  to  the  U.S.A.  (U.S.  Geological 

Survey  for  Columbia)  for  basins  between  1,000  and  24,000  sq.  miles  is 

Q  =  1400  a0*476 
o 

which  is  in  reasonable  agreement  with  Equation  5.32 

(Q  =  2,300  A*  )  for  basins  greater  than  10,000  sq.  miles. 


6 . 2  Procedure  to  Determine  Hydrograph  Shape  for  Ungauged  Watersheds 

Depending  on  the  amount  and  accuracy  of  the  information  required, 
the  results  of  this  study  can  be  used  to  determine  the  hydrograph  for  an 
ungauged  watershed  in  several  ways :  - 

1.  i)  Determine  required  depth  of  excess  precipitation  from  meteo¬ 
rological  records 

ii)  For  the  watershed,  determine  RAIN,  TDA,  FA,  MCL,  RE,  and 
MCS  according  to  the  definitions  of  Table  1. 
iii)  Determine  rates  of  rise  and  drop  (R  and  K)  and  baseflow 

discharge  (BF)  from  Equations  4.11,  4.12  and  4.13  using  the 
map  coefficients  appropriate  to  the  basin  location. 


iv) 


Substitute  the  above  values  in  equation  5.7  and  solve  for 


2. 


3. 


4. 


using  an  iterative  process, 
i)  to  iii)  as  above 

iv)  Determine  BA(Y)  and  solve  the  appropriate  equation  of 
Table  5  to  find  Q  ,  interpolating  as  necessary 
i)  to  iii)  as  above 

iv)  with  the  appropriate  value  of  R  and  K,  determine 

from  Fig.  12.  Having  determined  Q  ,  and  knowing  R  and  K, 
the  hydrograph  may  be  drawn  to  determine  the  time  distri¬ 
bution  of  runoff. 

If  only  the  peak  flow  is  required,  R  and  K  can  be  estimated 
from  a  knowledge  of  basin  characteristics,  and  determined 
from  Fig.  12  directly.  In  such  a  case,  R  and  K  need  only 
be  determined  as  high,  intermediate/mixed  or  low. 


6 . 3  Conclusions 

A  method  has  been  derived  to  determine  the  shape  of  a  hydro¬ 
graph  for  an  ungauged  watershed  in  the  Alberta/Yukon  region  of  Canada, 
for  a  given  depth  of  excess  precipitation.  The  slopes  of  the  rising 
and  falling  limbs  are  determined  from  the  physical  and  meteorological 
characteristics  of  the  watershed,  and  provide  two  numerical  parameters 
for  watershed  classification. 

The  assumption  of  constant  rates  of  rise  and  drop  was  shown 
to  be  reasonable  by  visual  inspection  of  hydrographs  derived  from  many 
years  of  flow  records.  This  assumption  implies  that  the  time  for  run¬ 
off  to  occur  increases  as  the  depth  of  excess  precipitation  increases 
and  differs  from  the  unit  hydrograph  theory  premise  of  a  constant  runoff 
time.  However,  the  unit  hydrograph  theory  assumption  that  peak  flow 
is  directly  proportional  to  excess  precipitation  was  substantiated. 
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For  small  watersheds,  meaningful  results  can  be  derived 
only  from  analysis  of  instantaneous  hydrographs;  too  much  detail 
being  lost  when  mean-of-the-day  discharges  are  computed. 

Only  a  moderate  linear  association  could  be  determined 
between  the  hydrograph  parameters  and  the  selected  watershed  characte¬ 
ristics  but  determination  of  map  coefficients,  to  reduce  the  effect  of 
variables  not  included  in  the  analysis,  improved  the  usefulness  of  the 
equations  for  predictive  purposes.  However,  the  most  important  para¬ 
meter,  peak  flow,  was  found  to  have  a  strong  linear  association  with 
drainage  basin  area.  The  regression  equations  allow  the  rates  of  rise 
and  fall  to  be  determined  with  sufficient  accuracy  for  peak  flow  values 
to  be  determined  for  a  range  of  hydrograph  shapes. 

6 . 4  Suggestions  for  Further  Study 

1.  If  sufficient  data  could  be  accumulated,  hydrograph  analysis  of 
watersheds  in  a  smaller  geographic  region  might  produce  stronger 
associations  with  basin  characteristics,  and  hence  increase  the 
predictive  value  of  the  regression  equations. 

2.  A  study  might  be  undertaken  to  correlate  flood  hydrograph  crest 
shape  with  basin  characteristics  to  make  the  triangular  hydrograph, 
used  in  this  analysis,  more  realistic. 

3.  A  comparison  might  be  made  of  the  values  of  R(av)  and  R(max) 

(Section  2.4)  to  determine  whether  they  approach  the  same  value, 
as  the  length  of  discharge  record  increases.  A  similar  analysis 
could  be  undertaken  for  K(av)  and  K(max) .  Such  an  analysis  could 
verify  the  assumption  of  a  limiting  slope  for  the  rising  and  falling 


hydrograph  limbs. 
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4.  The  computer  programs  developed  for  this  study  could  be  used  to 
determine  R  and  K  for  different  times  of  the  year.  Time  of  the 
year  could  then  be  introduced  into  the  regression  analysis  and  the 
resulting  equations  used  to  determine  peak  flows  at  different 
seasons . 
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APPENDIX  I 


COMPUTER  PROGRAMS 
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C  PROGRAM  11  SLOPE11 

C 

C  TO  COMPUTE  CO-ORDINATES  OF  THE  STEEPEST  POSSIBLE  HYDROGRAPH  FROM 

C  EXISTING  DISCHARGE  RECORDS  AT  ANY  GAUGING  STATION  WHERE  THE 

C  GIVEN  DISCHARGE  VALUES  ARE  SPACED  AT  REGULAR  TIME  INTERVALS. 

C 

C  Q=I N PU T  VECTOR  OF  DISCHARGE  VALUES. 

C  ALPH=HEADI NG  FOR  OUTPUT  (ALSO  USED  FOR  HYDROGRAPH  PLOT  TITLES.) 

C  X  AND  Y  =OUTPUT  VECTORS  OF  HYDROGRAPH  CO-ORDINATES. 

C 

REAL  N  SPA  N  ,  NGRAD 

DIMENSION  Q  (3  000) ,  Q PAIR  (2,3000) , QS  (2,3000) , 

♦NSPAN  (500)  , NGRAD  (5  00)  ,  ALPH  (20)  ,  Y  (3000)  ,R  (500)  ,D(500)  , 

♦DROPS  (3000)  , RISE  (3000)  ,Y2  (3000)  , QP  (2 , 3000 ) , X  (3000 ) 

C 

c 

C  READ  DISCHARGE  DATA  (MONTHS  N1-N2  INCL.  ONLY) 

C  NM  IS  NUMBER  OF  MONTHS  OF  DATA  ON  FILE 

C  IF  DATA  READ  FROM  FILE  WRITTEN  BY  "ADD",  ADD=1.0 

READ  (5,10)  NM,N  1 , N2 ,  ADD 

10  FORMAT  (316  ,F6 .0) 

NST ART= 1 

WRITE  (3,1)  NSTART 
1  FORMAT  (IX,  16) 

JP9  =  0 

JP10=0 

J=2 

Q  ( 1 ) =-999. 

IF  (ADD  .EQ.  1.0)  GO  TO  3 
DO  25  N= 1 , NM 
DO  26  11=1  ,2 
JP9=J+9 

READ (4 ,11) L,M,  (Q  (I) ,I  =  J,JP9) 

11  FORMAT  (11X,I2,I1,  (10F6.0) ) 

IF  (N.EQ.1  .AND.  II  .  EQ.  1  .AND.  M  .NE.  1)  GO  TO  13 
GO  TO  15 


13  WRITE  (6,14) 

14  FORMAT  ('  '  , ' DATA  DOES  NOT  START  AT  BEGINNING  OF  MONTH') 
GO  TO  970 


15 

IF ( N 1  .LE.  N  2  .AND. 

(L  .LT. 

N 1  .OR. 

L  .GT. 

N2)  ) 

GO 

TO 

26 

I F ( N 1  .GT.  N2  .AND. 

J  =  JP9+  1 

(L  .LT. 

N 1  .AND. 

L  .GT. 

N2)  ) 

GO 

TO 

26 

26 

CONTINUE 

JP10=J+10 

READ  (4,12)  L ,  (Q  (I)  , 

I=J,  JP10) 

12 

FORMAT  (11X,I2,1X,  (11F6.0)  ) 

IF  ( N 1  .LE.  N  2  .AND. 

(L  .LT. 

N 1  .OR. 

L  .GT. 

N2)  ) 

GO 

TO 

24 

IF  ( N 1  .GT.  N2  .AND. 

(L  .LT. 

N 1  .AND. 

L  .GT. 

N2)  ) 

GO 

TO 

24 

J=JP10  +  1 
GO  TO  25 

24  Q  ( J)  =-9  99. 

J=J«-1 

25  CONTINUE 

Q  (J)  =-999. 

Q  (U+1) =-999. 

GO  TO  4 

3  READ  (4  ,16)  K 

16  FORMAT  (17) 

J=K+1 

REA  D (4 , 20)  (Q  (I)  ,1  =  2, J) 

20  FORMAT  (1X,10F7.0) 

Q(J-H)  =-999. 

J=J+1 

4  IF  (J-3000)  2,2,8 

8  WRITE  (6,9) 

9  FORMAT  ('  ' , ' NO.  OF  DATA  POINTS  EXCEEDS  3000.') 

GO  TO  970 

C 

C  WRITE  HEADING 

2  READ  (5,5)  (ALPH  (I)  ,  1=1 , 20) 

5  FORMAT  (20A4) 

WRITE  (6,6)  (ALPH(I)  ,1=1,  12) 
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7  4 
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94 

95 

96 

97 

98 

99 

ICO 

101 

102 

103 

104 

105 

106 

107 

108 

109 

no 

in 

112 

113 

1 1  4 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

1  33 
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138 

139 

1  40 

141 

1  42 

143 
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6  FORMAT  (1H  1 ,6X,20A4) 

WRITE  (3,7)  (ALPH(I)  ,1=1,20) 

7  FORMAT  ( •  1  ,  20 A4) 

C 

C  TAKE  COMMON  LOGS  OF  DISCHARGE  VALUES 

DO  46  1=1 , J 
IF  (Q  ( I)  )  46,45,  45 

45  IF  (Q  (I)  .  LE.  1 .0)  Q  (I)  =1.001 
Q  (I )  =  ALOG  1 0  (Q  (I)  ) 

46  CONTINUE 

21  FORMAT  (»  »,F10.4) 

C 

C  DELETE  MISSING  DATA 

C  LIST  PAIRS  OF  DISCHARGES  ON  CONSECUTIVE  DAYS 

M=0 

DO  70  1=2, J 
IF  (Q  ( I) )  60,50,  50 
50  M  =  M  + 1 

QP(1,H)=Q(I) 

QP (2 , M ) =Q  (1+1) 

GO  TO  70 

60  IF  (Q  (1-1)  .GE.O. 0)  M=M-1 

70  CONTINUE 

DO  75  1=1, M 
QS  (1,I)=QP  (1,1) 

QS(2,I)=QP  (2,1) 

75  CONTINUE 

WRITE  (6,36)  M 

36  FORMAT  ('  '//6X,'NO.  OF  PAIRS  OF  CONSEC.  DISCHARGE  VALUES=«,I6) 
125  FORMAT  ( *  ' //1 6 , 2  FI  0 . 4 ) 

C 

C  SEPARATE  PAIRS  INTO  +VE  (QP)  AND  -VE  (QPAIR)  SLOPES 

JJ=0 
MM  =  0 

DO  99  1=1, M 

IF  (QS  (2,1)  .GT.  QS  (1 ,1)  )  GO  TO  98 
JJ  =  JJ+  1 

QPAIR  ( 1  ,JJ) =QS  (1,1) 

QPAIR (2, JJ) =QS (2,1) 

GO  TO  99 

98  MM  =  MM+  1 

QP(1  ,MM)  =  QS  (1  ,1) 

QP  (2, MM)  =  QS  (2,1) 

99  CONTINUE 
C 

C  FIND  LOCAL  MINIMAS  8  DISCONTINUITIES 

N=0 
N2  =  0 
N  3=0 

DO  100  1=2, M 
A=QS(1,I) 

B=QS  (2,1) 

C=QS (1,1-1) 

E=Q  S  (2,1-1) 

IF  (  (A  .  EQ.  E)  .  AND.  (B.GT.  A)  .  AND.  (C.GT.  A)  )  GOTO  101 
IF  (  (A  .  NE.  E)  .  AND.  (B.GE.  A)  )  GO  TO  102 
IF  (  (A.  NE.  E)  .  AND.  (B.  LE.  A)  )  GO  TO  103 
GO  TO  100 

101  N=  N  +  1 
Y2  (N)  =  A 
GO  TO  100 

102  N2  =  N2+  1 
R  (N  2)  =  A 
GO  TO  100 

103  N3=N3+ 1 
D  ( N  3 )  =  E 

100  CONTINUE 
C 

C  RECESSION  LIMB 

C 

C  ORDER  QPAIR  ON  1ST  ROW 

CALL  SORTRO  (QPAIR, 2, JJ, 1) 
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c 

C  DETERMINE  POINTS  TO  BE  USED  IN  SEARCH 

DO  108  1=1, N 
QS(1,I)=Y2  (I) 

QS  ( 2, 1 )  =  Y  2  (I) 

108  CONTINUE 

DO  1 C  4  1=1, JJ 

QS  (1  ,N  +  I) = QPAIR  (1,1) 

QS  (2,  N+I)  =  QP A IR  (2, 1) 

104  CONTINUE 
NI=N*I 

IF  (N3)  105,105,107 
107  DO  105  J=  1 , N 3 

QS(1,NI+J) =D(J) 

QS ( 2 , N 1+ J)  =D  (J) 

105  CONTINUE 
NJ  =  NI  +  N  3 

CALL  SORTRO  (QS,2 ,NJ,1 ) 

C 

C  COMPUTE  RANGES  SPANNING  EACH  DISCHARGE  VALUE 

L1  =  0 

DO  200  1=1 ,NJ 
K=0 

DO  150  J=1,JJ 
Q1  J  =  QP  AIR (  1, J) 

Q2J=QP AIR  (2, J) 

Q1I  =  QS  (1,1) 

Q2I=QS  (2,1) 

IF  (Q1J.LT. Q1I)  GO  TO  129 
IF  (Q1J.GE.Q1I.AND.Q2J.LT.Q1I)  GOTO  160 

128  FORMAT  (216) 

GO  TO  150 

160  K  =  K  + 1 

NSPAN (K) =Q  1J 
NGRAD (K)=Q2J-Q1 J 
150  CONTINUE 
C 

C  SELECT  STEEPEST  GRAD.  FOLLOWING  DISCHARGE  VALUE 

129  IF  (K)  130,200,130 

13C  GRI  =  Q2 1-Q 1 1 

TE  M  P  =  G  PI 
L 1  =  L  1  +  1 
DO  2 1 C  J=1  ,K 

IF  (NGRAD  (J)  .LE.  TEMP)  GO  TO  300 
GO  TO  210 
300  KK= J 

TEMP  =  NGRAD  (J) 

210  CONTINUE 
C 

C  LIST  DESCENDING  DISCHARGE  VALUE  S&  STEEPEST  FOLLOWING  GRAD. 

IF  (LI  .EQ.  1)  GO  TO  250 
IF  (TEMP-X  (LI-1)  )  220,230,240 
220  X (L 1 ) =  TEM  P 

Y  (LI)  =  NSP A  N  (KK) 

GO  TO  200 

230  Y  (LI)  =  Y  (LI-1) 

X  (L  1 )  =  X  (LI-1) 

GO  TO  200 
240  X  (L  1 )  =  TEMP 

Y  (LI)  =  Q1I 
GO  TO  200 

250  Y  ( 1 ) =Q  1 1 
X  (1)  =TEMP 
200  CONTINUE 
C 

C  CONDENSE  LIST  OF  ORDINATES 

Q(1)=Y  (1) 

DROPS  (  1)  =X  (1) 

L2=  1 

DO  310  1=2, LI 

IF  (Y  (I)  .  EQ.Q  (L2)  )  GO  TO  310 
L2  =  L2+  1 
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Q(L2)=Y(I) 

DROPS  (L2)  =X  (I) 

310  CONTINUE 
131  FORMAT  (F10 . 4) 

WRITE  (6,3  15)  L 2 

315  FORMAT  ( '  • //6 X, • ORIGINAL  NO.  OF  POINTS  ON  RECESSION  LIMB=',I6) 
C 

C  RISING  LIMB 

C  ORDER  QP  ON  2ND  F.OW 

CALL  SORTRO  (QP,  2  , MM, 2) 

C 

C  DETERMINE  POINTS  TO  BE  USED  IN  SEARCH 

DO  405  1=1  ,  N 
QS  (1,1)  =Y2  (I) 

QS  (  2 ,  I )  =Y  2  (I) 

405  CONTINUE 

DO  400  1=1, MM 
QS  (1  ,N+I)  =  QP  (1 ,1) 

QS  (2,N  +  I)  =QP  (2,1) 

400  CONTINUE 
NI=N*I 

IF  (N2)  401,401,402 
402  DO  401  J=1  ,N2 

QS (1,NI+J)  =R  (J) 

QS (2 , NT  +  J)  =R  (J) 

401  CONTINUE 
NJ=NI+  N2 

CALL  SORTRO  (QS,2,NJ,2) 

C 

C  COMPUTE  RANGES  SPANNING  EACH  DISCHARGE  VALUE 

LI  =  0 

DO  500  1=1, NJ 
K=0 

DO  650  J=1,MM 
QP2J  =  QP  (2  ,  J) 

QP1  J=QP  (1 ,  J) 

QP1T=Q  S  (1 , 1) 

QP2I  =  Q S  (2,  I) 

IF  (QP 2 J. LT. QP2I )  GO  TO  653 

IF  (QP2J. GE.QP2I. AND. QP1J.LT. QP2I)  GOTO  560 
GO  TO  650 
560  K=K  ♦  1 

NSPAN  (K)  =QP2J 
NGRAD ( K ) =QP2J-QP1 J 

650  CONTINUE 

651  FORMAT  (16  ,  (2F10. 4)  ) 

C 

C  SELECT  STEEPEST  GRAD.  PRECEEDING  DISCHARGE  VALUE 

653  IF  (K)  652,500,652 

652  GRI=QP2I-QP1I 
TEMP=GRI 

L 1 =L 1 +  1 
DO  710  J= 1 , K 

IF  (NGRAD  (J)  .  GE.  TEMP)  GO  TO  800 
GO  TO  710 
800  KK= J 


710 

C 

C 

740 

720 

760 

730 


TEMP=NGRAD  (J) 
CONTINUE 


LIST  DESCENDING 
IF  (LI  .  EQ.  1)  GO 
IF  (TEMP-X  (LI- 1) 
X  (LI)  =TEMP 
Y  (LI)  =  NSP  A  N  (KK) 
GO  TO  500 
X  (LI)  =  TEM P 
Y (LI)  =  QP2 1 
GO  TO  500 
Y (LI)  =  Y  (LI  -1 ) 

X  (LI)  =X  (LI-1) 

GO  TO  500 
Y (1)  =Q  P2I 


DISCHARGE  VALUES  S  STEEPEST  PRECEEDING  GRAD. 
TO  730 

) 720, 760, 740 
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X  (1)  =  TEMP 
500  CONTINUE 
C 

C  CONDENSE  LIST  OF  ORDINATES 

Y2  (1 )  =  Y  (1 ) 

RISE  (1)  =X  ( 1) 

L3=  1 

DO  750  1  =  2  , LI 

IF  (Y  (I)  .  EQ.Y2  (L3)  )  GO  TO  750 
L3  =  L3+  1 
Y2  ( L 3 )  = Y (I) 

RISE  (L3)  =X  (I) 

750  CONTINUE 

WRITE  (6,755)  L3 

755  FORMAT  ('  '//6X, • ORIGINAL  NO  OF  POINTS  ON  RISING  LIMB=',I6) 

C 

C  PREPARE  DATA  FOR  PLOTTING 

X(1)=0.0 

Y  (1)  =  Y  2  (L 3)  -RISE  (L3) 

L=  1 

DO  810  1=1, L3 
L=L+  1 

Y  (L)  =Y2  (L3-I+1) 

X  (L)  =X  (L-1  )♦  (  (Y  (L)  -  Y  (L-1)  )  /RISE  (L3-H-1)  ) 

810  CONTINUE 

Y (L+1) =Q  ( 1 ) 

IF  (Y(L)-Q(I))  820,830,  840 
8  20  X  (L+1)  =X(L)  *■  (  (Y  (L+1)  -Y  (L)  )/RISE  (1)  ) 

GO  TO  850 
830  X(L+1)=X(L) 

GO  TO  850 

84  0  X  (L+1 )  =X(L)+  {  (Y  (L+  1)  -Y  (L)  )  /DROPS  (1)  ) 

850  L4=L+ 1 

DO  860  1  =  2  , L2 
L4=L4+  1 
Y(L4)=Q(I) 

X  (L4)  =X  (L4-1)  ♦  (  (Y  (L4)  -Y  (L4-1)  )  /DROPS  (1-1)  ) 

860  CONTINUE 

Y  (L4  +  1)  =Y  (L4)  +  DROPS  (L2) 

X (L4+ 1 ) =X ( L4 ) +1 . 0 

L5  =  L  4+  1 

WRITE  (6,891)  L5 

891  FORMAT  (•  • //6X, • ORIGINAL  NO.  OF  DATA  POINTS=«,I6) 

C 

C  FIND  RANGE  OF  LOG-DISCHARGE  VALUES 

R  (1)  =  A  MIN  1  (Y  (1 )  ,  Y  (L5)  ) 

R  (2)  =  AM  AX  1  (Q(1)  ,  Y  2  ( 1 )  ) 

HITE  =  R  (2)  -R  (1 ) 

WRITE  (6,55)  R  (2)  ,R  (1)  ,  HITE 

55  FORMAT  (’  '//6X, 'RANGE  OF  LOG-DISCHARGE  V ALU ES= '  ,  F7 . 4  ,  •  MINUS', F7. 4 
*,'=', F7. 4) 

D (1) =AINT  (R  ( 1 ) ) 

D  (2)  =AINT  (R  (2)  ) +1. 0 
WRITE(6,74  6)  D  ( 1 )  ,  D  (2) 

7U6  FORMAT  (•  • //6 X , ' OR DINATE  RANGE  IS  FROM  10**',F4.1,'  TO  10**',F4.1) 
WRITE  (3,747)  D  ( 1 )  ,  D  ( 2) 

747  FORMAT  ('  •  ,2F 10. 4) 

C 

C  REDUCE  NO.  OF  PLOT  POINTS  TO  500 

DIST=  X  (L5) -X  (1) 

D£N= 10C  0. 0 
N=1 
L6  =  L5 
GO  TO  895 
896  L6=N 
N  =  1 

895  IF  (L6-500)  950,  950,  910 
910  Q  ( 1 )  =  Y  (1) 

Y2  ( 1)  =  X  ( 1) 

DEN  =  DE  N/2 . 0 
DO  920  1=2, L6 

IF  (ABS  (Y2  (N)  -X  (I)  )  .  LE.  DIST/DEN.  AND.  ABS  (Q  (N)  -Y  (I)  )  .  LE .  HITE/DEN) 
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*GO  TO  920 
N=  N  1 

Y 2  (N)  =  X  (I) 

Q(N)=Y  (") 

920  CONTINUE 

921  FORMAT  ( 1  ',16) 

DO  940  1=1 , N 
X(I)=Y2(I) 

Y  (I)  =Q  (I) 

940  CONTINUE 

IF  (N. GT. 500)  GO  TO  896 
L5  =  N 

950  CONTINUE 
C 

C  FOR  DATA  PREPARED  BY  "ADD",  CONVERT  BACK  TO  WHOLE  DAYS 

IF  (ADD  .EQ.  0.0)  GO  TO  951 
DO  952  1=1, L5 
X  (I )  =X  (I)  /2.0 
952  CONTINUE 
C 

C  TAKE  ANTI-LOGS  OF  DISCHARGE  TO  PLOT  ON  SEMI-LOG  SCALE 

951  DO  870  1=1  ,L5 

Y  (I)  =(  10.0**Y  (I)  )  +0 . 0  2 
870  CONTINUE 

C 

WRITE  (3,960)  L5 
KRITE(6,890)  L5 

890  FORMAT  (1H  //6X,5 NUMBER  OF  DATA  POINTS  TO  BE  PLOTTED= •  , 1 6) 
880  FORMAT  (1H  //I  OX,  •  TIME  (DAYS)  DISCHARGE  (CFS)  ' ) 

WRITE  (6,880) 

WRITE  (  6,96  0)  (  (I  ,  X  (I )  ,  Y  (I)  )  ,  1=  1 ,  L5) 

WRITE  (3,960)  (  (I,X  (I)  ,  Y  (I)  )  ,1  =  1  ,L5) 

960  FORMAT  (‘  '  ,1  6  ,  F  1 2 . 4  ,  FI  4 . 0) 

965  FORMAT  (■  *,I6,2F12.4) 

GO  TO  975 
970  STOP  4 
975  STOP 
END 


1 

C 

PROGRAM  "ADD" 

2 

C 

3 

c 

TO  FABRICATE  A  DISCHARGE  RECORD  BY  INTERPOLATING  DISCHARGE  VALU 

4 

c 

AT  HALF-DAY  INTERVALS  FROM  THE  CO-ORDINATES  OF 

UP  TO 

FIVE  SUB- 

5 

c 

PERIOD  HYDROGRAPHS. 

6 

c 

7 

c 

X  AND  Y  =  SUB-PERIOD  HYDROGRAPH  CO-ORDINATES. 

8 

c 

QQ=OUT  PUT  DISCHARGE  VALUES  INTERPOLATED  AT  HALF 

-DAY 

INTERVALS. 

9 

c 

10 

DIMENSION  X  (5  00)  , Y  (50 0) , QQ  (1000)  , Q1  (500) 

1 1 

5 

FORMAT  (///I7) 

12 

10 

FORMAT  (7X, FI  2. 4, FI  4.0) 

13 

READ  (8  ,  1 2)  N 

14 

12 

FORMAT  (24  X ,16) 

15 

J=0 

16 

DO  40  11=1  ,  N 

17 

GO  TO  (6, 7, 8, 9,  13)  ,  II 

18 

6 

READ  (1  ,5)  N1 

19 

READ  (1  ,  10)  (X  (I)  fY  (I)  ,1=1, N1) 

20 

70 

FORMAT  ('  '  ,2F10. 4) 

21 

GO  TO  15 

22 

7 

READ  (2, 5)  N1 

23 

READ  (2,10)  (X  (I)  ,Y  (I)  ,1=1, N1) 

24 

GO  TO  15 

25 

8 

PEAD(3,5)  N 1 

26 

READ  (3 ,10)  (X (I)  ,Y  (I)  ,1=1 , N1) 

27 

GO  TO  15 

28 

9 

READ  (4 , 5)  N  1 

29 

RE  AD ( 4  ,  1 0 )  (X  (I)  ,Y  (I)  ,1=1  ,  N1) 

30 

GO  TO  15 

31 

13 

READ  (9,5)  N1 

32 

READ  (9 , 10)  (X (I)  ,Y  (I)  ,1  =  1 ,N1) 

33 

C 

34 

c 

INTERPOLATE  VALUES  AT  HALF-DAY  INTERVALS 

35 

15 

DO  20  1=1,  N1 

36 

X  (I)  =2.0*X  (I) 

37 

20 

CONTINUE 

38 

ND  =  INT  (X  (N  1)  )  ♦  1 

39 

IF  (ND  .LT.500)  GO  TO  24 

40 

WRITE  (6,21)  II 

4  1 

21 

FORMAT  ('  • /2CX, ' SUB-PERIOD  NO.  ',12,'  :NO.  OF 

DATA 

PTS.  >500') 

42 

GO  TO  60 

43 

24 

CALL  DAY  (X ,Y,Q1 , ND,K) 

44 

WRITE  (6,25) II, K 

45 

25 

FORMAT  ('  ' /20X, *  SUB-PERIOD  NO.  ', 12 ,'  :NO.  OF 

DATA 

PTS .  = '  ,13) 

46 

DO  30  I:=  1  ,  K 

47 

J=J+1 

48 

IF  (J  .LE.  999)  GO  TO  29 

49 

WRITE  (6,27) 

50 

27 

FORMAT  (•  '/20X, 'TOTAL  NO.  OF  DATA  PTS.  >1000') 

51 

GO  TO  60 

52 

29 

CQ(J)=Q1(I) 

53 

30 

CONTINUE 

54 

J  =  J  +  1 

55 

QQ  (J)  =-999. 

56 

40 

CONTINUE 

57 

WRITE (6,26)  J 

58 

26 

FORMAT  ('  '/20X, 'TOTAL  NO.  OF  DATA  PTS.  =  ',I4) 

59 

WRITE  (7,45)  J 

60 

45 

FORMAT  ('  '  ,16) 

61 

WRITE  (7,50)  (QQ  (I)  ,1=1,  J) 

62 

50 

FORMAT  ('  •  ,  10F7. 0) 

63 

GO  TO  65 

64 

60 

STOP  4 

65 

65 

STOP 

66 

END 

STOP  0 

12:09.15  .571  RC=0 
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C  SUBROUTINE  "SUBDAY" 

C 

C  TO  INTERPOLATE  HYDROGRAPH  CO-ORDINATES  AT  HALF-DAY  INTERVALS. 

C 

SUBROUTINE  D A Y  (X  ,  Y , Q , ND , K ) 

DIMENSION  X  (500)  ,Y  (500) ,Q  (500) 

Q(1)=Y  (1) 

A  =  0 . 0 
K=  1 
J=1 

DO  10  I  =  2f  ND 
A  =  A+  1 . 0 

14  IF  (A-X  (J)  )  20,20,15 

15  J=J+1 

GO  TO  14 
20  K=K+ 1 

Q  (K)  =Y  (J-  1 )  ♦  (  (A-X  (J-  1)  )  *  (Y  (J)  -Y  (J-1)  )  /(X(J)  -X  ( J-1)  )  ) 

10  CONTINUE 

RETURN 
END 
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1 

C 

SUBROUTINE  "SORTROW" 

2 

C 

TO  REARRANGE  ALL  THE  COLUMNS  OF  A  TWO- 

3 

c 

ARRANGING  A  SELECTED  ROW  IN  DESCENDING 

4 

c 

5 

c 

X  =  IN PUT  ARRAY  TO  BE  REARRANGED  : NO.  OF 

6 

c 

NO.  OF 

7 

c 

NROW=NO.  OF  THE  ROW  TO  BE  ARRANGED  IN 

8 

SUBROUTINE  SORTRO (X , IROW , ICOL,NROW) 

9 

DIMENSION  X(2 ,5000)  ,TEM  (2) 

10 

L=0 

1 1 

DO  300  1= 1 ,ICOL 

12 

L=L  1 

13 

TEM  P=X  (NRO  W , I ) 

14 

DO  200  J=L , ICOL 

15 

IF  (X  (NROW  ,J)  .GE.TEMP)  GO  TO  100 

16 

GO  TO  200 

17 

100 

TEM  P  =  X  (NRO  W,  J) 

18 

M= J 

19 

200 

CONTINUE 

20 

250 

DO  300  JJ  =  1 , IROW 

21 

TEM  (JJ)  =X  ( J  J ,  M) 

22 

X  (J  J ,  M )  =X  (JJ,  I) 

23 

X  (JJ,I )  =  TE  M  (JJ) 

24 

300 

CONTINUE 

25 

RETURN 

26 

END 
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C  PROGRAM  "DRAW" 

C 

C  TO  PLOT  HYDROGRAPH  ON  S EM  I- LOG R I TH M I C  AXES  TO  SCALES  OF:- 

C  5 . 2 1  NS /  LOG  CYCLE  FOR  DISCHARGE, 

C  1.0INS/  DAY  FOR  TIME. 

C  MAX.  PLOT  LENGTH  =3  6  INS. 

C 

DIMENSION  ALPH  (20)  ,R(2)  ,E  (2)  ,X(500)  ,Y  (500) 

READ  (3,6)  NSTART 

READ  (3 ,3)  (ALPH  (I)  ,1=1, 20) 

3  FORMAT  (IX, 20A4) 

READ  (3 , 5)  D  ( 1 )  ,  D  (2) 

5  FORMAT  (1X,2F10. 4) 

R E  A D  ( 3 , 6)  L5 

6  FORMAT  (IX,  16) 

READ  (3 , 20)  (  (X  (I)  ,  Y  (I)  )  ,1  =  1 , L 5) 

20  FORMAT  (7  X  ,F  1  2 . 4  ,  F  1  4 .  0 ) 

J=  N ST A  RT-  1 
C 

C  REDUCE  Y-VALUE  RANGE  TO  5  LOG  CYCLES 

IF  (  (D  (2) -D  (1 )  )  .  LE.  5.0)  GO  TO  24 
D  ( 1 )  =  D  (2)  -5.0 
K  =  0 

DO  30  I=NSTART,L5 
IF  (Y  (I)-10.0**D  (1)  )  30,32,32 
32  K=K  +  1 

Y(K)=Y  (I) 

X(K)=X  (I) 

30  CONTINUE 

L5  =  K 

WRITE  (6,34) 

34  FORMAT  ('  V///20X, 'ORIGINAL  Y-AXIS  LENGTH  EXCEEDS  261  NS : '/2 1 X , • Y- A 

*XIS  REDUCED  TO  5  LOG-CYCLES') 

GO  TO  55 

24  L5  =  L  5- J 

IF  (NSTART. EQ. 1 )  GO  TO  55 
DO  50  1=1, L5 
J=J+  1 
X(I)=X  (J) 

Y  (I)  =Y  (J) 

50  CONTINUE 

C 

C  REDUCE  TIME  SPAN  TO  36  DAYS:  DELETE  LATER  VALUES 

55  DO  40  1=1, L5 

IF  (  (X  (I) -X(1)  )  .  GT.  36.0)  GO  TO  45 
40  CONTINUE 

GO  TO  25 
45  L5=I-1 

25  HC=  (X  (L5)  -X  (  1)  )  ♦  1.  0 
VB  =  D (2) -D  (1) 

VC=5.  2  * V B 
H A  =  A  I  N  T  (X  ( 1)  ) 

CALL  CGPL(X,Y,Y,L5,-1,1,2,4,1,HA,1.0,HC,D(1)  ,VB,VC,ALPH,6) 

STOP 

END 
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vO 

o 

o 

O 

CO 

o 

1 - 1 

CM 

CM 

o 

CO 

0) 

60  ^ 
CO  > 
H  c0 
a)  V— r 

4  K 


CO 

r-'- 

in 

CM 

CO 

CM 

O 

CN 

CN 

<T 

o 

\ — 1 

CM 

CM 

in 

in 

r-. 

CM 

rH 

00 

00 

vO 

in 

vo 

CM 

rH 

*Cf 

i — i 

in 

in 

co 

CO 

o 

in 

<1) 

px 

o 

iH 

CO 

'i 

•H 

rj 

60 

C 

•H 

CO 

•H 


•  • 

4-J 

CD  , — s 
<D  X 

VO 

1 — 1 

00 

rH 

ON 

On 

CM 

vO 

in 

CO 

o 

00 

CO 

00 

CO 

in 

vO 

in 

o 

O' 

in 

i — i 

cx  p 

ON 

oo 

o 

CM 

<r 

o 

CO 

vO 

o 

o 

o 

o 

<D  B 

C- 

o 

o 

o 

VO 

o 

NT 

CO 

n 

o 

o 

o 

o 

d)  NH 

1 — 1 

o 

o 

o 

o 

o 

t — 1 

o 

CM 

o 

o 

o 

o 

X)  psd 

oo 

<D 

60  ^ 

VO 

in 

in 

CM 

cd 

> 

vo 

On 

in 

vO 

<1* 

vo 

o 

r-'. 

co 

LO 

cd 

NT 

in 

t — 1 

o 

CO 

CO 

in 

CO 

vo 

^d- 

o 

0)  V_/ 

CO 

o 

o 

\ — 1 

1 — 1 

o 

rH 

1 — 1 

CM 

o 

O 

o 

o 

>  C3 
<3 

O 

CO 


G 

o 

•H 

4J 

cO 

4-> 

CO 


o 

53 


co 

CM 

co 

00 

On 

in 

ON 

o 

t — i 

< — i 

CM 

in 

CM 

O 

o 

o 

i — 1 

o 

o 

CM 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

<3 

PQ 

u 

q 

Q 

frj 

rJ 

rQ 

S3 

PQ 

u 

w 

w 

< 

<3 

C 

<3 

<d 

PQ 

PQ 

u 

u 

o 

a 

u 

in 

in 

in 

in 

LO 

in 

in 

in 

in 

in 

in 

in 

in 

O 

53 


CM 


co 


m 


VO 


00 


On 


CNI 


CO 
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co 

tH 

o 

CO 

CP 

i — 1 

CM 

CO 

o 

CD 

CO 

CO 

CO 

CD 

r^- 

p 

CP 

<r 

O 

1 - 1 

cD 

m 

in 

CM 

rH 

1 — 1 

CM 

CP 

CM 

<T 

P 

P 

cp 

m 

re 

re 

re 

re 

re 

re 

re 

re 

re 

re 

re 

re 

cP 

CO 

oo 

<T 

CM 

co 

00 

O 

r- 

r- 

in 

m 

p 

1 - 1 

CM 

CM 

CO 

rH 

CM 

P 

X — 1 

CM 

CO 

CM 

rH 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

m 

m 

i — ! 

<t 

o 

CM 

rH 

CM 

CM 

- 

w 

CD 

CO 

CM 

co 

c n 

<T 

CO 

CO 

X — 1 

r- 

CM 

co 

co 

CO 

co 

CO 

cP 

o 

o 

o 

in 

CP 

1 - 1 

CO 

p 

in 

CO 

CM 

CM 

CO 

o 

CM 

o 

o 

x — 1 

X — 1 

CP 

CM 

o 

CM 

r". 

1 - 1 

p 

CP 

<t 

r- 

x — 1 

00 

in 

CP 

P 

CO 

CM 

CO 

CD 

co 

p 

p 

p 

r-- 

CO 

<t 

co 

CM 

<fr 

f"> 

in 

x — 1 

• 

CO 

00 

CO 

1 — 1 

CO 

1 - 1 

m 

CM 

CP 

CM 

n- 

CD 

X — 1 

O 

r- 

in 

o 

o 

<d" 

o 

o 

O 

p 

CM 

i-'- 

p 

m 

o 

o 

o 

o 

o 

o 

in 

P 

o 

CD 

X — 1 

• 

co 

C 

o 

o 

o 

o 

o 

CM 

m 

O 

1 — 1 

O 

r-~ 

CM 

CM 

1 — I 

<d" 

m 

CM 

CP 

CP 

CM 

m 

m 

r'- 

P 

00 

m 

CO 

m 

CO 

x — 1 

P 

o 

p 

o 

O 

C M 

CO 

o 

1 — 1 

1 — 1 

o 

o 

o 

CD 

CO 

rH 

CO 

• 

• 

rH 

CN 

rH 

x — 1 

X - 1 

i — 1 

1 — 1 

CO 

CP 

x — 1 

X - 1 

CM 

rH 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

PQ 

o 

<3 

w 

<3 

p 

w 

<3 

p 

pq 

p 

Eh 

P 

p 

w 

w 

PH 

fH 

IH 

o 

< 

<3 

<3 

<3 

<3 

m 

in 

m 

in 

m 

in 

in 

in 

CD 

r- 

m 

VP 

CO 

p 

o 

1 — 1 

CM 

co 

m 

cp 

\ — i 

rH 

1 — 1 

x — | 

rH 

1 — 1 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

£ 

o 

i-i 

IX  E 
o 
g 

G  4-1 

0) 

P4 


W 

CO 

<u  = 

O  rH 
X  ^ 
w  o 

O' 

• 

CO  •• 

C  4-i 
*H  (X 
P 


P 

P 

a) 

CO  •• 
cc3  £ 
pq  o 


0) 

tx 

o 

I- 1 

co 

'E 

•H 

P 

00 
G 
•H 
rH 
r — i 

cc) 

P 


4-1 

CO  /*  s 

<U  X 

CX  G 

<D  E 

d) 

+j  i4 

co 


0) 

00 

CO  > 
5-i  cO 
a)  w 
>  t*3 
<3 


a) 

Pi 

o 

tH 

CO 

H 

•H 

P 

00 

G 

•H 

CO 

•H 

Pd 


4-1 

CO  *->, 

0)  X 
CX  CO 

a)  E 

<D 

4-t  pd 
CO 


CL) 

00  N 

G  > 
G  cO 
d) 

>  P 
<3 


G 

•  o 

C_>  P  • 

•  4-1  O 
CO  CO  53 

•  4-1 
ts  CO 


o 

33 
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55 

o 

t — ! 

PM 

6 

Excess 

o 

m 

on 

783 

co 

00 

On 

266 

148 

<3 

VO 

rH 

114 

1 - 1 

VO 

C> 

ON 

o 

vO 

vO 

m 

vO 

245 

oo 

o 

00 

o 

• 

0\ 

#v 

tr\ 

r* 

r\ 

€\ 

r\ 

r» 

e» 

ri*5 

03  •• 

rH 

on 

n- 

LO 

i — 1 

!3 

CO 

ID 

rH 

00 

m 

1 — 1 

cd 

0) 

PH 

<44 

.  in 
Ppt 

CM 

m 

tH 

CM 

1 — 1 

O 

tH 

CO 

lO 

1 — 1 

Pq 

FQ 

(l) 

C/D  •• 

cd  £ 
FQ  O 


O 

CM 


o 

o 

o 

o 

CM 

1 — 1 

o 

m 

O 

o 


o 


m 

CM 


<u 

(X 

o 

T - I 

CO 

'i 

•H 

P 

bO 

a 

•H 

t — I 
rH 
Cd 
Pm 


Steepest 

K(max) 

m 

oo 

ON 

in 

oo 

vO 

CO 

CO 

m 

CM 

i — 1 

o 

CM 

VO 

m 

vo 

tH 

CO 

LO 

00 

O 

m 

r- 

o 

CM 

00 

CO 

in 

O 

o 

o 

o 

O 

o 

co 

o 

O 

! - 1 

tH 

O 

o 

0) 

bO 

cd  > 
H  cd 
0)  v^- 

5  w 


CM 

CM 

CO 

<r 

CM 

ON 

CO 

00 

LO 

tH 

m 

co 

o 

r-. 

co 

o 

co 

vo 

CM 

CM 

CM 

vo 

m 

CM 

CO 

O' 

CO 

CM 

<U 

P- 

o 

rH 

CO 

Is 

•H 

i-3 

60 

C 

•H 

CO 

•H 

Pd 


4-1 

03 

o)  X 
CH  cd 
QJ  0 

4)  ' — ✓ 
4J  pd 
CO 


CM 

CO 

CM 

n- 

CO 

CO 

VO 

t — i 

CM 

n- 

ON 

vO 

vd- 

CO 

ON 

O 

o 

1 — 1 

oo 

VO 

CN 

CM 

1 — 1 

CM 

O 

o 

m 

LO 

m 

tH 

i — 1 

CM 

o 

o 

tH 

O 

o 

c- 

l — 1 

CM 

O 

o 

O 

o 

o 

CU 

bO  ^ 
cd  > 
H  cd 

0)  v_^ 

>  Fd 

<3 


n- 

vo 

in 

VO 

CM 

co 

ON 

vo 

m 

VO 

m 

CM 

00 

00 

1 - 1 

o 

in 

in 

ON 

CO 

CM 

ON 

1 — 1 

1 — 1 

CM 

o 

o 

r-- 

CO 

1 — 1 

tH 

tH 

o 

o 

a 

•  o 

o  *H 

•  4-1 

co  cd 

•  -U 

ts.  co 


o 

53 


1 — 1 

CO 

CM 

rH 

tH 

t — l 

m 

t — l 

tH 

1 — 1 

CM 

CO 

i — l 

o 

o 

O 

O 

O 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

O 

o 

o 

o 

Pd 

FQ 

o 

PH 

*o 

Pd 

Pd 

P 

w 

pH 

r-rt 
t— t 

<3 

o 

FQ 

FQ 

FQ 

P Q 

FQ 

FQ 

FQ 

u 

o 

O 

o 

Pd 

Pd 

r>. 

n- 

r"- 

r- 

o 

53 


r-~ 

CM 


oo 

CM 


CTn 

CM 


O 

CO 


CO 


CM 

CO 


CO 

CO 


CO 


m 

co 


vo 

co 


i"' 

CO 


oo 

CO 


O'. 

CO 
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cn 
co 
d)  = 


U  t-H 

CO 

CM 

vf 

Ox 

T— I 

o 

CO 

o 

CO 

o 

X  ^ 

P 

rv 

o 

m 

CM 

vO 

VO 

rv 

in 

t-H 

P  O 

co 

Ox 

c 

o 

c 

<3 

c 

co 

CM 

m 

CM 

vO 

rv 

P-< 

6  O' 

r\ 

#» 

«N 

“V. 

ex 

A 

ex 

ex 

ex 

ex 

o  • 

MD 

1 - 1 

CO 

?3 

121 

CO 

Ox 

i — 1 

co 

Ox 

Ci  co  •• 

CM 

r— 1 

CO 

vO 

o 

vO 

rH 

CD 

44  C  41 

rH 

0) 

PL, 


CL 

P-l 


CD 

PQ 

o  o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

co 

•  • 

o 

o 

m 

m 

o 

o 

m 

o 

o 

o 

o 

cD 

Es 

m 

I — 1 

CM 

rH 

m 

m 

CM 

vr 

o 

o 

CM 

PQ 

o 

rH 

P*H 

ex 

rH 

ex 

rH 

ex 

rH 

CD 

CL 

o 

r— I 

co 

’i 

P 

60 

a 

•H 

t— I 

i — ! 

CC5 

Cl, 


4J 

CO  / — - 

CM 

m 

a)  X 

Cl,  qj 

CM 

in 

CO 

co 

ox 

cd  e 

CD  v^ 

iH 

co 

m 

m 

O 

m 

CO 

rv 

co 

ax 

4-»  Ld 

CO 

rv 

ox 


vr 

LO 

vO 


CM 

O 


* 

OX 

vO 

vl- 

Vf 

o 

O 

CO 

Ox 

m 

CM 

rv 

m 

0) 

60  /*  N 

CD  > 

5l  c D 

a) 

>  w 

c 


CO 

CO 

vf 

CO 

vr 

Ox 

vO 

VI¬ 

* 

Vi¬ 

co 

rv 

CO 

o 

Vf 

tv 

CO 

<T 

VO 

VO 

co 

Ox 

co 

rv 

ox 

rv 

Ox 

OX 

OX 

rv 

UO 

vO 

m 

rv. 

VO 

CD 

CL 

O 

T— I 

m 

H 

•H 

p 

60 

a 

•H 

CO 

•H 

P4 


•U 

CO  k-s 
CD  * 
CL  CD 
CD  £ 
0) 

4J  PD 

CO 


Ox 

l — 1 

CM 

co 

CM 

co 

CO 

04 

VO 

CO 

OX 

CO 

CO 

CM 

CM 

CO 

m 

rH 

m 

m 

i — i 

rv 

m 

rv 

in 

l — 1 

O 

m 

co 

vo 

o 

o 

o 

co 

CM 

a) 

60 

CD  > 
Cl  CD 
CD  v^ 

3* 


04 

oo 

CM 

o 

04 

m 

ox 

00 

Ml 

CM 

O 

CM 

ox 

P 

Ox 

vi¬ 

OX 

in 

m 

P 

<r 

m 

vf 

i — 1 

rv 

vO 

m 

vo 

rv 

vf 

CM 

CM 

uo 

m 

m 

CM 

co 

rv 

00 

i — 1 

<r 

, — 1 

1 — 1 

1 — 1 

CM 

CM 

CM 

CO 

a 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 

O 

O 

O 

• 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 

O 

o 

O 

o 

•H 

• 

P 

< 

<3 

pq 

o 

u 

w 

o 

<3 

u 

O 

Q 

<3 

• 

4-1 

o 

<3 

<3 

<3 

<3 

<3 

< 

c 

« 

« 

n 

n 

w 

co 

CD 

13 

rv 

oo 

Ox 

Ox 

ox 

OX 

OX 

ox 

OX 

Ox 

OX 

OX 

Ox 

•  4J 

rs  co 


o 

J3 


O 

vf 


CM 

Vf 


CO 

VT 


Vt 

Vf 


IP 

Vf 


VO 

vr 


tv 

vi- 


co 

vf 


Ox 

Vf 


o 

m 


m 


CM 

m 
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CO 

CO 

0) 


£ 

O  rH 

00 

CM 

OV 

00 

rH 

\D 

CN 

i — 1 

o 

x 

r-~ 

CNJ 

o 

^el- 

00 

LO 

\ — 1 

oo 

iH 

w  o 

<1- 

CM 

rH 

oo 

in 

m 

00 

i — 1 

» — i 

c^ 

P4 

e 

O' 

r» 

•> 

#\ 

#\ 

r\ 

#\ 

r\ 

o 

• 

m 

00 

00 

r-* 

CM 

CO 

c-- 

m 

G 

CO 

oo 

CH 

m 

CT\ 

-3" 

CO 

oo 

rH 

oo 

ctf  4-1 

<u 
P-. 


g  4-* 

•H  Pu 
PM 


pH 


<D 

PQ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CO 

•  • 

o 

o 

o 

o 

o 

m 

o 

o 

1 — 1 

a) 

£ 

o 

vO 

CM 

o 

m 

i — i 

in 

o 

PQ 

o 

rH 

pH 

rv 

CM 

#1 

rH 

#v 

rH 

#» 

CM 

<u 

P. 

o 

r— I 
CO 

'i 

•H 

»-3 

00 

a 

•H 
rH 
t — I 

CCS 

fH 


CO  ✓ — s 

0)  x 

Pu  g 
a)  g 

<U  w 
4-J  Si 
CO 


t-- 

CO 

vO 

co 

m 

rH 

1 — 1 

r^. 

VO 

m 

CM 

CO 

o 

o 

o 

m 

m 

m 

in 

m 

vO 

CO 

o 

<D 

OO 

aj  > 
G  cd 
a)  s^ 

>  w 

•3 


o 

r^» 


-d- 

cr> 


f'' 

co 

i-- 

CM 

co 

CO 

o 

o 

CN 

<3* 

co 

CO 

r-'- 

vO 

vO 

vO 

vo 

co 


-d" 

oo 


QJ 

(X 

o 

rH 

CO 

•H 

rQ 

to 

G 

•H 

CO 

•H 

pi 


•u 

CO  ,-s 

as  x 
(X  as 

<u  e 

as  w 
x>  pi 

co 


vO 

CM 

CN 

VO 

CO 

CO 

o 

in 

vO 

CM 

-ci- 

in 

~4- 

vO 

CO 

O 

OV 

oo 

CM 

-d" 

in 

in 

vO 

i — 1 

o 

o 

CM 

CO 

O 

o 

o 

CM 

CN 

CM 

o 

o 

as 

00  ,-N 

CCS  > 
G  cd 
a)  v- ✓ 
>  Pi 

<3 


rH 

in 

CM 

CM 

VO 

r- 

CM 

CO 

r^- 

CM 

o 

o 

1 - 1 

rH 

in 

o 

CO 

CO 

in 

r- 

CN 

o 

• 

rH 

1 — 1 

CO 

rH 

1 — 1 

tH 

rH 

i — 1 

in 

rH 

G 

o 

o 

o 

O 

o 

O 

O 

O 

o 

CO 

• 

o 

o 

O 

o 

O 

o 

O 

O 

o 

o 

O 

o 

•H 

• 

<j 

PQ 

<3 

PQ 

u 

PQ 

PQ 

<1 

<d 

< 

• 

u 

o 

< 

<3 

w 

W 

w 

O 

W 

< 

< 

CO 

a) 

o 

o 

o 

O 

o 

O 

O 

o 

rH 

i — 1 

• 

-t-i 

rH 

1 — 1 

rH 

>H 

rH 

rH 

rH 

rH 

t — I 

rH 

3s  co 


o 

53 


oo 

40 


-d- 

m 


m 

m 


vO 

m 


r^- 

m 


oo 

m 


ov 

m 


o 

vO 


VO 


CM 

VO 
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CO 

CO 

0) 

* 

£ 

o 

1 — 1 

o 

X 

rH 

w 

o 

Pm 

£ 

ON 

o 

rH 

CO 

On 

1-"- 

NO 

o 

O 

• 

<t 

ON 

i — 1 

On 

CO 

CO 

CM 

CO 

5m 

CO 

•  • 

m 

i — i 

CM 

CO 

rH 

NO 

CM 

CO 

UM 

c 

4-1 

(U 

•H 

pH 

PH 

Pm 

PQ 

<15 

CO  •• 

CCJ  £ 

PQ  o 


0) 

ex 

o 

i — i 

'i 

*rl 

rJ 

tO 

c 

•H 

rH 

rH 

dj 

Pm 


4-) 

CO 

V  X 

PX  oJ 
0)  £ 
<U  w 

to 


0) 

to  ^ 
o5  > 
5m  CCJ 
cu  w 

>,  « 

< 


<1* 

NO 

m 

NO 

m 

1 — 1 

CM 

rH 

CO 

o 

00 

NO 

o 

O 

o 

o 

o 

o 

o 

o 

CD 

PX 

o 

f— I 
CO 

£ 

•H 

i-3 

to 

c 

*H 

CO 

•H 

PSh 


4J 

CO  r-v 
QJ  * 
PX  CO 

<U  £ 

a) 

+j 

CO 


a) 

tO  f— s 
CO  > 
5m  cO 
a)  v_^ 
>  Pd 
<3 


CO 

CO 

CM 

I — 1 

rH 

CO 

oo 

00 

CM 

LO 

O 

o 

CO 

rH 

CM 

o 

NO 

• 

• 

• 

• 

• 

• 

• 

• 

o 

o 

o 

o 

o 

O 

o 

o 

n 

•  o 

O  iH  • 

•  X>  O 
CO  CO  ^ 

•  X) 

JS  CO 


o 

• 

5m 

• 

• 

£ 

• 

U 

• 

Pm 

>M 

o 

5m 

5m 

• 

C_> 

o 

EH 

U 

5m 

CJ 

• 

5m 

(U 

Pm 

O 

C 

Pn 

Pm 

0) 

CJ 

Xi 

O 

o 

Q) 

5m 

C 

c 

5m 

CO 

O 

PX 

O 

d) 

u 

£ 

CO 

XI 

p 

PX 

O 

PX 

p 

o 

a) 

05 

y-1 

PQ 

S 

CO 

PQ 

EH 

PQ 

CO 

m 

NO 

r^. 

oo 

ON 

o 

VO 

NO 

NO 

NO 

NO 

NO 

NO 

o 

P5 


I 


TABLE  I 1-2  HYDROGRAPH  PARAMETERS  FOR  ADDITIONAL  SMALL  BASINS 


APPENDIX  III 


DRAINAGE  BASIN  CHARACTERISTICS 


101 


p 

p 

p 

o 

i — i 

c 

♦ 

>■4 

o 

O 

o 

o 

o 

m 

p 

LO 

o 

cO 

Ml 

CO 

n 

CO 

CO 

o 

o 

CN 

o 

00 

r-\ 

O 

p 

Cm 

p 

s 

C 

a 

• 

CO 

• 

• 

O 

in 

p 

1 — 1 

t — 1 

i — 1 

• 

• 

• 

<1 

TO 

o 

e» 

o 

o 

CO 

i — 1 

#v 

#v 

00 

p 

rH 

in 

« — 1 

g 

CO 

t — 1 

00 

CN 

o 

P 

Kf 

CN 

o 

p 

p 

CH 

oo 

a 

4-4 

o 

o 

• 

O 

o 

o 

4-1 

Mi 

TO 

p 

CO 

O'- 

p 

o 

O 

o 

o 

in 

o 

in 

CO 

pH 

p 

}-l 

c 

a 

• 

O 

in 

• 

• 

• 

o 

CN 

o 

CO 

p 

• 

• 

• 

<3 

as 

cO 

CO 

o 

p 

cO 

cO 

o 

o 

in 

co 

o 

cO 

P 

p 

i — 1 

in 

4-1 

4-1 

n 

r\ 

#» 

r» 

CO 

cO 

CO 

p 

CO 

p 

in 

& 

« 

p 

Jo 

15 

cO 

CO 

o 

OO 

o 

.  pq 

c 

O 

CO 

o 

CO 

cO 

o 

00 

o 

o 

o 

•o 

in 

CO 

Pm 

o 

as 

C2s 

e 

o 

1 — 1 

in 

cO 

cO 

p 

O' 

o 

in 

o 

in 

rH 

1 - 1 

p 

in 

<1 

i — 1 

• 

• 

• 

in 

p 

in 

in 

CO 

• 

• 

• 

LO 

4-1 

cO 

1 — 1 

o 

o 

e> 

rv 

«s 

CO 

o 

i — 1 

4-1 

u 

1 

CN 

i — 1 

CO 

rH 

m 

p 

• 

<N 

CO 

o 

p 

CN 

o 

£ 

a 

o 

00 

p 

I"'' 

o 

o 

o 

o 

in 

CO 

p 

CO 

eg 

o 

• 

cO 

o 

o 

in 

o 

• 

• 

• 

o 

o 

CO 

• 

• 

• 

Ph 

< 

1 — 1 

Mi 

1 - 1 

P 

1 — 1 

o 

o 

00 

i — 1 

p 

p 

co 

CO 

p 

rH 

m 

t — 1 

C 

O 

1 

r« 

*\ 

r» 

r» 

•H 

p 

p 

co 

p 

co 

rs 

•• 

p 

p 

• 

CJ 

p 

CO 

CO 

CN 

c 

• 

- 

« — i 

1 - i 

o 

cO 

Mi 

a 

tH 

p 

o 

o 

o 

o 

o 

in 

in 

00 

rH 

< 

e 

a 

o 

in 

o 

in 

o 

• 

o 

• 

o 

o 

o 

o 

CN 

• 

• 

• 

P^ 

< 

'b 

Ml 

in 

P 

CN 

p 

in 

p 

1 — 1 

o 

o 

1 — 1 

p 

i — 1 

in 

r~ 1 

TO 

1 

r> 

ri 

rc 

n 

o 

T — 1 

00 

o 

CO 

in 

cO 

rH 

P 

• 

o 

/"N 

X 

to 

ta 

• 

• 

• 

• 

• 

• 

Sz; 

o 

CJ  cJ 

o 

<3 

• 

• 

p 

pH 

• 

co 

H 

p 

w 

• 

•  o 

a) 

P 

<3 

<3 

cj 

CO 

<3 

<3 

w 

w 

w 

pL] 

w  - 

CJ 

•  <3 

o 

o 

CO  *H 

+J 

H 

P 

pH 

g 

p 

pq 

M 

p 

o 

g 

S 

PQ 

g 

P 

cm 

a 

IS 

•  4-1 

aj 

Cm 

cd 

4-1 

tz; 

4-J 

CO 

CO 

1  Ml 

as 

• 

Ml  4-1 

o 

1 — 1 

CN 

CO 

<r 

in 

CO 

P 

CO 

p 

o 

rH 

CN 

CO 

in 

CO  Cl) 
CM  6 

p 

p 

t — 1 

P 

1 — 1 

1 — 1 

1 — 1 

Foothills;  P  =  Prairies;  A  =  Aspen  Parkland;  B  =  Boreal  Forest;  YT  =  Yukon  or  N.W.T. 

table  iii-i  drainage  basin  characteristics 


102 


3 - 1 

T3 

<S> 

O 

a) 

o 

O 

P3 

-d 

CN 

1 — 1 

pq 

• 

4-1 

m 

O 

CN 

oo 

r-. 

o 

o 

o 

m 

i — 1 

O 

O 

a) 

P3 

d 

CN 

• 

• 

• 

• 

o 

o 

o 

O 

• 

• 

<3 

LO 

rH 

o 

Os 

o 

1 — 1 

i — 1 

o 

o 

m 

o 

m 

m 

vO 

00 

CN 

1 - 1 

4-1 

3-4 

s 

CO 

m 

4J 

03 

m 

CO 

m 

CN 

CO 

•H 

03 

Q 

pi 

1 - 1 

P 

03 

o 

o 

o 

o 

in 

o 

o 

03 

• 

P 

o 

CN 

CN 

o 

o 

o 

o 

o 

r- 

CO 

1"- 

CN 

as 

o 

a) 

Pi 

"d 

m 

o 

CO 

• 

• 

o 

o 

crv 

o 

1 — 1 

o 

• 

• 

• 

pH 

<3 

Q 

fi 

d 

as 

m 

rH 

r\ 

r\ 

#» 

#v 

1 — 1 

t - 1 

o 

d 

o 

CO 

o 

vO 

m 

CO 

in 

in 

x) 

CO 

0) 

Pd 

03 

• 

'd 

o 

P5 

5o 

<N 

• 

co- 

o 

o 

o 

O 

o 

i — 1 

oo 

o 

a, 

P 

P 

o 

O 

o 

o 

o 

CN 

o 

CO 

o 

o 

n- 

i-3 

(13 

d 

03 

VO 

o 

vO 

co 

• 

• 

o 

vO 

CO 

VO 

CN 

m 

• 

• 

• 

pH 

pq 

a) 

*d 

CD 

m 

1 — 1 

*s 

#v 

e* 

#v 

m 

o 

i — i 

m 

-d 

rH 

O 

CO 

o 

P- 

in 

oo 

c 

1 - 1 

rH 

03 

"d 

T3 

as 

o 

5o 

o 

o 

o 

• 

(0 

CD 

in 

n- 

r- 

o 

o 

o 

O 

in 

o 

<r 

r- 

o 

£ 

Pt 

P 

O 

o 

VO 

• 

• 

• 

o 

O 

o 

o 

CN 

• 

• 

• 

PH 

hJ 

-d 

a 

0) 

Os 

O 

o 

o 

■H 

VO 

VO 

o 

rH 

pq 

P0 

"d 

#v 

r. 

A 

r\ 

e» 

m 

m 

•H 

rH 

CO 

CO 

O 

VO 

w 

<3 

• 

'• 

1 - 1 

• 

Pi 

4-1 

a 

CCS 

in 

K 

o 

CO 

vO 

o 

a 

03 

CN 

ffl 

o 

ca> 

d 

<T 

CN 

O 

cn 

n- 

o 

o 

o 

o 

o 

(N 

o 

m 

< 

(0 

•H 

• 

o 

OV 

• 

• 

• 

o 

o 

o 

o 

m 

• 

• 

p-t 

m 

3-i 

C3 

o . 

o 

o 

o 

00 

CN 

CO 

VO 

CO 

CO 

o 

1 — 1 

ai 

•H 

1 — 1 

C\ 

A 

#N 

1 - 1 

m 

pq 

T3 

CN 

CN 

CN 

03 

5S 

• 

o 

'x 

53 

• 

d 

<3 

• 

• 

hJ 

Pq 

CO 

H 

53 

H 

O 

M 

-S 

o 

O 

< 

<3 

u 

CO 

<3 

<3 

w 

pq 

pq 

pq 

W 

o 

<£ 

<3 

O 

O 

•  o 

CO  *H 

*H 

4-J 

Q) 

H 

pH 

a 

pq 

pq 

pq 

Q 

o 

s 

2 

pq 

£ 

2 

g 

Is 

•  -u 

ctf 

§ 

13  03 

4-1 

53 

•M 

CO 

C/3 

1  p 

03  03 

• 

3-i  4J 

o 

1 — 1 

CN 

oo 

m 

VO 

r- 

CO 

ov 

o 

i — 1 

CN 

CO 

m 

d  a) 

5Z 

rH 

rH 

tH 

1 — 1 

rH 

rH 

Ph  b 

TABLE  1 1 1 - 1  (CQNT)  DRAINAGE  BASIN  CHARACTERISTICS 


103 


• 

0  w 

CM 

r— ^  H 

o 

CO  (1) 

o 

o 

o 

o 

o 

o 

m 

m 

o 

cO  chj  to 

00 

CM 

m 

CO 

o 

o 

o 

o 

o 

o 

m 

cO 

CM 

P 

o 

CO  0 

C0 

• 

00 

• 

• 

o 

o 

<r 

• 

• 

• 

P^H 

Q 

3 

O 

i — 1 

#v 

rv 

#\ 

CM 

rH 

LO 

•  cO 

rH 

C0 

CO 

1 — 1 

m 

S  to 

rH 

•  • 

C*3  0) 

CO 

Pi  K 

rH 

CL) 

O 

o 

o 

o 

o 

o 

<r 

o 

P  •  0 

1 - 1 

m 

m 

1 — 1 

r-~ 

o 

m 

o 

m 

in 

o 

CO 

o 

0  Pi  P 

CM 

o 

C0 

o 

• 

• 

• 

m 

i — 1 

o 

00 

m 

• 

• 

• 

fH 

\ — 1 

PQ 

g  0  £3 

rH 

C0 

o 

o 

e\ 

r\ 

r\ 

CM 

CM 

i — i 

Q 

1 - 1 

oc 

CO 

rH 

CO 

m 

LO 

cO  to 

d)  P 

t — 1 

i — 1  o 

u  o 

Old 

• 

Pd  TO 

LO 

0 

00 

o 

TO  CO 

O 

o 

o 

o 

o 

o 

CO 

CM 

cO 

CO 

o 

0  i — 1 

'o 

r-~- 

o 

CO 

0- 

o 

o 

o 

o 

CO 

• 

•  * 

* 

Pm 

iH 

w 

0  ©  0 

<!■ 

• 

m 

CM 

• 

• 

• 

cO 

cO 

cO 

o 

o 

1 - 1 

rH 

o 

0)  0) 

#v 

o 

1 - 1 

o 

o 

o 

rv 

€S 

e\ 

r\ 

m 

uo 

CO  Pd 

o 

Pd 

1 — 1 

1 — 1 

1 

CO 

CM 

CO 

rH 

CO 

Cr. 

r 

0) 

CM 

CO  i — 1 

o 

1 — 1  •  t — 1 

o 

00 

o 

rH  P  CD 

m 

o 

00 

CsJ 

m 

r-» 

m 

O 

m 

m 

o 

cO 

CM 

w 

•H  0  P 

ac 

• 

• 

• 

00 

CM 

CO 

♦ 

• 

• 

PM 

1 - 1 

CJ 

.0  0 

m 

o 

o 

CM 

1 — 1 

CO. 

i — 1 

cO 

00 

i — 1 

\ - 1 

m 

0)  0 

i 

#N 

rs 

r\ 

in 

d)  P 

CO 

CM 

CO 

i — 1 

nj 

0  Q 

• 

Pi 

CL) 

d)  rH 

CM 

rH 

o 

C  rH 

<T 

m 

in 

t - 1 

f-. 

n- 

o 

O 

o 

o 

O 

cO 

CO 

m 

i — 1 

o 

•H  •  *H 

in 

• 

in 

m 

• 

• 

• 

m 

o 

CO 

CD 

co 

• 

• 

• 

<3 

o 

U  Pi  > 

r^. 

o 

CO 

o 

o 

CO 

o 

-<r 

CM 

CM 

CM 

rH 

o 

•H  0  P 

T - 1 

i. 

n 

«s 

r» 

m 

m 

to  a 

CO 

CO- 

CO 

CO 

d)  w 

S 

• 

• 

• 

o 

X 

Fo 

53 

• 

• 

• 

• 

P 

O 

• 

0 

<3 

• 

• 

P 

Fh 

CO 

H 

H 

H 

o  0 

O 

PI 

<3. 

<3 

u 

CO 

<3  . 

<3 

w 

w 

w 

W 

o 

<3 

<c 

O 

• 

•  O 

•H  d) 

H 

P 

p 

S 

pp 

pp 

pp  ■ 

Q 

o 

§ 

2 

PP 

2 

P 

2 

W 

o 

CO  *H 

p  6 

2 

•  p 

aj  c3 

32  03 

4-J 

p 

CO 

CO 

1  u 

CO  CL) 

• 

Pi  p 

o 

rH 

CM 

CO 

in 

VO 

n- 

00 

00 

o 

i — 1 

CM 

CO 

m 

cO  a) 
PM  0 

1 — 1 

rH 

P 

rH 

rH 

p 

TABLE  III-I  (CONT)  DRAINAGE  BASIN  CHARACTERISTICS 


104 


20 

5FE001 

Battle  R. 

nr . 

Unwin 

O 

o 

1 — 1 

o 

1 — 1 

m 

20 

354 

12.5 

o 

• 

1 — 1 

O 

3,010 

1,750 

3,125 

1,375 

i 

2,195 

2.0 

52.95 

r-' 

• 

rH 

<3 

• 

P 

3 

u 

o 

i — 1 

p 

in 

o 

03 

• 

p 

o 

CO 

00 

O' 

r-' 

o 

m 

o 

m 

in 

CO 

OC 

o 

3 

p 

a) 

OC 

t - 1 

1 — 1 

00 

• 

• 

• 

■H 

CO 

CN 

m 

• 

• 

• 

c 

rH 

Q 

o 

3 

to 

cn 

p". 

o 

o 

CO 

o 

CO 

m 

CN 

CO 

CN 

Ph 

p 

T3 

C\ 

rv 

r» 

r\ 

m 

UO 

CO 

w 

CN 

CN 

CN 

CN 

•H 

P<3 

rH 

3 

o 

03 

• 

P3 

cO 

oo 

■  o 

rH 

P 

O 

rH 

o 

CN 

CO 

r-' 

r'- 

o 

O 

O 

o 

m 

CO 

co 

oo 

1 — 1 

<3 

P 

3 

3 

t — 1 

• 

m 

m 

• 

• 

• 

o 

in 

O 

m 

o 

• 

• 

• 

c 

fH 

P 

O 

r- 

o 

CO 

o 

o 

CN 

cO 

CN 

in 

CN 

CN 

i — 1 

m 

CO 

P-: 

l 

r\ 

ev 

r« 

in 

FP 

CO 

CN 

CO 

CN 

pi 

03 

1 - 1 

00 

C3 

1 - 1 

o 

00 

o 

o 

m 

o 

m 

o 

cO 

co 

1 — 1 

O 

• 

•H 

o 

• 

cn 

Ol 

• 

• 

m 

o 

in 

cr> 

• 

• 

• 

<3 

r-'- 

o 

*H 

P 

> 

CN 

<3* 

i — 1 

Ol 

cO 

rH 

o 

CO 

oc 

CO 

O 

CN 

CO 

1 — 1 

1 — 1 

o 

i — 1 

3 

3 

r\ 

1 — 1 

in 

w 

•H 

3 

CN 

CN 

I — I 

CN 

CN 

w 

0 

3 

. 

in 

p 

S 

03 

’  >' 

• 

•  * 

pi 

• 

3 

ac 

rH 

a 

PS 

o 

on 

1 - 1 

o 

O 

o 

o 

o 

O 

r-. 

co 

cO 

o 

o 

• 

3 

1 — 1 

• 

LO 

cO 

• 

• 

o 

o 

o 

o 

rH 

• 

• 

• 

<3 

rH 

o 

03 

P 

3 

CO 

CN 

CN 

CN 

t — 1 

o 

CO 

o 

r^- 

CN 

CO 

co 

l — 1 

<1 

00 

3 

CO 

•N 

i-H 

rH 

rc 

C\ 

r\ 

m 

w 

P 

rH- 

CN  • 

CN 

CN 

CN 

m 

3 

• 

P 

p 

C/5 

[H 

• 

✓“s 

o 

X 

13 

• 

• 

• 

V — 

• 

• 

S 

O 

3 

<3 

• 

• 

pi 

m 

H 

H 

H 

o  3 

O 

Q 

C 

<3 

o 

CO 

<J 

w 

w 

W 

W 

W 

V 

C 

O 

o 

•  o 

CO  -H 

•H 

4_J 

3 

H 

PI 

U-i 

S 

FP 

FP 

CP 

Q 

o 

§ 

s 

FP 

S3 

PI 

p2 

S3 

•  P 

cd 

IS  03 

p 

2 

P 

in 

m 

1  P 

03  <13 

• 

P  P 

o 

>H 

04 

co 

m 

cO 

r- 

00 

OC 

o 

rH 

CN 

CO 

NT 

in 

03  0) 

iz; 

rH 

rH 

rH 

rH 

1 — l 

rH 

FP  0 

TABLE  III- I  (CONT)  DRAINAGE  BASIN  CHARACTERISTICS 


105 


Pi 

pi 

CO 

CM 

0) 

o! 

O 

o 

o 

o 

o 

m 

O' 

O 

T3 

> 

S-i 

1 - 1 

m 

O 

CM 

O' 

o 

o 

CO 

CO 

O' 

o 

to 

O 

O 

o 

34 

O 

• 

at 

at 

• 

• 

• 

m 

x — 1 

m 

at 

• 

• 

• 

IH 

CM 

f-H 

03 

.40 

03 

#\ 

o 

00 

o 

o 

r\ 

rt 

n 

r\ 

rH 

CO 

X — 1 

<1 

.-3 

03 

-O 

1 - 1 

to 

CO 

00 

in 

CO 

m 

O' 

O 

E 

g 

W 

Pi 

03 

1 — 1 

cti 

o 

m 

o 

o 

o 

o 

o 

00 

o 

O 

a 

? — 1 

to 

o 

CO 

o 

r-' 

o- 

o 

m 

o 

m 

at 

x — 1 

co 

CM 

o 

CO 

CH4 

03 

at 

• 

to 

o 

• 

• 

• 

CO 

CM 

CO 

o 

00 

• 

• 

• 

Pm 

•  CM 

0 

cO 

>-i 

o 

1 — 1 

CO 

o 

o 

rv‘ 

r\ 

#\ 

#N 

in 

CO 

X - 1 

<3 

rQ 

4J 

CO 

1 

CM 

CO 

CM 

at 

CO 

x — I 

in 

oj 

C 

rH 

X — 1 

rC 

w 

4-1 

c 

c 

03 

•H 

T3 

X — 1 

C 

a 

o 

M 

• 

o 

O' 

CO 

o 

!-i 

CO 

1 — 1 

CO 

o 

o 

at 

O' 

i"-. 

O 

o 

o 

o 

o 

X — 1 

1 — 1 

CM 

0 

03 

C 

03 

l — 1 

• 

to 

to 

• 

• 

• 

O 

o 

m 

• 

• 

• 

< 

Pi 

03 

at 

CO 

o 

o 

o 

m 

to 

X - 1 

O' 

CO 

1 — 1 

03 

PQ 

1 

#\ 

<T 

in 

fl 

CO 

CO 

at 

to 

CO 

• 

CO 

03 

• 

Pi 

to 

Pi 

o 

03 

o 

CM 

o 

3-4 

AJ 

to 

o 

o 

o 

o 

o 

to 

CM 

Q 

03 

<S> 

03 

m 

o 

m 

O' 

o 

in 

o 

m 

CM 

CM 

CO 

<T 

<3 

> 

h-1 

r» 

CO 

in 

• 

• 

• 

CO 

to 

m 

00 

at 

• 

• 

• 

PQ 

VO 

03 

m 

r— 1 

o 

o 

r» 

4\ 

r\ 

r\ 

CO 

X — 1 

0) 

03 

1 

CM 

rH 

CM 

rH 

m 

PQ 

1 - 1 

O 

u 

• 

3m 

CJ 

34 

CO 

O 

r-' 

o 

3-4 

• 

4-> 

m 

l — 1 

CO 

r-' 

O' 

o 

o 

O 

o 

in 

to 

at 

O' 

X - 1 

o 

O 

34 

•H 

to 

• 

o 

rH 

• 

♦ 

• 

o 

o 

o 

o 

at 

• 

• 

• 

P4 

CM 

<3 

4-1 

ti 

a 

m 

CM 

<r 

O' 

o 

o 

to 

CM 

to 

CO 

O' 

rH 

rH 

o 

•H 

o 

r» 

r\ 

m 

to 

c 

g 

CM 

CM 

CM 

CM 

o 

g 

• 

o 

X 

3o 

£5 

Jz; 

• 

• 

• 

• 

• 

33 

O 

a 

<3 

• 

• 

IH 

C/3 

H 

H 

w 

o  c 

o 

a 

<3 

<3 

o 

CO 

<3 

<3 

W 

w 

w 

w 

O 

<3 

<3 

o 

O 

•  o 

in  *h 

•H 

4-> 

03 

H 

r-Q 

Pm 

g 

PQ 

PQ 

PQ 

« 

o 

s 

3 

PQ 

g 

g 

32 

•  -h 

cd 

Ls  o3 

4-> 

33 

4-> 

CO 

CO 

1  U 

03  0) 

• 

34  4-1 

o 

i— 1 

CM 

CO 

m 

to 

o- 

CO 

at 

o 

rH 

CM 

CO 

m 

03  0) 

32 

rH 

rH 

rH 

rH 

rH 

x — 1 

PM  0 

TABLE  III-I  ( CONT)  DRAINAGE  BASIN  CHARACTERISTICS 


106 


Pi 

rH 

a) 

• 

03 

o 

o 

•H 

Ph 

Td 

o 

LO 

O' 

OH 

CO 

o 

pH 

d 

rH 

o 

CO 

CO 

LO 

• 

• 

• 

LO 

o 

LO 

LO 

o 

1 - 1 

Pm 

•H 

•H 

LO 

cd 

oo 

CD 

o 

o 

o 

<3* 

O' 

CO 

O' 

• 

• 

• 

pq 

pq 

0) 

d 

tH 

1 

CD 

<1* 

LO 

CO 

CD 

LO 

t - 1 

r^. 

PH 

w 

r\ 

#» 

#\ 

LO 

PH 

1 - 1 

OH 

CM 

w 

Pi 

O) 

03 

<13 

o 

o 

O 

o 

O 

LO 

LO 

o 

d 

•H 

m 

LO 

o 

OH 

O' 

o 

o 

o 

LO 

LO 

CM 

LO 

LO 

CD 

o 

•H 

CS) 

> 

LO 

• 

LO 

o 

• 

• 

• 

LO 

o 

00 

00 

o 

• 

• 

• 

pq 

o^ 

u 

rO 

Ph 

r» 

rH 

CD 

LO 

o 

o 

•s 

r\ 

e» 

c\ 

0\ 

cO 

<T 

1 - 1 

CQ 

s 

03 

<r 

co 

O' 

OH 

O' 

LO 

CO 

LO 

<V 

PH 

Pi 

CO 

rX 

i — 1 

.  o 

CJ 

• 

03 

• 

o 

•H 

Ph 

Po 

o 

LO 

LO 

LO 

•  o 

o 

o 

o 

o 

LO 

LO 

rH 

CO 

00 

pq 

4J 

d 

4-) 

O' 

• 

• 

cO 

• 

• 

o 

o 

o 

o 

o 

• 

• 

• 

pq 

CM 

pq 

C/3 

CO 

cO 

cO 

cO 

LO 

1 — 1 

LO 

CD 

O' 

oo 

CO 

1 — 1 

o- 

03 

rv 

rs 

r\ 

LO 

O 

CO 

CM 

CO 

OH 

)d 

Pi 

• 

Ph 

CO 

1 — 1 

£ 

O 

O 

o 

o 

o 

o 

o 

CD 

1 — 1 

LO 

o 

0) 

o 

rH 

LO 

o 

O') 

o 

o 

cO 

LO 

CD 

o 

• 

• 

• 

pq 

o 

d 

rH 

r — 1 

• 

00 

00 

• 

• 

o 

LO 

00 

O 

O' 

1 - 1 

CO 

I - 1 

O') 

<3 

•H 

a) 

03 

r\ 

o 

OC 

t — 1 

CD 

rH 

*\ 

r. 

r\ 

<D 

1 - 1 

LO 

pq 

.  03 

PQ 

od 

3 - 1 

CM 

O' 

OH 

O'  . 

LO 

CO 

i-' 

0 

03 

CL) 

Pm 

Ph 

• 

•  • 

0) 

d 

•H 

• 

o 

Pi 

03 

rH 

•H 

cO 

o 

c 

d 

OH 

CO 

o 

o 

o 

o 

o 

o 

LO 

CO 

LO 

cO 

o 

03 

• 

•H 

CD 

o 

CD 

o 

• 

• 

o 

o 

o 

o 

o 

o 

• 

• 

• 

pq 

Cs) 

£nj 

0 

Ph 

03 

CO 

CD 

rH 

cO 

1 — 1 

oo 

o 

00 

00 

CO 

1 — 1 

<fr 

1 - 1 

<tj 

<1) 

d 

03 

t - 1 

r* 

r\ 

#\ 

rc 

rH 

LO 

r-' 

a> 

<3 

CO 

OH 

CO 

1 — 1 

CM 

pH 

Ph 

• 

4-1 

pH 

• 

. 

■ 

/-N 

o 

X 

Po 

a 

• 

• 

• 

'wh' 

• 

• 

a 

o 

d 

<3 

• 

• 

►d 

Pm 

• 

• 

• 

r\ 

• 

• 

• 

CO 

H 

H 

H 

o 

a 

<3 

<3 

CJ 

•C/3 

<3 

<3 

w 

w 

pq 

Pd 

w 

u 

< 

o 

o 

•  o 

CO  *H 

•rH 

4-J 

03 

H 

hd 

Pm 

S 

pp 

pq 

pq 

Q 

o 

a 

3 

pq 

a 

hd 

3 

g 

P3 

•  4-1 

aj 

|2  03 

4-» 

a 

4-) 

CO 

CO 

1  Ph 

o3  a) 

• 

PH  4-> 

o 

rH 

O) 

CO 

<fr 

LO 

CO 

o~ 

oo 

CD 

o 

rH 

CM 

CO 

LO 

Cd  (U 

pa 

1 — 1 

rH 

i — i 

rH 

rH 

1 - 1 

Pd  0 

TABLE  III-I  (CONT)  DRAINAGE  BASIN  CHARACTERISTICS 


107 


3 

•H 

• 

p 

pc; 

•H 

1 — 1 

cd 

o 

•H 

•  P 

LO 

o 

■p 

P  Pn 

O 

o 

o 

O 

o 

m 

rH 

CO 

w 

•H 

3 

m 

m 

00 

o 

o 

o 

O 

o 

r^. 

CO 

I"- 

CM 

o 

PC 

3 

CO 

o 

Ol 

CN 

• 

• 

O 

m 

VO 

ON 

n- 

• 

• 

• 

Pi 

r- 

c0 

33 

rv 

1 — 1 

1 — 1 

rv 

rv 

rv 

r\ 

rv 

o 

m 

rH 

12 

3 

rH 

in 

CM 

CM 

m 

3 

O 

pi 

p 

rH 

<u 

o 

p 

P 

o 

o 

o 

o 

o 

O 

r- 

o 

3 

3 

o 

r'- 

00 

CO 

o 

o 

o 

o 

o 

o 

CM 

in 

VO 

Ov 

CO 

o 

£ 

<B;  0- 

00 

• 

Ov 

oc 

• 

• 

o 

00 

CO 

OV 

1 — 1 

CM 

• 

• 

• 

PQ 

CJ 

p 

3 

rv 

,-P 

CM 

rH 

rv 

rv 

rv 

rv 

CO 

vO 

o 

r-» 

3 

P 

1 — 1 

i 

1 — 1 

rH 

t — 1 

1 — i 

m 

a> 

Q 

T— 1 

i — 1 

CJ 

pi 

3) 

m 

X 

5 - 1 

o 

2 

3 

vO 

CO 

o 

.  co 

•  *H 

o 

O 

CO 

r>. 

r- 

o 

o 

o 

o 

■  o 

1 - 1 

CM 

CO 

P 

<D 

P  m 

CO 

• 

00 

CM 

• 

• 

• 

o 

CM 

o 

CO 

Ov 

• 

• 

• 

PQ 

PQ 

P 

3  P 

o 

rH 

ctv 

t - 1 

1 — 1 

o 

o 

o 

CT\ 

o 

o 

CM 

vO 

m 

1 - 1 

rP 

3 

rv 

t - 1 

t — 1 

n 

rv 

rv 

rv 

rv 

m 

3 

PC 

rH 

1 — 1 

CM 

CM 

cd 

CO 

CO 

• 

Pi 

3 

1 — 1 

3 

Pi 

o 

> 

3 

o 

o 

o 

o 

o 

O 

o 

Csl 

o 

cd 

3) 

o 

o 

o 

r>. 

n- 

o 

in 

o 

m 

o 

CM 

CO 

CM 

CO 

p 

1 — I 

<3J 

CO 

CM 

00 

00 

vO 

• 

• 

CO 

CO 

CO 

o 

• 

• 

• 

PQ 

PQ 

CO 

3 

1 — 1 

o 

o 

rv 

rv 

rv 

rv 

•v 

CM 

m 

1 - 1 

I''- 

> 

m 

i 

CM 

1 — 1 

CM 

CM 

m 

• 

3 

rP 

CO  . 

• 

rH 

Pi 

• 

o 

P  o 

CO 

1 — l 

o 

3 

3  3 

CNJ 

in 

o 

o 

o 

o 

o 

m 

I''. 

CO 

CM 

CO 

•o 

cd 

3 

o 

00 

« 

o 

• 

o 

o 

VO 

oo 

• 

• 

• 

PQ 

PQ 

IS 

3 

CTv 

VO 

m 

tH 

Ov 

-Cl 

CM 

LO 

m 

I - 1 

n- 

CO 

•H 

CM 

LO  • 

1 - 1 

CM 

CO 

• 

o 

a 

X 

P 

a 

a 

• 

• 

• 

• 

• 

o 

u  3 

•  O 

C 

• 

• 

P 

P 

• 

CO 

H 

H 

M 

, 

o 

3 

Q 

<3 

<3 

o 

CO 

<3 

<3 

p 

w 

w 

s 

w 

CJ 

<3 

<5 

o 

o 

CO  *H 

4_J 

H 

p 

P 

a 

PQ 

PQ 

PQ 

Q 

o 

§ 

PQ 

a 

P 

pi 

S 

a 

•  4-J 

a3 

Ph 

13  3 

p 

!2i 

4-t 

CO 

CO 

1  5-i 

cd  a) 

• 

O 

rH 

CM 

CO 

m 

vO 

P' 

00 

Ov 

o 

rP 

CM 

co 

m 

3  a; 

PM  0 

33 

rH 

1 - 1 

rH 

rH 

1 — l 

1 — 1 

TABLE  III- I  (CONT)  DRAINAGE  BASIN  CHARACTERISTICS 


108 


c 

o 

•H 

• 

3 — 1 

CN 

Pd 

rH 

o 

• 

•H 

O 

o 

co 

m 

o 

u 

u 

£ 

CM 

LO 

o 

C'- 

00 

o 

in 

03 

a 

u 

<d“ 

• 

CO 

o 

3 - 1 

r^. 

uo 

o 

m 

CO 

o 

m 

r^- 

<3- 

l-) 

Po 

03 

*\ 

CO 

OV 

o 

• 

• 

• 

av 

»v 

co 

r\ 

• 

• 

• 

pq 

r^. 

o 

> 

CM 

rH 

o 

o 

1 — 1 

00 

3 - 1 

vO 

rH 

co 

CO 

o 

o 

pq 

• 

m 

4-1 

Pm 

P i 

C 

00 

i — I 

•H 

a 

O 

O 

o 

m 

m 

m 

CM 

o 

£ 

•H 

t - 1 

CM 

CSV 

00 

o 

CM 

Ov 

OV 

ov 

o 

0) 

<BJ 

a 

CO 

CM 

00 

CO 

CO 

o 

m 

m 

o 

3 - 1 

• 

• 

• 

pq 

CO 

o 

Pd 

C3 

r\ 

CT\ 

m 

1 — I 

• 

o 

vO 

vO 

o 

jn 

•H 

03 

3 — 1 

1 — 1 

3 - 1 

CO 

3 - 1 

CO 

CM 

CM 

m 

4-1 

s 

. 

O 

Pd 

Pd 

CO 

• 

03 

.  o 

4-1 

3-1 

CM 

00 

o 

3-1 

a 

£ 

LO 

CO 

o 

o 

O 

o 

o 

O 

CM 

3 - I 

o 

CO 

<5 

03 

03 

LO) 

• 

O 

t"- 

• 

• 

o 

o 

CM 

o 

00 

r-^ 

• 

• 

• 

pq 

03 

Pd 

3 — 1 

OV 

LO 

1 — 1 

CO 

m 

vO 

o 

vO 

3 - 1 

Pd 

CM 

1 

r\ 

#\ 

rv 

m 

CM 

rH 

CM 

CM 

Pd 

CM 

to 

• 

o 

O 

o 

O 

O 

o 

VD 

o 

03 

3-i 

Jo 

co 

LTl 

o 

o 

o 

o 

o 

m 

CM 

OV 

r-'. 

o 

Pd 

C 

d 

3 - 1 

• 

00 

o 

• 

o 

• 

o 

vD 

o 

r^~ 

• 

• 

• 

pq 

CO 

Pd 

O 

o 

r\ 

1 - 1 

<Tn 

vO 

VO 

rH 

r\ 

r\ 

rv 

r\ 

CT\ 

m 

o 

o 

£ 

C-J 

3 - 1 

in 

rH 

m 

CM 

CM 

m 

r->- 

CO 

• 

• 

f-J 

Pd 

a 

3 — 1 

03 

•H 

o 

o 

o 

o 

o 

o 

vO 

o 

4-1 

• 

> 

CM 

o 

o 

o 

vO 

m 

m 

rH 

CM 

CO 

o 

4-1 

3-i 

od 

av 

o 

Cv 

CM 

• 

• 

o 

vO 

o 

CO 

• 

• 

• 

pq 

Pm 

03 

a 

o 

OO 

140 

o 

3 - 1 

r\ 

«N 

#v 

rv 

m 

m 

3 - 1 

o 

C 

o 

3 — 1 

3 - 1 

3 - 1 

m 

rH 

m 

CM 

CM 

3 - 1 

m 

O 

o 

£ 

•H 

CO 

• 

/'“N 

X 

Jo 

Pd 

• 

• 

• 

• 

• 

jz; 

O 

£3 

<3 

• 

• 

hJ 

m 

• 

CO 

H 

H 

H 

V-* 

o 

Q 

<J 

<3 

O 

CO 

<3 

<1 

w 

w 

pq 

W 

W 

o 

<3 

O 

o 

CO  -H 

•H 

4_J 

Cl) 

H 

hJ 

Pm 

S 

pq 

pq 

pq 

Q 

o 

§ 

3 

pq 

s 

3 

§ 

Pd 

•  4-> 

Ctf 

Ph 

13  o3 

4-> 

23 

4-1 

CO 

CO 

1  u 

03  CD 

• 

3-1  4-> 

o 

rH 

CM 

CO 

in 

vO 

r-~ 

00 

CV 

o 

i — 1 

CM 

co 

<3- 

m 

03  a) 

jz; 

rH 

rH 

3 - i 

rH 

rH 

rH 

PM  6 

TABLE  III- I  (CONT)  DRAINAGE  BASIN  CHARACTERISTICS 


109 


3— 

CO 

• 

£ 

P3 

cO 

CO 

O 

O 

•H 

£ 

•H 

S3 

(*3 

<24 

O 

O 

0 

m 

m 

0 

1 — 1 

Cj 

00 

1 — 1 

O 

0 

n- 

0 

m 

CO 

m 

<J 

rC 

_vj 

m 

• 

m 

O 

• 

• 

• 

m 

CM 

CM 

O 

\ — 1 

• 

• 

• 

H 

O'! 

n- 

CM 

00 

0 

0 

r\ 

r» 

rt 

rc 

CM 

O 

0 

>4 

£ 

1 — 1 

m 

CM 

co 

CO 

co 

CM 

CO 

H 

44 

S3 

O 

• 

0) 

Pd 

CO 

rH 

O 

•H 

•  0 

i.  1  r' 

O 

Oc 

CM 

CO 

0 

0 

0 

m 

m 

m 

m 

<0 

m 

0 

<1 

cr> 

a  a) 

4-J 

•H 

cO 

• 

cO 

T — 1 

• 

• 

0 

m 

1^. 

CM 

CM 

• 

• 

• 

H 

& 

cO 

cc 

CM 

CO 

1 — l 

0 

rH 

m 

#v 

m 

1 — I 

*s 

CM 

CM 

co 

1 — 1 

rc 

P 

re 

CM 

CTc 

0 

co 

0 

>4 

H 

pi 

03 

CO 

P 

U 

00 

O 

O 

co 

O 

O 

•  42 

Jh  03 

t - 1 

p 

0 

0 

0 

0 

b 

0 

Oc 

v£> 

in 

cr> 

• 

rH 

O 

• 

• 

0 

m 

<0 

CO 

m 

• 

• 

• 

H 

9 

*H 

LO 

(XI 

CO 

CM 

O 

p 

m 

CO 

co 

m 

<r 

CM 

0 

0 

>4 

CO 

rH 

0 

1 

#s 

CM 

#s 

CM 

co 

r\ 

rc 

CM 

cO 

a 

pi 

ly> 

•  r4 

O 

m 

CT\ 

1 — 1 

O 

O 

ac 

Oc 

1 — 1 

CO 

CO 

0 

co 

CM 

<3- 

0 

oc 

m 

• 

in 

• 

• 

O 

O 

co 

co 

CO 

• 

• 

• 

EH 

O 

40 

40 

C  P 

44 

<3 

CO 

<r 

OC 

CM 

0 

0 

0 

1 

<3* 

CM 

CM 

CM 

CO 

r» 

cO 

1 — 1 

CO 

rv 

CM 

CO 

0 

CO 

O 

>4 

S3 

i-l 

• 

• 

(2 

CO 

0 

0 

03 

44 

O 

0 

CO 

OC 

m 

p 

X- 

0 

O 

m 

m 

m 

CM 

m 

p 

m 

P 

cO 

<U 

TO 

Pi 

P 

< 23 

CO 

0) 

0 

• 

• 

• 

• 

• 

0 

OO 

CM 

ac 

CM 

• 

• 

• 

EH 

CX| 

r» 

CM 

cO 

KT 

t — ! 

cO 

0 

l 

0 

1 

m 

ac 

r\ 

CD 

rc 

co 

1 — 1 

r\ 

CM 

0 

cO 

0 

>4 

CO 

CO 

1 — 1 

m 

CM 

00 

cO 

a 

N 

•H 

03 

cO 

P 

S2 

• 

0 

S3 

X 

S3 

• 

pj 

• 

• 

Vw' 

• 

• 

;z 

0 

0  c 

b 

<3 

• 

• 

p 

fH 

• 

CO 

H 

p 

p 

• 

•  0 

•H 

03 

Q 

<3 

<3 

0 

CO 

<3 

<3 

w 

w 

pJ 

W 

u 

<d 

<3 

0 

0 

CO  *H 

44 

H 

p 

P-) 

s 

PQ 

M 

PQ 

Q 

0 

PQ 

S 

P 

p2 

w 

S3 

•  4-J 

03 

3 

12  cfl 

4-1 

S3 

4-> 

CO 

CO 

1  3-4 

CO  cu 

• 

U  4-) 

O 

rH 

CM 

CO 

m 

cO 

co 

OC 

0 

rH 

CM 

CO 

m 

cO  <U 
P-.  0 

S3 

rH 

P 

P 

rH 

p 

1 — 1 

TABLE  III- I  (CONT)  DRAINAGE  BASIN  CHARACTERISTICS 


110 


pi 

CM 

4J 

O 

O 

d 

o 

o 

o 

o 

o 

m 

ON 

o 

o 

r3 

CEJ  RJ 

o 

o 

m 

O' 

O' 

o 

o 

NO 

NO 

CM 

in 

m 

m 

o 

£ 

a 

1 — 1 

• 

CM 

o 

• 

• 

• 

o 

04 

co 

1 — 1 

m 

• 

• 

H 

Q 

<u 

* 

i — 1 

CO 

00 

NO 

o 

o 

#N 

r\ 

*N 

r\ 

CM 

CO 

o 

>-l 

ON 

4-J 

CM 

CM 

i 

1 

NO 

CM 

O' 

m 

CM 

NO 

CO 

T - 1 

• 

d 

Pd 

<U 

CM 

> 

O 

o 

m 

o 

m 

m 

ON 

O 

•H 

O' 

O' 

O' 

O' 

o 

ON 

O' 

O' 

o 

ON 

On 

NO 

O'* 

O 

r— 1 

cej  pd 

NO 

• 

CM 

CO 

• 

• 

• 

m 

CM 

rH 

o 

• 

• 

• 

Eh 

o 

t — 1 

* 

r — i 

CO 

o 

m 

o 

o 

r\ 

*\ 

#N 

O' 

t — 1 

o 

>m 

pq 

cu 

CO 

O' 

CM 

1 

m 

CM 

O' 

m 

CO 

NO 

ON 

PM 

CO 

o 

pc; 

d 

d 

<u 

1) 

t - 1 

> 

> 

o 

•H 

•H 

o 

o 

-  O 

o 

o 

m 

ON 

00 

o 

pd 

pc; 

o 

oo 

l — I 

O' 

O' 

o 

m 

O' 

O' 

NO 

m 

ON 

NO 

< 

CEJ 

00 

• 

ON 

NO 

• 

• 

• 

m 

CO 

i — 1 

CO 

• 

• 

• 

EH 

pq 

CO 

co 

rv 

1 — 1 

00 

O' 

1 — I 

o 

o 

ri 

r\ 

#N 

r\ 

m 

rH 

o 

ON 

CO 

CO 

Osl 

i — ! 

< — 1 

1 

1 

NO 

CM 

O' 

m 

CM 

NO 

o 

o 

Pd 

Pd 

pi 

d 

co 

rH 

o 

•  Pd 

o 

o 

O 

o 

o 

m 

00 

o 

£ 

d  CJ 

CO 

ON 

O' 

O' 

o 

O 

m 

m 

NO 

o 

CO 

m 

O' 

o 

iH 

d  RJ 

vO 

• 

oo 

• 

• 

• 

ON 

00 

CO 

m 

CM 

• 

• 

• 

H 

o 

RJ 

s 

r\ 

o 

vO 

O' 

o 

o 

#v 

r\ 

r\ 

r\ 

NO 

1 — 1 

o 

<! 

CO 

H 

04 

i — 1 

1 

CM 

t — 1 

O' 

m 

CM 

NO 

ON 

R) 

OO 

•  CJ 

pq 

• 

pi 

rH 

d 

o 

o 

o 

o 

o 

o 

00 

o 

d 

o 

oo 

CnI 

O' 

O' 

o 

NO 

CM 

CM 

m 

NO 

NO 

o 

•H 

H  rH 

• 

in 

CM 

• 

• 

• 

o 

CM 

1 — 1 

ON 

m 

• 

• 

• 

H 

w 

rH 

d  co 

r\ 

CO 

NO 

NO 

CM 

o 

o 

r\ 

r\ 

r. 

CO 

o 

o 

>-• 

< 

CO 

CU 

rH 

t - 1 

1 

in 

CM 

O' 

<3* 

CM 

NO 

ON 

0) 

H 

i — 1 

H 

• 

/*s 

>"S 

X 

(M 

• 

• 

• 

• 

• 

[z; 

O 

d 

<n 

• 

• 

t-q 

P-< 

C/D 

H 

H 

H 

o 

Q 

<J 

< 

u 

C/D 

<3 

< 

w 

W 

w 

W 

W 

a 

<d 

< 

O 

O 

CO  *H 

•H 

±j 

H 

d) 

PM 

a 

pq 

pq 

pq 

Q 

o 

a 

pp; 

pq 

a 

rJ 

W 

2 

•  -u 

ctf 

§ 

RJ 

4-J 

4-J 

CO 

CO 

1  u 

n3  0) 

• 

d  4-J 

o 

rH 

CM 

CO 

in 

NO 

O' 

CO 

ON 

o 

i — 1 

CM 

CO 

m 

RJ  a) 

1 — 1 

rH 

1 - 1 

1 — 1 

rH 

rH 

PM  0 

TABLE  III-I  (CONT)  DRAINAGE  BASIN  CHARACTERISTICS 


Ill 


pi 

P3 

CO 

CSl 

4-1 

m 

o 

4-1 

3 

O 

o 

o 

O 

o 

O 

CO 

o 

CO 

O 

oo 

CM 

CM 

O'. 

r-^ 

o 

o 

O 

o 

o 

CM 

in 

vO 

to 

<3 

r-H 

S 

CM 

• 

00 

VO 

• 

■ 

• 

o 

LO 

co 

00 

m 

• 

• 

• 

H 

W 

pin 

*v 

O 

rH 

vO 

o 

o 

cv 

r\ 

•\ 

CM 

i — 1 

O 

>4 

o 

1 - 1 

CO 

1 

<T 

rH 

00 

VO 

CM 

i — 1 

vO 

P3 

• 

i-P 

1 - 1 

CO 

(3 

o 

o 

LO 

o 

m 

m 

m 

O 

<D 

• 

O 

vo 

o 

rH 

o 

o 

VO 

00 

CM 

CM 

O'. 

vO 

o 

O 

CO 

LO 

• 

LO 

00 

• 

o 

• 

o 

CM 

in 

• 

« 

• 

H 

m 

P9 

c 

c 

4J 

CM 

r^- 

1 — 1 

1 — 1 

#v 

* 

n 

r\ 

rH 

o 

o 

>4 

<3 

cd 

cd 

1 

CM 

r^- 

LO 

CM 

VO 

o 

V4 

IS 

i — 1 

Ph 

• 

00 

3 

i — 1 

•H 

o 

o 

O 

O 

o 

O 

o 

03 

<Hl 

>< 

o 

o 

O'. 

00 

O'. 

o 

m 

CO 

o 

S-l 

LO 

LO 

vO 

1 — 1 

o- 

o- 

o 

O'. 

i"- 

r'- 

CM 

o 

vO 

m 

<j 

cd 

n 

• 

CO 

CT. 

• 

• 

• 

r. 

ev 

n 

rv 

• 

• 

H 

<3 

•H 

0) 

CM 

1 - 1 

1 - 1 

co 

o 

o 

VO 

rH 

in 

CM 

oo 

o 

o 

o 

a i 

rH 

1 

vO 

r— 1 

a 

pi 

• 

u 

03 

CO 

a) 

03 

cd 

o 

O 

o 

o 

o 

O 

o 

> 

1 — 1 

o 

o 

o 

o 

O 

CM 

o 

•H 

O 

N 

o 

o 

LO 

<t 

<r 

o 

o 

1 — 1 

CO 

Ov 

m 

o 

m 

LO 

<3 

03 

HO 

a 

* 

co 

CO 

• 

♦ 

o 

e\ 

r\ 

r\ 

r\ 

• 

• 

• 

H 

w 

(3 

03 

03 

CO 

1 — 1 

•<r 

1 — 1  • 

<3- 

1 — 1 

VO 

m 

1 - 1 

r— 

o 

>4 

c r. 

O 

C3 

1 

vO 

rH 

O 

PQ 

W 3 

• 

a 

P3 

•H 

CM 

o 

o 

o 

O 

o 

m 

O 

4-1 

o 

o 

IT) 

VO 

i — 1 

1 — 1 

CO 

o 

V4 

CH J 

4-> 

LO 

oc 

LO 

00 

LO 

o 

00 

CO 

in 

m 

OV 

CO 

in 

H 

1 - 1 

Q 

03 

u 

•  • 

oo 

? — 1 

» 

>4 

LO 

Q 

0) 

co 

o 

CO 

o 

o 

vO 

» — 1 

m 

CM 

co 

CO 

o 

Q 

•03 

tS 

1 — 1 

i 

1 

vO 

O'. 

4-1 

03 

CO 

4-1 

CO 

• 

o 

z 

X 

Po 

z 

a 

• 

• 

• 

• 

Z 

o 

co  c 

•  o 

<c 

• 

• 

hJ 

PH 

• 

• 

• 

• 

• 

• 

CO 

H 

H 

M 

o 

03 

Q 

<< 

<3 

CO 

CO 

<3 

<3 

w 

w 

pL] 

w 

M 

CO 

<3 

< 

CO 

o 

C/D  -H 

1J 

Eh 

hJ 

pH 

s 

PP 

PP 

PP 

Q 

o 

P3 

PP 

S 

rP 

p2 

•  -H 

cd 

Pm 

£2  cd 

4-1 

Z 

4-1 

CO 

CO 

1  U 

cd  <13 

• 

H  4-1 

o 

rH 

CM 

CO 

LO 

VO 

o- 

00 

cr> 

O 

rH 

CM 

CO 

in 

03  d) 
p4  a 

S3 

I — 1 

rH 

rH 

rH 

1 - 1 

rH 

TABLE  III-I  (CONT)  DRAINAGE  BASIN  CHARACTERISTICS 


112 


• 

• 

34 

rH 

pp; 

a) 

O 

o 

o 

o 

o 

o 

O 

o 

> 

o 

o 

m 

vO 

i — 1 

m 

a\ 

o 

i — 1 

o 

p3 

CM 

l - 1 

vO 

o 

m 

O' 

O' 

m 

CM 

O' 

m 

CM 

m 

00 

m 

H 

<d 

QJ 

rQ 

O 

• 

CM 

vO 

0) 

03 

o 

O 

1 - i 

CM 

o 

o 

m 

t — 1 

O' 

VO 

CM 

vO 

in 

o 

o 

PM 

a 

rH 

VO 

r— 1 

03 

C_> 

pi 

rH 

0) 

o 

o 

o 

o 

o 

00 

o 

a 

• 

qd 

co 

m 

m 

o 

CO 

O' 

O' 

o 

o 

o 

o 

o 

O' 

3 - 1 

in 

H 

o 

o 

u 

4-1 

o 

• 

CO 

• 

• 

• 

o 

o 

OV 

m 

o 

• 

• 

• 

pq 

4-1 

P3 

33 

rv 

o 

l — 1 

in 

o 

o 

#v 

#\ 

O' 

<T 

o 

Pd 

CO 

O 

VO 

O' 

00 

oo 

CM 

vO 

o 

03 

23 

t - 1 

<u 

Pd 

P^' 

03 

1 - 1 

rX 

dd 

O 

o 

o 

m 

m 

O 

03 

• 

4-1 

CO 

CM 

VO 

o 

o-n 

o 

o 

o 

o 

m 

vO 

CM 

vo 

VO 

H 

o 

rH 

Pi 

d> 

■  CO 

CM 

VO 

CO 

• 

• 

o 

in 

m 

m 

o 

in 

• 

• 

• 

pq 

£ 

P3 

O 

#s 

1 — 1 

CO 

1 — 1 

o- 

C\ 

1 - 1 

CM 

o 

O 

o 

23 

CO 

1 

CM 

CM 

VO 

o 

i — 1 

1 — 1 

i — 1 

•H 

rs 

pi 

3-i 

•H 

O 

o 

o 

o 

O 

o 

in 

t — 1 

C3 

a) 

o 

m 

m 

vO 

CO 

r~- 

O' 

o 

o 

o 

m 

m 

1 - 1 

CM 

vo 

H 

o 

a 

> 

a 

O'! 

• 

o- 

CO 

• 

• 

• 

o 

m 

rH 

CM 

o 

• 

• 

• 

o 

03 

o 

03 

rv 

o 

CM 

VO 

o 

o 

co 

rv 

r\ 

ft 

vO 

l — 1 

o 

LO 

o 

qq 

CO 

CM 

i 

VO 

a\ 

00 

CM 

VO 

w 

03 

03 

33 

I - 1 

o 

23 

1 — 1 

1 — 1 

03 

• 

O 

CO 

• 

CO 

Pd 

1 - 1 

1 - 1 

•H 

0) 

1 — 1 

a 

a) 

Pm 

O 

o 

o 

o 

o 

in 

CO 

o 

P3 

> 

uo 

3 - 1 

m 

00 

O' 

O' 

O' 

o 

o 

o 

o 

t — 1 

VO 

vo 

H 

o 

03 

o 

03 

vO 

• 

O 

• 

• 

• 

o 

vo 

1 — 1 

m 

vo 

• 

• 

• 

>H 

LO 

pq 

-C 

rO 

•H 

n 

1 - 1 

CM 

o 

o 

r\ 

rs 

#> 

r\ 

#\ 

CO 

1 — 1 

o 

w 

03 

03 

a 

in 

1 

vo 

1 — 1 

ov 

O' 

CM 

vo 

o 

23 

•H 

rH 

00 

• 

34 

CO 

•H 

> 

• 

o 

• 

/""S 

X 

Pn 

23 

p! 

• 

'w' 

• 

• 

2; 

o 

U  qj 

O 

<3 

• 

• 

►J 

IH 

CO 

H 

IH 

H 

• 

•  o 

•H 

CD 

Q 

<3 

<3 

CO 

CO 

<3 

<3 

w 

w 

pq 

pq 

w 

o 

<d 

<3 

O 

o 

CO  *H 
•  4-» 

U 

aj 

H 

t-J 

23 

pq 

pq 

pq 

Q 

e> 

s 

2 

pq 

23 

q 

p2 

s 

IS  a3 

4-1 

2: 

4-i 

CO 

CO 

1  3-1 

cfl  <U 

• 

34  4-> 

o 

rH 

CM 

CO 

mT 

m 

vo 

O' 

00 

CTv 

o 

rH 

CM 

CO 

m 

03  03 

23 

rH 

iH 

rH 

1 — 1 

rH 

rH 

P-l  0 

TABLE  III- I  (CONT)  DRAINAGE  BASIN  CHARACTERISTICS 


113 


65 

Moore 

Cr. 

8.1 

0.08 

61 

6*8 

3.3 

o 

• 

i — 1 

1 

0*0 

5,900 

3,100 

5,900 

2,800 

4,280 

420 

80‘  09 

9*0 

YT 

s 

o 

4_> 

• 

04 

oo 

o 

o 

o 

o 

o 

o 

3-1 

av 

04 

O'- 

VO 

oo 

o 

o 

m 

o 

m 

o 

p- 

m 

vO 

H 

vO 

CU 

CO 

« 

• 

OV 

• 

• 

• 

• 

CO 

m 

1 — 1 

n 

cr. 

• 

ex 

VO 

O 

04 

1 — 1 

rH 

o 

r\ 

r* 

rv 

rv 

04 

o 

o 

ex 

P 

i 

i 

04 

04 

CM 

v£> 

CO 

vO 

CM 

tn 

m 

rH 

04 

vO 

o 

o 

o 

o 

o 

O 

m 

CO- 

VO 

vO 

H 

CO 

P 

« 

• 

• 

CO 

• 

• 

• 

• 

o 

m 

o 

in 

04 

• 

• 

• 

vO 

o 

3-i 

VO 

o 

CT. 

I"'- 

o 

rH 

o 

04 

m 

CM 

00 

co 

o 

o 

3 

CJ 

vO 

\ - 1 

r\ 

r\ 

«s 

rv 

00 

vO 

P 

m 

04 

m 

CO 

04 

• 

3-i 

T3 

rH 

03 

C 

co 

> 

pq 

oo 

O 

•H 

O 

O 

00 

p- 

o 

O 

o 

o 

in 

CTv 

CO 

< 

pi 

(El  • 

CO 

o 

1 — 1 

• 

• 

• 

o 

CO 

o 

p~ 

CO 

• 

• 

• 

p-i 

04 

<3 

rH 

vO 

CM 

1 — 1 

o 

o 

m 

p» 

cr. 

T - 1 

<1* 

p^ 

co 

vO 

i — 1 

Pd 

4-1 

* 

04 

cv 

r\ 

ev 

rv 

r\ 

-  rH 

\ — 1 

a 

04 

oo 

04 

oo 

vO 

co 

•H 

4-1 

s 

a 

. 

■  H 

n 

S-i 

m 

03 

03 

o 

> 

> 

O 

o 

O 

o 

O 

m 

m 

o 

•H 

•H 

vO 

oo 

p- 

p- 

o 

04 

o 

00 

Ov 

O'. 

rH 

CO 

1 - 1 

<j 

Pi 

<21  Pi 

O 

o 

04 

CO 

• 

• 

• 

in 

Cv 

m 

00 

• 

• 

• 

PM 

vD 

<1 

r\ 

I - 1 

p^ 

o 

o 

r» 

r\ 

r> 

rv 

p- 

O'. 

rH 

Pd 

Pd 

1 — 1 

rH 

00 

CO 

oo 

m 

CO 

1 - 1 

l — 1 

1 - 1 

•H 

•H 

a 

a 

• 

o 

X 

to 

P 

a 

• 

• 

• 

s — ' 

• 

• 

O 

u  a 

o 

<3 

• 

• 

P 

CO 

H 

M 

M 

• 

•  o 

•H 

03 

P 

<1 

<3 

o 

CO 

c 

<3 

w 

w 

w 

w 

w 

u 

C 

<3 

o 

O 

cn  -h 

4-1 

H 

1-1 

pM 

a 

pq 

PQ 

pq 

p 

o 

3 

pq 

a 

P 

3 

w 

2 

•  4-J 

cti 

3 

s 

32  aj 

4-1 

p 

-U) 

CO 

co 

1  n 

a3  03 

• 

3-1  4-1 

o 

rH 

04 

co 

m 

vO 

p- 

oo 

cr. 

o 

rH 

CM 

CO 

in 

03  CL) 

rH 

rH 

1 - 1 

rH 

rH 

i — i 

P-<  0 

TABLE  III- I  (CONT)  DRAINAGE  BASIN  CHARACTERISTICS 


114 


C 

CM 

o 

o 

• 

CM 

r- 

CM 

1 - 1 

r^. 

r-~ 

O 

o 

CM 

o 

m 

t — 1 

rH 

vO 

H 

o 

CO 

3-i 

• 

CO 

1 - 1 

• 

• 

• 

• 

o 

1 - 1 

n- 

• 

• 

• 

>H 

r-~ 

4-1 

CJ 

O'- 

• 

av 

CO 

m 

o 

o 

cr. 

CM 

O'. 

VO 

m 

CO 

o 

o 

03 

CO 

1 — 1 

1 - 1 

1 

#\ 

#v 

r» 

in 

vO 

IS 

co 

CM 

CO 

rH 

CM 

34 

• 

O 

cr\ 

a) 

34 

in 

vO 

n- 

o 

O 

o 

m 

CO 

VO 

H 

VO 

CU 

CJ 

• 

o 

CTv 

• 

• 

• 

• 

o 

m 

o 

o 

00 

o 

• 

• 

>4 

PQ 

1 — 1 

O'. 

t — 1 

CM 

o 

o 

r-« 

m 

in 

o 

o 

o 

1 

c\ 

r\ 

CM 

CM 

vO 

CM 

CM 

CM 

CM 

O 

m 

Cv 

n- 

CM 

r-. 

O 

m 

O 

m 

o 

CO 

CO 

VO 

H 

00 

£ 

• 

vO 

CM 

co 

• 

• 

• 

• 

O 

CT. 

CO 

00 

m 

• 

• 

• 

VO 

o 

34 

« — 1 

« 

av 

00 

m 

o 

o 

1 - 1 

CO 

1 — 1 

IO 

m 

o 

o 

H 

a 

rH 

CM 

1 

i 

rv 

r\ 

r\ 

CM 

vO 

"3- 

CM 

m 

CM 

CM 

4-1 

o 

C 

• 

-4- 

o 

i — 1 

r-. 

o 

o 

O 

o 

o 

m 

<3* 

CM 

vO 

H 

u 

34 

• 

i — 1 

Ov 

• 

• 

• 

• 

o 

o 

o 

o 

CM 

o 

• 

• 

>4 

vO 

3 

CJ 

• 

o^ 

CO 

CO 

o 

o 

o 

1"- 

o 

CO 

i — 1 

CM 

o 

o 

PQ 

o 

1 

r» 

ev 

rv 

e\ 

vO 

. 

CM 

1 — 1 

co 

• 

3-i 

<u 

• 

o 

o 

rH 

vO 

o 

34 

• 

CO 

r-- 

• 

o 

O 

o 

o 

m 

1 - 1 

vO 

E-4 

vO 

£ 

CJ 

o 

• 

r-* 

VO 

• 

• 

• 

o 

o 

n- 

r-'. 

1 — 1 

O'! 

• 

• 

>4 

vO 

<U 

VO 

o 

1 - 1 

o 

o 

CO 

00 

vO 

00 

VO. 

Cv 

o 

O 

a 

i 

rv 

rv 

c\ 

#v 

VO 

CO 

vO 

CM 

vO 

CO 

CO 

• 

o 

33 

X 

3o 

<3 

a 

• 

• 

• 

• 

;z 

O 

cj  a 

o 

< 

• 

• 

rJ 

fH 

CO 

H 

H 

H 

• 

•  O 

•H 

CU 

Q 

C 

<i 

CJ 

CO 

C 

C 

w 

w 

M 

w 

w 

CJ 

<d 

<3 

CJ 

o 

CO  *H 

4-1 

pi 

H 

rJ 

pH 

S 

PQ 

PQ 

PQ 

Q 

CJ 

2 

PQ 

S 

2 

W 

S3 

•  U 

cd 

2 

33  03 

4-1 

!3 

4-> 

CO 

CO 

1  34 

as  cu 

• 

34  4-1 

O 

rH 

CM 

CO 

m 

VO 

f" 

00 

CTv 

o 

i — 1 

CM 

CO 

<r 

m 

03  0) 

33 

1 — 1 

rH 

rH 

1 — 1 

rH 

1 — 1 

PQ  £ 

TABLE  III-I  (CONT)  DRAINAGE  BASIN  CHARACTERISTICS 


APPENDIX  IV 


MAP  COEFFICIENTS  FOR 


R(av)  ,  K(av) ,  BF  and  PEAK  FLOW 
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FIG.  IV— 2 


MAP  COEFFICIENTS  FOR  R  (av)  FOR  USE  WITH  EQUATION  4-|| 
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FIG.  IV— 3 


MAP  COEFFS  FOR  K  (av)  FOR  USE  WITH  EQUATION  4-12 
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FIG.  IV-4 


MAP  COEFFICIENTS  FOR  K  (av)  FOR  USE  WITH  EQUATION  4*12 
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FIG.  IV-5  MAP  COEFFS  FOR  BASEFLOW  (BF)  FOR  USE  WITH  EQUATION  4' 13 
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FIG.  IV —6 


HAP  COEFFICIENTS  FOR  BASEFLOW  (BF)  FOR  USE  WITH  EQUATION  4- 13 
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FIG.  IV-7 


MAP  COEFFICIENTS  FOR  PEAK  FLOW  (Qo)  FOR  USE  WITH  EQUATIONS 
5-10  TO  5-15  (FOR  0-1  TO  10  INS  OF  EXCESS  PRECIPITATION) 
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FIG.  IV  —8 


MAP  COEFFICIENTS  FOR  PEAK  FLOW  (Qo)  FOR  USE  WITH  EQUATIONS 
5-IO  TO  5-15  (0-1  -  10  INS  OF  EXCESS  PRECIPITATION) 


